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ABSTRACT 


Data  from  the  field  and  laboratory,  dealing  with 
the  behaviour  of  coarse-bed  rivers,  are  analyzed  under  the 
following  topics:  the  sampling  and  analysis  of  coarse  bed- 
material,  the  threshold  of  motion  condition,  the  amount  of 
material  in  transport,  the  resistance  to  flow,  the  design 
of  stable  channels,  and  the  depth  of  .scour  at  river  bends. 

The  adequate  sampling  of  a  coarse  river-bed  should 
consist  of  a  subsurface  scoop  sample,  borings  and  a  surface 
grid  sample.  Analysis  of  a  wide  range  of  field  data  in¬ 
dicated  that  the  largest  stones  in  a  mixture  start  to  move 
at  lower  velocities  than  the  same  size  stones  in  a  uniform 
bed  of  rip  rap.  The  analysis  of  flume  experiments  using 
dimensional  analysis  resulted  in  the  conclusion  that  the 
relative  depth  is  a  significant  parameter  in  assessing  the 
rate  of  bed-load  transport.  A  comparison  of  the  measured 
bed-load  transport  in  the  Elbow  and  North  Saskatchewan 
Rivers  to  computed  transport  indicated  that  the  existing 
formulas  are  inadequate.  The  measurement  of  the  actual 
protrusion  height  of  stones  on  a  river-bed  resulted  in  the 
development  of  a  flow  resistance  equation,  for  rigid  beds, 
which  is  somewhat  dissimilar  from  the  formulas  of  Keulegan 
and  Kellerhals .  Formulas  relating  stable  channel  width, 
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depth,  and  slope  to  discharge  and  size  of  bed-material  were 
developed  for  coarse-bed  channels.  The  sounding  of  river 
bends  before,  during,  and  after  floods  revealed  that  scour 
holes  at  bends  were  relatively  stable;  filling  of  the  scour 
hole  during  recession  of  the  high  flows  did  not  take  place. 
Analysis  of  a  number  of  forced  and  free  bends  led  to  the 
development  of  tentative  design  formulas  for  assessing  the 
maximum  depth  of  scour.  This  maximum  depth  of  scour  was 
related  to  the  width,  radius  and  internal  angle  of  the 
river  bend. 

The  importance  of  more  adequate  field  data  was 
stressed  and  recommendations  regarding  field  and  laboratory 
investigations  were  formulated. 
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CHAPTER  I 


INTRODUCTION 

Engineers  today  are  being  confronted  with,  in¬ 
creasingly  complex  problems  in  river  engineering.  In 
Particular /  the  training  of  rivers  in  mountainous  areas 
has  called  for  increased  attention  because  of  the  con¬ 
struction  of  vast  hydro-electric  and  diversion  schemes, 
pipelines,  and  recreation  facilities.  Unfortunately, 
most  of  the  river  training  experience  that  engineers  have 
accumulated  is  related  to  sand-bed  rivers  and  is  inadequate 
for  dealing  with  coarse-bed  mountainous  rivers.  The  term 
coarse-bed  river  refers  to  a  river  having  an  immobile  or 
a  mobile-bed  of  non-cohesive  material  larger  than  2  mm.  in 
size.  This  bed-material  is  generally  classified  into 
gravel  (2  mm.  to  2.5  inches ), cobbles  (2.5  inches  to  10 
inches) ,  and  boulders  (over  10  inches) . 

Although  the  behaviour  of  coarse-bed  rivers  was 
discussed  in  the  nineteenth  century  by  DuBoys  (1879) ,  few 
field  investigations  have  been  conducted  to  date.  An  in¬ 
vestigation  by  Lane  and  Carlson  (1953,  1954)  yielded  use¬ 
ful  information  on  the  design  of  canals  in  this  material. 
Through  a  cooperative  effort  by  the  University  of  Alberta, 
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the  Alberta  Department  of  Highways,  the  Water  Resources 
Division  and  the  Alberta  Research  Council,  investigations 
were  carried  out  to  assess  the  character  of  rivers  in  Alberta 
and  British  Columbia  (Qureshi  1962,  Kellerhals  1963,  and 
Van  Der  Giessen  1966) .  The  study  by  Kellerhals  extended 
the  knowledge  in  regard  to  stable  channel  design  and  re¬ 
sulted  in  new  design  equations.  Subsequent  investigations 
by  Hollingshead  (1968a,  1968b)  yielded  information  on  bed¬ 
load  discharge.  Investigations  on  sixteen  rivers  in 
Pennsylvania  by  Brush  (1961)  resulted  in  information  on 
the  relation  of  the  river  channel  to  the  geologic  character¬ 
istics  of  '  the  basin.  Fahnestock  (1963)  studied  unstable 
coarse-bed  streams  flowing  from  glaciers  and  obtained  in¬ 
formation  on  the  shifting  of  channels  and  the  stability 
of  the  bed-material.  Investigations  on  rivers  in  Europe, 
in  particular  by  Meyer-Peter  (1949)  and  by  Ramette  and 
Heuzel  (1962) ,  were  concerned  primarily  with  the  applic¬ 
ability  of  the  Meyer-Peter  and  Muller  bed-load  formula. 

As  part  of  the  continuing  river  research  program 
in  Alberta  the  North  Saskatchewan  River  at  Drayton  Valley 
was  investigated  during  the  summers  of  1965  and  1966 
(Galay  1967a) .  Information  regarding  the  shifting  of  the 
main  channel  and  the  movement  of  large  bed-material  during 
near-bankfull  flows  was  obtained.  Subsequent  longitudinal 
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soundings  on  the  North  Saskatchewan  River  downstream  of 
Edmonton,  the  Oldman  River  and  the  Athabasca  River  resulted 
in  valuable  data  concerning  scour  at  river  bends  (Hollings- 
head  and  Schultz,  1968) .  The  location  of  the  investigations 
is  shown  in  FIGURE  1.  A  limited  field  program  was  carried 
out  on  Wilson  Creek  in  Manitoba  in  1968  to  assess  the  flow 
conditions  necessary  to  move  coarse  bed-material. 

The  Alberta  information  combined  with  data  from 
other  investigations  will  be  analyzed  in  depth  to  yield 
formulas  that  may  be  of  assistance  to  engineers.  Specific 
problems  receiving  attention  are: 

(1)  What  bed  sampling  techniques  are  best  for  immobile 
or  mobile-bed  rivers? 

(2)  Given  the  size  distribution  of  the  bed- material , 
under  what  flow  conditions  will  it  begin  to  move? 

(3)  What  is  the  best  method  for  estimating  the  amount 
of  bed-material  in  transport  under  known  flow 
conditions? 

(4)  What  is  the  resistance  to  flow  for  an  immobile  or 
a  mobile-bed  of  coarse  material? 

(5)  Given  the  expected  flow  conditions  and  the  size 
distribution  of  the  bed-material,  what  width. 
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depth  and  slope  will  result  in  a  stable  channel? 

(6)  What  is  the  maximum  depth  of  scour  that  can  be 
expected  at  a  river  bend? 

It  is  emphasized  that  the  lack  of  data  prevents  the 
answers  to  these  problems  from  being  conclusive  and  compels 
the  exclusion  of  the  consideration  of  several  important 
problems,  for  example,  the  design  of  training  works.  After 
the  analysis  of  existing  data,  specific  recommendations  are 
presented  to  guide  future  investigations. 
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CHAPTER  II 


SAMPLING  AND  ANALYSIS  OF  BED-MATERIAL 

2.1  INTRODUCTION 

The  purpose  of  sampling  and  subsequently  analyzing 
the  bed-material  in  coarse-bed  rivers  is  three-fold: 

(1)  To  assess  the  size  of  bed-material 
that  will  just  begin  to  move  under 
known  flow  conditions ,  (threshold 

of  motion)  ,  as  well  as  ultimate  degradation; 

(2)  To  estimate  the  volume  of  material 
transported  along  the  bed  under 
known  flow  conditions; 

(3)  To  assess  the  roughness  or  the  resis¬ 
tance  due  to  friction  that 

(a)  an  immobile  bed  offers  to  the  flow, 

(b)  a  mobile  bed  offers  to  the  flow. 

In  all  the  above  cases  it  is  necessary  to  know  a 
representative  size  as  well  as  the  variation  in  size  en¬ 
countered  on  the  river-bed.  This  chapter  will  discuss  the 
basic  types  of  bed-material  samples  that  can  be  obtained, 
the  various  sampling  techniques  with  their  appropriate 
methods  of  analysis,  and  the  relation  between  the  purpose 
of  sampling  and  the  available  sampling  techniques. 

2.2  TYPES  OF  BED-MATERIAL  SAMPLES 

A  review  of  the  literature  on  coarse-bed  rivers 
reveals  that  there  are  many  techniques  used  in  sampling 
the  river-bed,  and  that  the  various  sampling  techniques 
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yield  three  basic  types  of  samples: 

(1)  Volumetric  (or  bulk) 
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(2)  Areal  (or  surface) 

(3)  Combination  of  volumetric  and  areal. 

The  basic  types  of  samples  along  with  the  associat¬ 
ed  sampling  techniques  are  summarized  in  TABLE  1.  These 
techniques  will  now  be  discussed  in  more  detail. 

2.3  VOLUMETRIC  SAMPLING  TECHNIQUES 
2.3(a)  Scoop  Sampling  Technique 

This  technique  is  rather  straightforward  with  an 
appropriate  sampler  being  used  to  scoop  up  material  from 
the  bed  of  a  river  (Inter-Agency  Committee  on  Sedimentation 
1963) .  Material  can  also  be  shovelled  from  open  bars  into 
containers.  This  results  in  a  sample  having  a  mixture  of 
surface  as  well  as  sub-surface  material;  this  technique 
can,  however,  be  refined  so  that  only  surface  or  sub-surface 
material  is  picked  up.  PHOTOGRAPH  1  shows  a  sub-surface 
scoop  sample  being  obtained  at  6  inches  below  the  surface 
of  an  open  bar. 

2.3(b)  Square-Depth  Sampling  Technique 

In  this  case  all  the  material  to  a  depth  equal  to 
the  largest  stone  within  a  small  enclosed  area  on  an  open 
bar  is  collected.  PHOTOGRAPH  2  shows  a  three  foot  square 
framework  prior  to  the  removal  of  the  enclosed  material. 
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SUMMARY  OF  BED-MATERIAL  SAMPLING  TECHNIQUES  AND  ANALYSIS 
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This  technique  also  results  in  a  sample  having  a  mixture 
of  surface  and  subsurface  material  (Van  Der  Giessen  1966, 

p.  37)  . 

2.4  AREAL  SAMPLING  TECHNIQUES 
2.4(a)  Grid  Sampling  Technique 


This  sampling  technique  was  initially  devised  by 

Wolman  (1954)  and  quoting  directly  consists  of: 

In  the  desired  reach  of  the  stream  a  grid 
system  is  established  either  by  pacing  or 
with  actual  lines.  The  size  of  the  grid 
is  determined  by  the  length  of  reach  which 
the  sampler  desires  to  describe.  ...  After 
establishing  the  grid,  ...  100  individual 
pebbles  are  picked  up  from  the  bed  (at 
designated  grid  points) .  The  sampling  is 
probably  less  subjective  when  lines  or 
tapes  rather  than  pacing  are  used  to  fix 
the  individual  sampling  points.  Randomness 
in  the  selection  of  each  pebble  can  only  be 
obtained  if  the  sampler  tries  not  to  look 
at  the  bed  as  he  picks  up  each  pebble..  The 
author's  practice  is  to  draw  each  pebble 
from  beneath  the  tip  of  the  toe  of  his  boot. 

This  sampling  technique  has  been  widely  used  in 

North  America,  in  some  cases  with  slight  modifications. 

Brush  (1961)  noted  that  60  pebbles  were  adequate  for  a 

sample.  Various  other  investigators  using  this  technique 

were  Qureshi  (1962)  on  the  Red  Deer  River,  Kellerhals 

(1963)  on  rivers  in  British  Columbia,  Fahnestock  (1963) 

on  the  White  River  and  Van  Der  Giessen  (1966)  on  streams 

in  Western  Alberta.  PHOTOGRAPH  3  illustrates  the  pacing 

technique  in  grid  sampling. 


■ 
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Several  investigators  (Kellerhals  1963)  have 
noted  that  the  sample  obtained  by  pacing  is  slightly 
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different  from  that  obtained  by  using  a  grid  with 
established  lines.  This  discrepancy  is  due  to  two 
factors : 

(1)  The  operator's  toe  usually  falls  on 
a  stone  projecting  from  the  bed  in¬ 
stead  of  in  the  space  between  the 
stones  resulting  in  a  larger  stone 
being  picked  up; 

(2)  In  picking  up  a  stone  beneath  an 
operator's  toe  the  finger  usually 
comes  in  contact  with  the  larger 
stones  before  the  smaller  ones, 
even  thbugh  these  larger  stones 
may  not  be  directly  below  the  toe. 

2.4(b)  Line  Sampling  Technique 

This  .  technique  is  essentially  a  minor  modificat¬ 
ion  of  the  grid  sampling  technique  in  that  pebbles  are 
picked  up  at  regular  intervals  along  an  outstretched 
tape  instead  of  at  grid  points .  The  sample  should  be 
identical  to  that  obtained  by  the  grid  sampling  technique 
and  is  well  suited  for  sampling  along  river  banks  or  on 
top  of  long  narrow  bars  (PHOTOGRAPH  4) . 

2.4(c)  Square-Surface  Sampling  Technique 

A  square-surface  sample  is  obtained  by  picking 
up  all  the  surface  pebbles  within  a  small  enclosed  area 
on  an  open  bar.  This  area, is  taken  to  be  representative 
of  the  whole  channel  bed.  This  technique  was  initially 


■ 
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used  by  Lane  and  Carlson  (1953)  in  their  studies  of  the 
San  Luis  Valley  Canals. 

A  variation  of  this  technique,  which  was  used  on 
the  North  Saskatchewan  River  (Galay  1967a) ,  is  to  photo¬ 
graph  the  stones  in  an  enclosed  area  with  some  scale  in  the 
photograph  (PHOTOGRAPH  5) .  The  stones  usually  orient  them¬ 
selves  so  that  the  long  axis  is  perpendicular  to  the  flow 
and  the  intermediate  axis  parallel.  These  two  axes  can  be 
measured  with  sufficient  accuracy  directly  from  the  photo¬ 
graph.  This  photographic  technique  has  also  been  used  by 
Pashinskiy  (1964) ,  Church  (1968) ,  and  Ritter  and  Helley 
(1969)  . 

2.5  COMBINATION  SAMPLING  TECHNIQUES 

The  sampling  techniques  in  this  section  could  be 
combinations  of  the  volumetric  and  areal  techniques 
previously  discussed.  There  are  several  possible  combinat¬ 
ions  but  only  that  developed  by  Wolman  (1955)  will  be  dis¬ 
cussed.  This  combination  technique  consists  of  obtaining 
four  scoop  samples  of  material  finer  than  2  inches  along  a 
transect  (line  across  a  channel)  with  all  the  pebbles  larger 
than  2  inches  being  measured  directly. 

2.6  ANALYSIS  OF  BED-MATERIAL  SAMPLES 

A  sample  of  river  bed-material  obtained  by  one  or 
several  of  the  previously  mentioned  techniques  has  to  be 
analyzed  and  described  in  a  reasonable  numerical  manner . 


■ 
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Engineers  are  generally  interested  in  a  representative 
size  of  the  sample  as  well  as  its  dispersion  or  sorting. 
2.6(a)  Volumetric  Samples 

A  scoop  sample,  a  square-depth  sample  and  the 
finer  portion  of  Wolman ' s  (1955)  combination  sample  are 
normally  sieved  and  weighed  in  order  to  arrive  at  a 
distribution  of  sizes.  The  results  from  a  sieve  analysis 
are  generally  plotted  on  a  graph  with  cumulative  %  finer 
(ordinate)  versus  the  sieve  size  (abscissa)  as  shown  in 

FIGURE  2.  In  some  cases  it  is  more  convenient  to  plot 

* 

the  sieve  results  on  logarithmic-probability  paper  as 
shown  in  FIGURE  3.  Natural  sands  often  plot  as  straight 
lines  on  this  type  of  paper  (Blench  1952,  Vanoni  et  al 
1960) ,  but  this  is  not  true  for  all  coarse  materials  as 
they  may  be  bimodal.  The  log-probability  plot  facilitates 

the  determination  of  percentile  values. 

* 

2.6(b)  Areal  Samples 

As  previously  mentioned  a  grid  sample  should  be 
identical  to  a  line  sample  and  for  the  forthcoming  dis¬ 
cussion  only  a  grid  sample  will  be  mentioned.  Also,  a 

grid  sample  obtained  by  pacing  will  usually  yield 
pebbles  having  a  larger  size  than  a  sample  obtained  from 
established  grid  points  (FIGURE  4).  Therefore,  wherever 
possible,  a  grid  sample  should  be  obtained  from  an 
established  grid  instead  of  by  pacing. 
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It  should  also  be  noted  that  a  grid  sample 
analyzed  by  number  will  yield  an  entirely  different  curve 
than  an  analysis  of  the  same  sample  by  weight  (FIGURE  5) . 
The  analysis  by  weight  is  greatly  influenced  by  any  large 
cobbles  in  the  sample.  There  does  not  appear  to  be  any 
good  reason  why  a  grid  sample,  which  is  areal  and  there¬ 
fore  two-dimensional,  should  be  analyzed  by  weight  which 
is  based  on  volume  and  therefore  three-dimensional. 
Although  the  grid  sample  has  been  analyzed  by  weight  by 
several  investigators  it  is  recommended  that  this 
analysis  procedure  be  discarded  for  grid  samples. 

•The  square-surface  sample  can  be  analyzed  by 
number  (taking  into  account  all  the  stones  within  the 
enclosed  area) ,  by  weight  and  by  percent  of  the  bed  area 
(from  a  photograph)  covered  by  a  certain  size  or  smaller 
The  resulting  distribution  curves  may  be  quite  different 
from  each  other.  This  may  be  better  understood  by  con¬ 
sidering  a  hypothetical  square-surface  sample  of  100 
spheres  of  various  sizes  as  shown  in  TABLE  2.  The  sizes 
to  be  considered  range  from  0.5  inches  to  5.0  inches  and 
there  are  10  spheres  of  each  size  in  the  sample.  This 
square-surface  sample  is  analyzed  in  the  three  different 
ways  mentioned  previously  and  plotted  in  FIGURE  6.  The 
plotted  points  are  joined  by  smooth  lines  instead  of 
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ANALYSIS  OF  SQUARE-SURFACE  SAMPLE  SPHERES 
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broken  lines  in  order  to  appear  similar  to  conventional 
distribution  curves.  The  analysis  by  number  yields  re¬ 
presentative  diameters  smaller  than  would  be  obtained  from 
the  weight  or  percent  of  bed  area  analysis.  The  latter  two 
analysis  techniques  yield  distribution  curves  that  are 
similar. 

In  grid  sampling,  which  is  most  common,  the  pro¬ 
bability  of  picking  up  a  certain  stone  is  dependent  upon  its 
surface  area  and  it  should  be  possible  to  demonstrate  that 
a  grid  sample  analyzed  by  number  is  similar  to  a  square- 
surface  sample  analyzed  by  percent  of  bed  area  covered. 

FIGURE  7  illustrates  that  a  grid  sample  analyzed  by  number 
does  yield  a  similar  curve  to  a  photographic  square-surface 
sample  analyzed  by  percent  of  bed  area.  It  is,  therefore, 
not  surprising  that  Kellerhals  (1967)  found  that  a  square- 
surface  sample  analyzed  by  weight  yields  very  similar  distri¬ 
bution  curves  to  a  grid  sample  analyzed  by  number.  Ritter 
and  Helley  (1969) ,  using  a  size  particle-size  analyzer,  came 
to  a  similar  conclusion.  A  paced  grid  sample,  from  the 
Chowchilla  River,  analyzed  by  number  agreed  closely  with  a 
square-surface  sample  (from  photograph)  analyzed  by  volume 
(FIGURE  8) .  The  volumetric  analysis  is  equivalent 
to  the  weight  analysis  if  all  the  particles  have  identical 
specific  weights.  The  volume  of  a  particle,  which  has  its 
intermediate  axis  determined  from  a  photograph,  can  be  cal¬ 
culated  by  multiplying  the  cube  of  the  intermediate  axis  by 


. 
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a  particle  shape  coefficient  (Ritter  and  Helley,  1969,  p.6). 

It  would  therefore  be  expected  that  a  photographic  sample 

(square-surface)  analyzed  by  volume  should  agree  closely  to 

a  sieve  analysis  of  the  same  surface  particles.  Ritter  and 

Helley  (1969),  however,  state: 

For  each  of  the  three  samples  (two  of 
coarse  sand,  one  of  ash  slough)  the 
median  diameters  determined  by  sieve 
analysis  are  smaller-almost  half  in 
two  samples  -  than  those  determined 
by  the  particle-size  analyzer. 

This  statement  may  be  in  error  as  the  particles  that  were 

analyzed  may  not  have  been  identical;  the  sieve  analysis 

may  have  been  carried  out  on  a  "depth"  sample  while  the 

volume  analysis  with  the  particle-size  analyzer  would  be  on 

a  surface  sample.  These  two  samples  could  be  distinctly 

different  from  each  other.  Ritter  and  Helley  (1969)  do  not 

clearly  describe  how  the  sieve  samples  were  obtained. 

Also,  the  grid  sample  analyzed  by  number  appears 

to  yield  the  most  useful  distribution  curve  as  it  conveys 

a  picture  of  river's  bed.  One  can  visualize  the  meaning 

of  the  following  statement:  "80%  of  the  river's  bed  is 

covered  by  pebbles  equal  to  or  smaller  than  2  inches". 
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2 . 6  (c )  Combination  Samples 

Wolman's  (1955)  combination  sample,  as  previously 
discussed  does  not  give  a  definite  picture  of  the  surface 
or  sub-surface  of  a  river-bed  and  will  not  be  dealt  with 
further. 

2.7  IMPORTANT  BED-MATERIAL  PARAMETERS 

In  studying  the  behaviour  of  the  river-bed  engineers 

are  generally  interested  in  a  representative  size  as  well  as 

the  dispersion  of  the  material  that  makes  up  the  bed. 

2.7(a)  Representative  Size 

There  are  several  measures  used  in  arriving  at  a 

representative  size  of  bed-material,  the  most  common  being 

the  median  or  Dr.  size.  The  median  size  is  the  50%  size  on 

o  u 

a  cumulative  frequency  curve  and  defines  the  size  separating 

the  sample  into  two  equal  halves.  The  median  size  is  the 

easiest  to  determine  but  since  it  is  not  influenced  by  the 
* 

distribution  of  particle  sizes  in  a  sample  its  use  as  a 
measure  of  representative  size  is  not  recommended  (Folk  1965, 
1966). 

Another  common  measure  is  the  arithmetic  mean  which 
is  obtained  by  adding  up  all  the  material  sizes  and  dividing 
by  the  number  of  values.  This  measure  is,  however,  difficult 
to  derive  for  material  finer  than  1/2  inch  and  a  number  of 
investigators  resort  to  the  arithmetic  mean  by  weight  which 
is  equal'  to  £ P where  Pi  =  %  in  class  interval  of  mean  diameter 
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Di.  This  latter  measure  is  also  referred  to  as  the  arith¬ 
metic  mean  computed  by  the  method  of  moments.  The  deter¬ 
mination  of  the  arithmetic  mean,  by  any  technique,  is, 
however,  time  consuming. 

Folk  (1965)  recommends  the  use  of  a  graphic  mean 
diameter  which  is  based  on  geologic  phi  units.  In  using 
percentile  values  directly  from  log-probability  graphs  this 
measure  would  become  the  geometric  mean  and  could  be  deter¬ 
mined  by: 


°G  =  V  D84  X  °50  X  °16  <2-X> 

Folk  (1965)  states  that  this  measure  corresponds  very 
closely  to  the  mean  as  computed  from  the  method  of  moments, 
yet  it  is  much  easier  to  use.  However,  the  study  of  dist¬ 
ribution  curves  for  coarse  materials  has  shown  that  the 
geometric  mean  is  equal  to  the  median  for  samples  which 
are  log-normal  (Galay  1967a) .  This  can  be  shown  by  consider¬ 
ing  the  distribution  curve  for  sample  A  in  FIGURE  9.  Using 
equation  2.1  the  geometric  mean  is: 


Dq  =  3J  2.69  x  l.Ox  0.37  =  1.0  inches 

which  is  identical  to  the  median  diameter.  The  arithmetic 
mean  is,  however,  equal  to  1.55  inches.  A  second  distri¬ 
bution  curve  for  sample  B  is  shown  in  FIGURE  9  to  illustrate 
that  both  the  median  and  the  geometric  mean,  being  equal  to 
1.0  inches  for  both  samples,  are  poor  measures  of  representative 


. 

4 


18 

size.  Howeverr  the  geometric  mean  does  approximate  the 
arithmetic  mean  fairly  closely  for  bimodal  distributions 
(Folk  1965) . 

A  convenient  measure  of  representative  particle  size 
would  be  (Galay  1967a) : 

D  =  10%  +  30%  +  50%  +  70%  +  90%  sizes  (2.2) 
m  - c - 


The  representative  size  as  obtained  from  the  above  relation 
is  reasonably  close  to  the  arithmetic  mean  as  shown  by 


comparing  the  values  for  the  two  samples  in  FIGURE  9 : 


Sample 


Arithmetic 

Mean  Dm  D50 

(in.)  (in.)  (in.) 


DG 

(in. ) 


A 


1.55  1.45  1.00 


1.00 


B 


1.03 


1.02  1.00  1.00 


More  investigation  into  representative  particle 


size  is  necessary  in  order  to  make  final  recommendations. 
Presently,  however,  most  of  the  data  as  reported  by  various 
investigators  lists  the  median  as  the  representative  par¬ 
ticle  size;  it  is  this  value  that  will  be  used  in  sub¬ 
sequent  analysis. 

The  mode,  which  is  the  most  frequently  occurring 
particle  size,  is  rarely  used  by  engineers  although  it  is 
useful  in  assessing  the  source  of  moving  sediment  in  natural 


rivers . 


’ 
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2.7(b)  Dispersion 


Trask's  sorting  coefficient  =  v/D75//d  ,  using 

25 

millimeter  sizes,  is  the  most  common  engineering  measure 
of  sorting  or  dispersion  but  it  measures  only  the  sorting 
in  the  central  part  of  the  distribution  curve.  A  better 
measure  as  noted  by  Simons  and  Richardson  (1966)  would  be 
the  standard  deviation  or  gradation  which  is  given  by  the 


formula : 


a,  = 


,  ,D50  ,  D84. 
h  ltttf  +  ftettJ 


D16 


D50 


(2.3) 


This  measure  would  be  especially  applicable  to 
log-normal  distributions  and  would  still  be  a  reasonable 
measure  of  dispersion  for  distributions  that  are  not  log¬ 
normal  . 


2.8  LOCATION  OF  SAMPLES 

The  knowledge  of  the  variation  in  sizes  of  bed- 
material  with  location  on  the  river-bed  is  of  importance. 
The  pool  and  riffle  sequence  in  coarse-bed  rivers  has  been 
well  documented  (Leopold,  Wolman  and  Miller  1964)  with  the 
bed  material  in  riffles  being  larger  than  that  of  pools. 

It  would  be  advantageous  to  obtain  samples  from  open  bars, 
from  pools  and  from  riffles  and  average  the  results  from 
these  three  samples  in  order  to  arrive  at  a  representative 
picture  of  the  river-bed.  Sampling  the  river-bed  in  pools, 
especially  during  high  flows,  may,  however,  be  somewhat 


hazardous . 


> 
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It  is  interesting  to  note  that  the  distribution  of 
material  in  transport  as  collected  by  a  bed-load  sampler 
resembles  the  distribution  of  a  sample  obtained 

by  the  square-depth  sampling  technique  (Hollingshead  1968a) . 

A  comparison  of  the  two  distribution  curves  is  shown  in 
FIGURE  10  .  However ,  more  samples  of  various  types  from  a 
variety  of  coarse-bed  rivers  would  be  necessary  in  order  to 
make  conclusive  statements. 

2.9  SAMPLING  AND  ANALYSIS  TECHNIQUES  IN  RELATION 
TO  THE  PURPOSE  OF  SAMPLING 

As  outlined  in  the  introduction  the  purpose  of 
sampling  the  material  on  the  bed  of  a  river  is  three-fold 
and  each  purpose  will  now  be  discussed  in  more  detail. 

2.9(a)  Threshold  of  Motion 

In  order  to  assess  what  particles  will  begin  to 
move  along  the  river-bed  under  known  flow  conditions  it  is 
essential  to  have  a  distribution  curve  of  the  particles 
making  up  the  surface  layer.  In  this  instance  the  sample 
could  be  obtained  by  a  grid  technique  and  analyzed  by  number. 
Sampling  along  river  or  channel  banks  is  also  of  importance 
in  assessing  bank  stability.  A  grid  sample  analyzed  by 
number  should  suffice  here. 

Ultimate  degradation  of  a  river  is  achieved  when 
the  small  bed  particles  have  been  removed  leaving  a  layer  of 
particles  that  are  too  large  to  move.  The  remaining  particles 
would  be  very  close  to  the  threshold  of  motion  state. 


i  . 


. 
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The  sampling  of  the  river-bed,  however,  to  assess  degradation, 
should  take  place  at  various  depths  below  the  bed  surface; 
core  samples  analyzed  by  weight  would  be  necessary  here. 
2.9(b)  Sediment  Transport 

Once  the  surface  layer  of  particles  is  in  motion  the 
sub-surface  particles  will  be  carried  off  in  a  rush  and  the 
moving  material  could  range  from  fine  sand  to  boulders.  A 
distribution  curve  of  sub-surface  material  would  be  necessary 
in  order  to  assess  the  volume  of  material  in  motion  when 
using  a  sediment  transport  formula.  This  material  could  be 
obtained  by  a  scoop  sample  at  some  depth  below  the  surface 
and  would  be  analyzed  by  weight.  A  square-depth  sample 
(from  an  open  bar)  analyzed  by  weight  should  also  yield  a 
representative  picture  of  the  moving  bed  material. 

2.9(c)  Resistance  to  Flow 

It  is  generally  known  that  the  resistance  to  flow 
of  an  alluvial  river  bed  depends  on  whether  a  bed  is  immobile 
or  mobile.  Each  of  these  aspects  will  now  be  discussed. 

2.9(c)i  Immobile  River-bed 

In  this  case  only  the  surface  material  is  relevant 
and  a  grid  sampling  technique  with  subsequent  analysis  by 
number  could  be  used.  It  is  important  here  that  a  clear 
picture  of  the  make-up  of  the  bed  surface  be  given  in  order 
that  a  representative  size  could  be  utilized  in  a  flow 


formula. 


. 
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2.9(c)ii  Mobile  River-bed 

The  resistance  to  flow  may  depend  on  the  amount  of 
material  moving  along  the  bed  if  bed  forms  are  generated 
and  may  vary  from  the  case  with  no  moving  particles.  Sub¬ 
surface  material  would  be  moving  along  with  surface  material 
and  sampling  should  be  by  scoop  with  analysis  by  weight. 

2.10  CONCLUSIONS 

After  reviewing  the  sampling  and  analysis  techniques 
in  relation  to  the  purpose  of  sampling  it  becomes  apparent 
that  river-beds  should  be  sampled  in  the  following  two  ways, 
unless  it  is  certain  that  the  bed  is  immobile: 

(1)  Scoop  or  core  sampling  at  some  depth 

below  the  surface,  and 

(2)  Grid  sample  obtained  by  a  taped  grid  or 

by  pacing. 

The  two  samples  so  obtained  should  be  adequate  for 
assessment  of  threshold  of  motion,  sediment  transport  and 
resistance  to  flow. 

The  scoop  sample,  which  is  volumetric,  should  be 
analyzed  by  weight  while  the  grid  sample,  which  is  areal, 
should  be  analyzed  by  number.  A  grid  sample  analyzed  by 
number  conveys  a  clear  picture  of  the  surface  material  on  a 
river  bed  as  the  probability  of  picking  up  a  certain  stone 
is  dependent  upon  its  surface  area. 


*• 
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After  reviewing  the  various  measures  of  represent¬ 
ative  size  it  appears  that  the  arithmetic  mean  is  the  most 
reasonable  one  but  is  time  consuming  to  use.  Equation  2.2 
has  been  used  by  the  author  and” found  to  be  convenient. 

The  standard  deviation,  equation  2.3,  should  be  used  as  a 
measure  of  dispersion.  Some  of  the  measures  noted  by 
sedimentary  petrographers  are  not  especially  useful  to  eng¬ 
ineers  . 

2.11  RECOMMENDATIONS 

Studies  should  be  conducted  to  arrive  at  the  most 
convenient  measure  for  representative  size  of  bed-material. 
The  relationship  between  the  median  size  by  weight  to  the 
median  size  by  number  should  be  investigated  further. 

The  variation  in  size  and  dispersion  of  bed-material 
with  the  type  of  morphological  feature  in  the  river  channel 
could  certainly  receive  more  attention.  The  various  types 
of  morphological  features  such  as  bars,  pools  and  riffjgs 
are  spatially  distinct  indicating  distinct  processes  which 
may  result  in  distinct  types  of  coarse  material  deposits. 


' 
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CHAPTER  III 

FIELD  OBSERVATIONS  ON  THE  THRESHOLD 
-  OF  MOTION  FOR  COARSE  BED-MATERIAL 

3.1  INTRODUCTION 

Recently  there  have  been  many  papers  dealing  with 
the  theoretical  and  experimental  aspects  of  the  threshold 
of  motion  of  particles  located  on  the  bottom  of  a  river  or 
a  channel  (ASCE  Task  Committee  on  Sedimentation  Manual  1966, 
Neill  1967,  Coleman  1967,  Egiazaroff  1967,  Neill  and  Yalin 
1969).  These  papers  have  presented  detailed  analysis  of 
the  phenomenon  as  investigated  in  laboratory  flumes.  An 
excellent  history  of  this  topic  has  been  presented  by 
Leliavsky  (1955)  . 

The  engineer  is  interested  in  the  threshold  of 
motion  of  coarse  bed-material  for  several  reasons  which 

can  be  grouped  into  two  categories: 

Category  A  -  Coarse  Mixture  on  Bed 

(1)  In  designing  stable  channels  the  channel  must  be 
designed  so  that  neither  its  bed  nor  its  banks  experience 
excessive  erosion. 

(2)  The  application  of  several  bed-load  transport 
formulas  in  rivers  requires  the  knowledge  of  a  critical 
shear  stress. 


.  ' 
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Category  B  -  Uniform  Bed-material 

(1)  In  designing  river  training  works  and  hydraulic 

structures  the  use  of  riprap  is  common.  The  riprap  must 
be  so  designed  that  the  surface  layer  of  stones  will  not 
move  under  design  conditions. 

The  design  of  riprap  protection  (Category  B)  has 
been  discussed  by  many  engineers  (Ramette  1963,  Blench  1969, 
Peterka  1963,  California  Dept,  of  Public  Works  1960,  Izbash 
and  Khaldre  1970)  and  will  not  be  dealt  with  here.  However, 
the  various  criteria  developed  for  the  threshold  of  motion 
for  coarse  mixtures  will  be  compared  with  criteria  that 
pertain  only  to  uniform  particles. 

In  this  chapter  an  attempt  will  be  made  to  arrive 
at  criteria  for  movement  of  coarse  mixtures  which  are  based 
on  flow  information  that  is  generally  available  to  river 
engineers  such  as  mean  flow  velocity  in  a  cross-section, 
surface  width,  mean-flow  depths,  bed  material  composition 
and  water  surface  slope.  Field  data  obtained  by  the 
writer  from  the  North  Saskatchewan  River  and  Wilson  Creek, 
as  well  as  other  investigators'  data  will  be  analyzed. 

3.2  THRESHOLD  OF  MOTION  FIELD  DATA 

During  the  summers  of  1965  and  1966  measurements 
of  the  largest  stones  moved  on  the  bed  of  the  North  Sask¬ 
atchewan  River  at  Drayton  Valley  were  undertaken  with  the 


. 
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hope  of  correlating  this  information  to  the  existing  flow 
parameters  (Galay  1967a) .  Subsequent  measurements  in  1968 
were  taken  on  Wilson  Creek  in  Manitoba  with  the  same  pur¬ 
pose  in  mind.  This  information  will  now  be  summarized 
along  with  similar  information  obtained  by  other  invest¬ 
igators  . 

3.2(a)  North  Saskatchewan  River  near  Drayton  Valley 

APPENDIX  2  presents  pertinent  geomorphic  and 
hydraulic  information  concerning  the  North  Saskatchewan 
River  near  Drayton  Valley.  Measurements  of  the  largest 
stones  moved  by  the  1965  peak-flow  were  carried  out  after 
the  flow  receded.  These  stones  were  located  on  open  bars 
adjacent  to  the  low  flow  channel  and  their  sizes  and 
locations  are  tabulated  in  TABLE  3.  All  the  s.tones  were 
considered  to  be  part  of  the  river-stone  population  and 
not  "erratics".  The  distribution  of  surface  material  is 
shown  in  FIGURES  4  and  7;  the  largest  stones  moved  cover 
approximately  0.01%  of  the  bar  surface.  It  would  appear 
that  this  1965  flood  was  capable  of  moving  virtually  all 
the  material  making  up  its  bed. 

Subsequent  to  the  stone  survey,  cobbles  having 
intermediate  diameters  of  2,  3,  4  and  6  inches  were  paint¬ 
ed  various  colours  and  placed  in  separate  "squares"  on  an 
open  bar,  (PHOTOGRAPH  6) .  This  procedure  was  carried  out 


. 
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TABLE  3 

SIZES  AND  WEIGHTS  OF  STONES 
MOVED  BY  1965  PEAK  FLOW 
NORTH  SASKATCHEWAN  RIVER  NEAR  DRAYTON  VALLEY 


SIZE 

WEIGHT 

Location 

Dominant 

(Intermediate  Axis) 

Approximate 

Lithology 

Average 

Maximum 

Average 

Maximum 

- 

(in.) 

(in.) 

(lb.) 

(lb.) 

U/S  of  #82 

Quartzite 

8.0 

10.0 

42 

75 

U/S  of  #86 

Quartzite 

9.8 

14.0 

60 

100 

U/S  of  #87 

Quartzite 

8.4 

11.0 

- 

- 

D/S  of  #88 

Quartzite 

7.5 

10.2 

35 

55 

U/S  of  #89 

Quartzite 

8.6 

11 .5 

43 

90 

#94 

Quartzite 

7.9 

11.0 

32 

50 

#96 

Quartzite 

11.0 

14.0 

75 

100 

#97 

Quartzite 

9.5 

10.0 

50 

55 

D/S  of  #101 

Sandstone 

10.5 

16.5 

80 

120 

' 
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before  the  1966  peak  flow,  with  the  hope  that  the  water 
would  rise  sufficiently  to  displace  some  of  the  cobbles 
from  the  bar  surface.  Velocity  measurements  directly  over 
the  squares  were  to  be  carried  out  with  the  hope  of  correlat¬ 
ing  the  maximum  measured  velocity  with  the  largest  cobble 
moved.  However,  the  peak  flow  in  1966  was  so  high  that 
more  than  just  the  top  layer  of  cobbles  were  moved;  in  fact, 
the  whole  bar  was  eroded  and  moved  downstream.  It  was 
apparent  that  stones  larger  than  six  inches  were  moved 
making  velocity  -  stone  size  correlations  impossible. 

More  data  on  the  threshold  of  motion  was  made  avail¬ 
able  when  a  sacked-concrete  revetment  at  x-sec  93  failed. 

The  river  flows  through  the  Pembina  Oil  Field  and  the  pur¬ 
pose  of  the  revetment  was  to  protect  a  water  injection  well. 
PHOTOGRAPH  7  a,  b,  and  c  shows  the  stages  of  failure  of  the 
revetment,  while  PHOTOGRAPH  8  shows  the  extent  of  the  channel 
shift.  After  the  high  flow  receded  many  of  the  concrete 
sacks  were  found  on  top  of  the  first  and  second  point  bars 
downstream  from  x-sec  93  as  shown  in  FIGURE  11  and  PHOTO¬ 
GRAPH  9.  They  were  moved  as  far  as  a  mile  and  were  located 
on  both  sides  of  the  main  flow  channel.  The  concrete  sacks 
were  measured  and  weighed  and  correlated  to  a  critical  mean 

velocity  as  obtained  by  Vm  =  Q/A  using  the  area  of  x-sec 

c 

94  as  a  reference  since  a  large  number  of  sacks  passed 
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through  this  section.  However,  the  concrete  sacks  would  be 
set  in  motion  whenever  they  were  eroded  from  the  revetment 
which  is  dissimilar  from  stones  on  a  river  bed.  Some  con¬ 
crete  sacks  were  not  moved  from  the  location  of  the  revetment; 
therefore  it  was  assumed  that  the  flow  was  just  sufficient 
to  move  the  sacks  (critical  state) .  The  data  are  tabulated 
in  TABLE  A-5,  APPENDIX  3. 

3.2(b)  Wilson  Creek,  Manitoba 

Wilson  Creek  is  an  experimental  catchment  on  the 
escarpment  of  the  Riding  Mountains  (MacKay  and  Stanton  1964) . 
In  the  summer  of  1968  a  beaver  dam  was  blown  up  in  the  upper 
reaches  resulting  in  a  flow  of  1,000  cfs  (100  yr.  frequency) 
passing  down  the  creek.  A  number  of  stations  were  establish¬ 
ed  to  obtain  the  following  data  (Newbury,  1968) ; 

(1)  The  peak  flow  level  and  the  local  water  surface 
slope  from  stage  rods  200  ft.  apart. 

(2)  The  local  surface  velocities  at  peak  flow  stage 
by  using  wooden  surface  floats. 

(3)  The  largest  stones  moved  by  the  peak  flow.  The 
largest  stones  were  located  accurately  in  the  200  ft.  reach 
and  painted.  After  the  flow  subsided  the  stones  were  check¬ 
ed  for  movement. 

(4)  Cross-sections  before  and  after  peak  flow  in  order 
to  check  for  aggradation  or  degradation  and  to  obtain  mean 
flow  depths. 


■ 
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TABLE  4  -  WILSON  CREEK 

AVERAGE  HYDRAULIC  AND  CROSS-SECTION 
PROPERTIES  AT  PEAK  FLOW 
STAGE  -  1968- 


STATION 

AVERAGE 
AREA  FLOW 

ft2(A) 

AVE. 

TOP  WIDTH 
ft2(bw) 

MEAN 

FLOW  DEPTH 
A/bw  (ft) 

W.S. 

SLOPE 

SURFACE  FLOAT 
VELOCITY 
(ft/s) 

T  2 
(lb/ f  t  ) 

No.  1 

89.5 

33.6 

2.7 

0.0350 

14.3 

5.92 

No.  4 

94.2 

30.8 

3.1 

0.0345 

11.6 

6.65 

No.  6 

62.8 

30.0 

2.1 

0.0055 

10.5 

0.72 
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TABLE  5  -  WILSON  CREEK 
DIMENSIONS  OF  MOVED  BOULDERS 
THAT  ARE  SMALLER  THAN  FLOW  DEPTH 


STATION 

MAJOR 

BOULDER  DIMENSIONS 

(ft.) 

INTERMEDIATE 

MINOR  AXIS 

No.  1 

2.3 

2.0 

1.6 

1.9 

1.8 

1.3 

2.2 

2.0 

2.0 

No.  4 

2.8 

2.5 

2.0 

2.0 

1.8 

1.5 

3.0 

3.0 

2.0 

3.7 

3.0 

2.6 

2.2 

1.9 

1.2 

2.5 

2.4 

2.4 

2.6 

2.4 

2.2 

3.2 

3.0 

3.0 

2.1 

2.0 

1.4 

3.1 

2.8 

2.2 

No.  6 

2 .2 

1.8 

•  1.8 

* 

2.4 

1.8 

1.6 

1.9 

1.8 

1.6 

1.9 

1.8 

1.6 

2.0 

1.8 

1.4 

1.6 

1.4 

1.4 

1.0 

0.9 

0.9 

1.7 

1.6 

1.4 

■  ' 
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TABLE  4  presents  the  average  hydraulic  and  cross- 
section  data  for  stations  1,  4  and  6,  as  they  were  the  only 
locations  having  stone  size  data  and  TABLE  5  gives  the 
sizes  of  moved  stones.  The  representative  boulder  size 
for  each  station  was  obtained  by  averaging  all  the  moved 
stones  that  were  smaller  than  the  depth  of  flow.  These 
were  then  ^tabulated  with  the  corresponding  mean  velocities 
in  TABLE  A- 5  (APPENDIX  3) .  The  mean  velocities  were 
obtained  by  multiplying  the  surface  float  velocities  by 
0.87  (after  Fahnestock  1963).  Several  stones  larger  than 
the  depth  of  flow  were  also  moved,  but  the  distance  moved 
was  very  short  suggesting  that  the  motion  was  by  "canting". 
The  bed  of  Wilson  Creek  contains  material  that  is  pre¬ 
dominantly  in  the  cobble  range  as  shown  on  the  bed- 
material  distribution  curves  (FIGURE  12) . 

3.2(c)  Middle  Fork  Eel  River  (Ritter  1967) 

The  Middle  Fork  Eel  River  flows  from  the  western 
slopes  of  the  Coastal  Ranges  in  Northwestern  California. 

A  study  was  undertaken  to  determine  the  maximum  size  of 
bed-material  that. will  be  moved  by  different  velocities 
of  flows.  At  intervals  of  10  -  20  feet  along  a  cross- 
section  at  a  gauging  station,  twelve  1-foot  squares  of 
bed-material  were  painted  and  located  with  reference  to 
the  right  bank.  A  flow  of  some  3,750  cfs  having  an 
average  velocity  of  6  ft.  per  second  removed  approximately 
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50%  of  the  stones  in  the  painted  squares.  The  stone  sizes 
and  the  corresponding  velocities  are  tabulated  in  TABLE 
A-5  (APPENDIX  3) .  The  bed-material  ranged  in  size  from 
sand  to  .boulders  and  averaged  1.0  inches. 

3.2(d)  Trucked  River/  Calif ornia-Nevada  (BirkeTand,  1968) 
Glacial  outwash  along  the  Truckee  River  was  found 
to  contain  large  boulders  and  gravel  bars  indicating  that 
exceptionally  large  flows  took  place  in  the  valley.  An 
attempt  to  estimate  the  mean  velocities  and  tractive 
forces  required  to  move  large  boulders  was  made  by  comput¬ 
ing  Manning's  roughness  coefficient  from  data  collected 
after  the  failure  of  the  St.  Francis  Dam  in  California. 
Manning's  n  was  computed  to  be  0.08  for  the  St.  Francis 
Dam  failure  and  0.05  to  0.03  for  present  day  flood  flows. 
Instead  of  using  n  =  0.08  a  roughness  value  of  0.06  was 
thought  to  be  reasonable  and  this  value  was  used  to  obtain 
velocities  as  tabulated  in  TABLE  A-5  (APPENDIX  3) .  The 
data  presented  here  may  not  be  as  reliable  as  some  of  the 
previously  mentioned  cases  as  no  direct  measurements  of 
actual  flows,  water  levels,  etc.,  were  possible. 

3.2(e)  White  River,  Washington  (Fahnestock,  1963) 

Data  on  the  velocities  required  to  transport 
coarse  materials  was  compiled  by  Fahnestock  in  a  study 
of  the  processes  of  valley  train  formation  by  a  proglacial 
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stream.  Velocities  were  measured  directly  by  a  current 
meter  or  with  floats  depending  on  the  depth  of  flow.  The 
corresponding  sizes  of  stone  in  motion  were  obtained  by 
catching,  the  stones  on  a  screen  or  by  sampling  on  nearby 
exposed  gravel  bars. 

3.2(f)  Rubicon  River/  California  (Scott  and  Gravlee,  1968) 

On  December  23,  1964  a  torrential  rainfall  on  the 
upper  part  of  the  Rubicon  River  caused  the  failure  of  the 
partly  completed  Hell  Hole  Dam.  The  resultant  surge  pro¬ 
duced  exceptionally  high  peak  discharges  which  eroded  the 
valley  walls  and  transported  a  vast  quantity  of  boulders. 
Water  level  and  depth  measurements  were  utilized  to 
compute  the  tractive  force,  and  the  mean  diameter  of  the 
10  largest  stones  at  the  corresponding  sites  were 
correlated  to  this  tractive  force.  The  data  are  tabulated 
in  TABLE  A- 5  (APPENDIX  3) . 

3.2(g)  Flooding  Rivers  in  Connecticut  (Wolman  and  Eiler  1957) 
A  flood  in  August  1955  inundated  a  large  portion 
of  the  valley  bottoms  in  Connecticut.  Erosion  of  channel 
and  valley  bottom  resulted  in  the  movement  of  cobbles  and 
boulders.  Reconnaissance  surveys  after  the  flood  re¬ 
sulted  in  a  compilation  of  the  sizes  of  boulders  moved 
in  channels  having  various  discharges,  slopes,  velocities, 
and  depths.  The  discharges  were  based  on  indirect  measure¬ 
ments  and  mean  velocities  were  obtained  by  dividing 
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discharge  by  cross-sectional  area.  The  slope  was  obtained 
from  topographical  maps.  There  were  some  twelve  differ¬ 
ent  sites  investigated  along  seven  different  rivers.  The 
data  are*  tabulated  in  TABLE  A- 5  (APPENDIX  3)  . 

3.2(h)  Coffee  Creek/  California  (Stewart  and  LaMarche  1967) 
On  December  22  and  23,  1964  a  flood  having  a 
magnitude  of  17,800  cfs,  unprecedented  in  the  110  year 
period  of  settlement,  occurred  on  Coffee  Creek  destroying 
large  areas  of  meadowland.  Cross-sections  surveys  were 
conducted  to  assess  the  erosion  and  deposition  features  of 
the  flood.  From  these  cross-sections  and  a  knowledge  of 
the  peak  discharge  a  mean  velocity  was  computed  and 
correlated  to  the  maximum  size  stones  moved  in  the  vicinity 
of  the  cross-section. 

3.2(i)  Neri  River,  Japan  (Oishi  1956) 

In  order  to  assess  the  transport  of  gravel  and 
boulders  in  mountainous  streams  a  number  of  stones  were 
charged  with  cobalt-isotope  with  the  intention  of  finding 
these  stones  after  recession  of  the  high  flows.  In 
1953,  a  typhoon  struck  the  area  causing  a  small  flood  to 
take  place.  The  movement  of  the  stones  was  somewhat 
irregular  during-  the  flood.  Measurements  of  the  mean 
velocity,  depth  and  slope  were  obtained. 
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3.2  ( j )  Gravel  Rivers  in  Western  Alberta  (Van  Per  Giessen  1966) 

Four  coarse-bed  rivers,  namely  the  Castle,  Sheep, 

Elbow  and  Drywood  Creek,  were  investigated  in  order  to 

assess  their  character.  Colored  river-bed  stones  were  used 

to  study  the  threshold  of  motion  condition  and  the  distance 

of  transport  of  bed-material  in  the  Elbow  and  Drywood  Creek. 

It  was,  however,  difficult  to  sort  out  the  data  presented  in 

this  thesis;  plots  of  particle  size  DQA  versus  tractive 

y  u 

force  are  shown,  but  the  computed  tractive  force  was  for  a 
representative  flow  of  a  smaller  magnitude  than  the  maximum 
flow  that  took  place  in  the  various  rivers.  The  data  were 
adjusted  and  are  tabulated  in  TABLE  A-5  (APPENDIX  3). 

3.2(k)  Elbow  River,  Alberta  (Hollingshead  1968a) 

Field  investigations  were  carried  out  on  the 
Elbow  River  near  Bragg  Creek  in  order  to  learn  more  about 
bed-load  transport  in  coarse-bed  rivers  and  to  assess  the 
behaviour  of  this  river.  Hollingshead  (1968a)  assessed  the 
threshold  of  motion  state  using  values  of  bed  shear  stress 
toc  as  determined  from  velocity  profiles.  The  values 
used  in  this  study  are,  however,  obtained  by  using  to  =  yd*S 
in  order  to  be  comparable  to  the  previously  mentioned  invest¬ 
igations  . 

3.2(1)  San  Luis  Valley  Canals  (Lane  and  Carlson  1953) 

A  number  of  stable  coarse-bed  canals  were  in¬ 
vestigated  in  order  to  establish  criteria  for  design  of 
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stable  canals.  In  this  case  the  bed-material  was  assumed  to 
be  at  the  threshold  of  motion  state  for  the  maximum  sustained 
flows . 

3.2 (m)  Coarse-bed  Rivers  in  British  Columbia  (Kellerhals  1963) 
Seven  stable  reaches  were  investigated  in  order  to 
obtain  regime  type  formulas.  Channel  dimensions,  slopes  and 
bed-material  data  were  collected  with  formative  or  bank-full 
discharges  being  estimated;  the  bed-material  was  assumed  to 
be  at  the  threshold  of  motion.  The  data  are  tabulated  in 
TABLE  A- 5  (APPENDIX  3) . 

3.3  ANALYSIS  OF  FIELD  DATA 

The  field  data  described  in  the  previous  section  in¬ 
volving  bed-material  mixtures  will  be  analyzed  in  the 
following  ways : 

(a)  Mean  critical  velocity  related  to 
representative  stone  size. 

(b)  Critical  shear  stress  related  to 
representative  stone  size. 

(c)  Shields  criteria  relating  a  tractive 
force  coefficient  to  the  particle  Reynolds 
number . 

(d)  Shields  tractive  force  coefficient  re¬ 
lated  to  the  relative  depth  of  the  channel. 

(e)  Froude  number  related  to  the  relative 
depth  of  the  channel. 
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The  relationships  obtained  from  the  various  treat¬ 
ments  will  then  be  compared  with  corresponding  relationships 
for  uniform  materials. 

3*3(a)  Mean  Critical  Velocity  Related  to  Stone  Size 

One  of  the  earliest  attempts  to  assess  the  threshold 
of  motion  for  particles  on  the  bed  of  a  channel  was  by  Brahms 
(1753)  who  presented  the  equation: 

Vmc  =  a  Wp1/6  (3.3.) 

where : 

Vmc  =  critical  velocity 
a  =  empirical  constant 
Wp  =  weight  of  the  particle 
This  is  the  simplest  approach  to  the  problem  but  it  dis¬ 
regards  several  significant  factors  such  as  the  depth  of 
flow.  The  equation  can  be  modified  by  assuming  a  spherical 
particle  which  leads  to  the  weight  being  proportional  to  the 
cube  of  the  diameter  and: 

Vmc  =  a  D1//2 

where : 

D  =.  particle  diameter 

A  plot  of  the  mean  critical  velocity  versus  the 
corresponding  stone  size  is  shown  in  FIGURE  13.  It  should 
be  noted  that  the  velocity  is  the  mean  velocity  for  a  cross- 
section,  not  the  mean  velocity  in  a  vertical  at  some  designat¬ 
ed  point  in  a  channel.  Also,  the  stone  size  is  generally  the 
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mean  of  the  largest  particles  that  have  been  moved;  in  some 
cases  the  stone  size  has  been  taken  as  the  D^q  size,  where 
90%  of  the  particles  are  finer  than  the  specified  size.  The 

t 

D90  particle  size  was  chosen  for  rivers  or  canals  that  had 
developed  a  stable  non-moving  bed. 

The  numbers  adjacent  to  the  plotted  points  are  mean 
depths  of  flow. 

A  tentative  equation  representing  the  field  data 
for  bed-material  mixtures  would  be: 

Vmc  =  8.0  D1/3  (3.2) 

The  scatter  shown  on ' the  plot  could  be  attributed  to  a  number 
of  different  factors: 

(1)  Difficulty  in  obtaining  accurate  flow 
data  during  high  flow  conditions. 

(2)  With  low  values  of  relative  depth  the 
flow  would  be  highly  turbulent  with  a 
possibility  of  high  concentrations  of  fine 
particles.  This  condition,  which  may  be 
termed  a  mudflow,  was  aptly  described  by 
Richardson  (1968)  : 

they  saw  the  moving  mass  of  water  and  rock 
coming  over  the  end  of  the  glacier  . . . 
described  the  material  as  similar  to  a 
huge  mixture  of  concrete  except  that  it 
was  darker  in  color.  They  stated  that  the 
force  was  so  great  that  immense  boulders 
were  thrown  from  ten  to  thirty  feet  into 
the  air  as  the  mass  moved  forward. 
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Under  these  circumstances  the  values  of 
flow  density  p,  specific  weight  y,  and 
dynamic  viscosity  y  would  be  different 
from  pure  water. 

(3)  The  representative  particle  size  was 

* 

not  consistent  in  all  cases  as  some  in¬ 
vestigators  listed  only  the  one  largest 
stone  moved  while  others  took  an  average 
of  a  number  of  large  stones. 

It  would  now  be  in  order  to  compare  equation  3 . 2 

with  similar  equations  for  uniform  size  particles.  Bhowmik 

and  Simons  (1970)  compared  the  equations  obtained  by  Smith 

and  Hallmark  (1965) ,  the  USBR  (Peterka  1963) ,  and  Izbash  (1936) 

as  shown  in  FIGURE  14.  All  the  equations  are  of  the  form 
1/2 

V  «  D  and  the  following  design  equation  was  proposed  by 

Bhowmik  and  Simons  (1970) : 

V  =  10.4  D1/2  (3.3) 

me 

The  equations  shown  in  FIGURE  14  are  based  on  uniform  part¬ 
icles,  however,  some  of  the  plotted  points  pertain  to  bed- 
material  mixtures. 

Equation  3 . 2  for  bed-material  mixtures  is  compared 
with  the  various  equations  for  uniform  materials  in  FIGURE  15. 
The  curve  for  the  mixtures  is  somewhat  flatter  and  intersects 
the  curves  for  uniform  material  around  the  one-foot  stone 
size.  It  would  appear  that  higher  velocities  are  required 
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to  initiate  motion  of  mixtures  in  the  gravel  size  range. 

The  converse  applies  in  the  boulder  range;  the  D50  size  will 
be  in  motion  at  lower  velocities  for  mixtures  as  compared  to 
uniform  particles . 

Neill  (1967)  also  compared  various  other  equations 
which  could  be  reduced  to  a  relationship  between  velocity 
and  stone  size  if  a  flow  depth  is  known.  The  equations 
compared  here  are  those  of  Mavis  and  Laushey  (1948)  ,  USSR/ 
Lane  (1955)  ,  Sundborg  (1956)  ,  California  Highways  (1960)  and 
Straub  (1953).  The  velocity  used  in  these  equations  is, 
however,  not  the  mean  velocity  in  a  cross-section  but  the 
mean  in  a  vertical  which  is  assumed  to  be  adjacent  to  the 
moving  stones.  Equation  3.2  for  bed-material  mixtures  is 
compared  to  the  various  equations  in  FIGURE  16  for  material 
having  a  specific  gravity  of  2.65  and  a  depth  of  one  metre 
(3.28  feet).  An  equation  by  Neill  (1968b)  is  also  presented 
in  this  plot. 

The  curve  for  bed-material  mixtures  does  not  agree 
closely  with  the  curves  from  other  investigators.  However, 
this  is  not  surprising  since  a  mean  cross-section  velocity 
was  used  in  arriving  at  equation  3.2  while  a  mean  velocity 
in  a  vertical  was  used  by  the  investigators  noted  on 
FIGURE  16.  The  slope  of  equation  3.2  does  agree  quite  well 
with  the  slopes  of  the  other  curves  (except  for  the 
California  Highways  curve) . 
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3.3(b)  Critical  Shear  Stress  Related  to  Stone  Size 

The  critical  shear  stress  has  been  used  as  the 
basis  for  stable  canal  design  in  coarse  materials  (Lane, 
1955).*  A  plot  of  the  critical  shear  stress  versus  stone  size, 
shown  in  FIGURE  17  reveals  considerable  scatter;  no  attempt 
was  made  to  obtain  an  equation  from  this  plot. 

Several  relations  obtained  by  other  investigators, 
namely  Lane  (1955)  ,  Kellerhals  (1963)  ,  Egiazaroff  (1965)  and 
Komura  (1967)  are  added  to  FIGURE  17.  It  is  somewhat 
difficult  to  make  meaningful  comparisons,  however,  as  most 
of  the  curves  are  related  to  a  size  while  the  plotted 

data  do  not  correspond  to  any  one  frequency  value.  It  “ 

would  appear  that  the  relation  established  by  Lane  (1955)  is 

the  most  reasonable.  The  data  indicates  that  the  1  on  1 

slope  of  the  various  curves  may  be  too  flat;  more  reliable 

data  would  be  required  to  firmly  establish  the  slope. 

3.3(c)  Shields  Tractive  Force  Coefficient  Related 
to  Particle  Reynolds  Number 

The  most  common  criteria  for  the  threshold  of  motion 

is  based  upon  the  Shields  diagram  (Shields  1936) .  Shields 

plotted  data,  from  flume  experiments,  on  a  dimensionless 

plot  having  a  "tractive-force  coefficient"  XC//y^D  versus  the 

particle  Reynolds  number  v*D/^  where  xc  =  the  critical  shear 

stress  at  the  bed,  y'  =  the  bouyant  weight  of  the  bed- 

s 

material,  D  =  particle  diameter,  v*  =  shear  velocity  and 
v  =  kinematic  viscosity.  The  "tractive-force 


. 


*■ 


V 


' 


I 


43 


coefficient"  has  also  been  termed  an  "entrainment  function" 
(Randkivi  1967)  ,  a  "mobility  number"  (Yalin  1965)  ,  and  a 
"dimensionless  critical  shear  stress"  (ASCE  Task  Committee 
on  Sedimention  Manual  1966).  Since  the  term  "tractive-force 
coefficient"  has  not  been  used  extensively  by  research  work¬ 
ers  and  the  mobility  of  a  particle  is  related  to  the  dimen¬ 
sionless  parameters,  the  term  "mobility  number"  will  be  used 
in  the  remaining  discussion. 

Shields  used  physical  reasoning  to  arrive  at  his 
dimensionless  parameters.  However,  it  is  possible  to  obtain 
these  same  parameters  by  using  dimensional  analysis  and 
assuming  that  the  threshold  of  motion  is  determined  by 
tc,  Y's ,  D,  p  and  y  where  p  is  the  fluid  density  and  y  is 
the  dynamic  viscosity  (ASCE  Task  Committee  on  Sedimentation 
Manual  1966) .  The  choice  of  variables  by  the  Task  Committee 
is  somewhat  incomplete;  a  more  complete  presentation  in- 
eluding  all  the  pertinent  variables  is  presented  in  section 
3.3(e).  Shields  original  diagram  is  shown  in  FIGURE  18  with 
a  band  through  the  plotted  points.  The  influence  of  the 
particle  Reynolds  number  becomes  negligible  for  v*D/^>  70 
(Yalin  1966)  and  it  would  appear  that  the  mobility  number 
would  have  a  value  ranging  from  0.04  to  0.06. 

The  field  data  is  plotted  on  FIGURE  19  and  shows 

2 
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a  wide  range  of  mobility  numbers,  — — *—  .  The  values  of  the 
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mobility  number  are  plotted  on  a  histogram,  FIGURE  20;  the 
most  frequently  occurring  value  (mode)  is  0.015  while  the 
mean  of  all  the  obtained  mobility  numbers  is  0.042.  The  use 
of  a  mobility  number  equal  to  o'.  03,  also  recommended  by 
Neill  (1968b)  for  uniform  particles,  would  appear  to  be 
reasonable . 

3.3(d)  Shields  Tractive  Force  Coefficient  Related  to 
the  Relative  Depth  of  the  Channel 

In  using  dimensional  analysis  to  formulate  equat¬ 
ions  for  the  threshold  of  motion  it  has  been  argued  that  the 
mean  depth  of  flow  d*  as  well  as  particle  density  ps  should 
be  included  as  important  variables  (Neill  and  Van  Der  Giessen 
1966,  Yalin  1965).  This  would  result  in  the  following: 

tc  =  f  i  (y's,  D,  p,  y,  ps,  d*)  (3.4) 

The  critical  shear  stress  tc  can  be  replaced  by  v*  from 
v*  =  /xc7p ,  which  leads  to: 


pvV_=  f  ,  PV*D  f  p^ 

Y'sD  D*”  p 


(3.5) 


Within  the  range  of  existing  experimental  data  the  density 
ratio  ^s/p  appears  to  be  insignificant  and  for  coarse  part¬ 
icles  the  influence  of  particle  Reynolds  number  may  not  be 
influential  ( pv*D-  >  70)  .  This  results  in: 

y 


px*L  =  f2  (^t) 

YS'D  ' 


(3.6) 


The  form  of  the  mobility  number  was  changed  slightly  by 
Neill  and  Van  Der  Giessen  (1966)  by  replacing  v*  with  the 
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mean  critical  flow  velocity  V  : 

J  me 


PV, 


v'  D 
'  s 


me  _  f 

3  D  ; 


(3.7) 


The  collected  field  data  were  plotted  according  to 
the  above  relationship  (FIGURE  21)  and  shows  considerable 
scatter.  Neill's  (1968)  equation: 


pVmc  _  9  n  ,d^l/3 

y'  D  Z,U  } 


(3.8) 


is  also  shown  on  the  plot  but  does  not  pass  through  the 

plotted  points.  This  is  to  be  expected,  however,  as  the 

definition  of  mean  velocity  is  not  identical  -  the  field  data 

is  based  on  mean  cross-section  velocities  while  Neill's 

(1968)  equation  is  based  on  mean  velocities  in  the  vertical. 

3.3  ( e )  Froude  Number  Related  to  the  Relative 
Depth  of  Flow 

In  using  dimensional  analysis  there  are  many 
possible  combinations  of  the  variables  involved.  The  follow¬ 
ing  variables  are  involved  in  the  threshold  of  motion: 

T  ,  v*  or  V  =  f  (p,y,D,p  ,cr  ,fc,fg,b,d*,Cw,g) 

(3.9) 


where  the  new  variables  are: 

-  shape  of  bed-material  distribution  curve 
fc  =  factor  defining  cross-sectional  shape  of  channel 

fg  -  factor  defining  plan  geometry  of  channel 
b  =  width  of  channel 

Cw  =  concentration  of  suspended  particles 


"r 


Using  ^mc  as  the  dependent  variable  conventional 
dimensional  analysis  leads  to: 

Vmc  =  fi  ( pvmrD /  b  ,  d*,  q  ,  fc,  fg,  Cw,  ps  )  .  . 

✓P7  (  y  d*  D  ^  .  “P  )  (3'10) 

All  the  flume  and  field  data  have  been  obtained  from 
channels  that  were  relatively  straight,  wide  (b>5d*)  and  of  a 
similar  cross-section  which  would  remove  the  b/d* ,  fc  and 
fg  parameters  from  equation  3.10. 


The  effects  of  the  following  parameters  a,  ,  and 

•U  P 


Cw  are  assumed  to  be  insignificant,  primarily  because  they 
have  not  been  adequately  studied  to  date  and  equation  3.10 
becomes : 


V 


me  =  f 2  (pVmcD,  d*) 
/gd*  ( 


y 


d  ) 


(3.11) 


This  arrangement  is  very  similar  to  an  arrangement 
proposed  by  Blench  (19  69)  : 


V 

me 

7gd„ 


=  f  ( 3 /vgD ,  d*) 
(  v  D  ) 


(3.12) 


for  conditions  of  bed-load  charge  approaching  zero.  The 

term  ?vgD  is  similar  to  the  particle  Reynolds  number  and  is 
v 

called  the  particle  Vig  Number.  It  may  be  argued  that  the 


Reynolds  number  is  insignificant  for  fully  rough  turbulent 
flows  which  results  in  equation  3.11  becoming: 


. 

'■  . 
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Blench  (1970)  deletes  the  g  term  from  the  Froude  number  and 
arrives  at  a  zero  bed  factor  F, 


bo 


2 

'mc/a 


=  V  /  ;  equation  (3.13) 


can,  therefore,  be  modified  to: 


V 


IT  =  Fbo  =  fl* 

*  D°  (D  ) 


0.14) 


The  field  data  for  threshold  of  motion  is  plotted 
with  F^  versus  d*/^  in  FIGURE  22  and  exhibits  the  same  degree 
of  scatter  as  the  plot  with  the  modified  mobility  number 
(FIGURE  21).  The  equation: 

Fbo  =  29  <D/d*)V2  (3.15) 

presented  by  Blench  (1967)  is  plotted  on  the  figure  and 

appears  to  represent  the  plotted  points  reasonably  well. 


3.4  CONCLUSIONS 

After  examining  the  plots  of  the  river  data  it 
would  appear  that  the  relationship  between  mean  velocity  and 
stone  size: 

V  =  8.0  D1/3  (3.2) 

me 

can  be  used  to  assess  when  large  stones  in  a  mixture  will 
begin  to  move.  The  equation  is  inadequate  academically  but 
represents  the  data  as  well  as  the  more  refined  plots  which 
include  more  variables.  A  stable  natural  coarse-bed  channel 
could  therefore  be  designed  so  that  the  resultant  mean  velo¬ 
city  will  be  less  than  that  indicated  by  equation  3.2,  using 
D50  as  the  representative  diameter. 

'  Equation  3.2,  for  bed-material  mixtures  does  not 
agree  closely  to  equations  developed  for  uniform  materials 


' 

»■ 


48 


which  are  of  the  form  Vmc  «  D1/2.  It  would  appear  that 
higher  velocities  are  required  to  move  a  Dr^  size  for  mix- 
tures  as  compared  to  uniform  materials  in  the  gravel  size 
range.  The  converse  is  the  case  in  the  boulder  range. 

The  critical  tractive  force  criterion  as  developed 
by  Lane  (1955)  also  appears  to  correspond  reasonably  well  to 
the  field  data  plotted  in  this  study. 

The  various  threshold  criteria  that  can  be  developed 
through  dimensional  analysis ,  such  as  Shields  entrainment 
function  or  mobility  number  and  the  regime  theory  zero  bed 
factor,  were  generally  reduced  to  a  relationship  between  a 
flow  intensity  parameter  and  relative  depth.  The  various 
plots  to  test  these  relationships  indicated  considerable 
scatter.  However,  an  equation  presented  by  Blench  (1967) 
appears  to  fit  the  data  reasonably  well: 

Fbo  =  29  (D/d*) 1/2  (3. 15) 

3.5  RECOMMENDATIONS 

The  procurement  of  more  field  data  on  the  threshold 
of  motion  under  well  controlled  conditions  is  desirable. 

The  radioactive  labelling  of  various  sizes  of  coarse  material 
in  straight  canals  such  as  the  San  Luis  Valley  canals  and  the 
observation  of  the  corresponding  state  should  yield  valuable 
data.  Hydraulic  parameters  would  be  relatively  easy  to 
obtain  and  the  movement  of  labelled  stones  could  be  followed 
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closely  by  appropriate  detectors. 

Laboratory  studies  to  assess  the  threshold  state 
for  bed-material  mixtures  would  also  be  of  value.  Almost 
all  laboratory  studies  conducted  to  date  have  dealt  with 
uniform  materials.  Threshold  studies  should  be  conducted 
with  the  bed-material  naturally  sorted  by  flowing  water  as 
is  the  case  in  nature.  The  bed  should  not  be  placed  and 
molded  by  templates  or  screeds.  However,  the  assessment  of 
the  threshold  state  would  probably  be  very  difficult  with 
mixtures  as  fine  material  in  motion  may  obscure  the  observat¬ 
ion  of  the  bed-material.  It  would  then  be  necessary  to 
label  certain  particles  by  paint  or  radioactive  tracers  and 
inspect  the  bed  after  each  test  run.  These  laboratory  tests 
should  be  carried  out  for  several  distinctly  different  bed 
mixtures  using  various  density  materials.  These  experiments 
would,  however,  be  very  complex  and  time  consuming. 


- 
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CHAPTER  IV 


THE  TRANSPORT  OF  COARSE 
BED-MATERIAL  IN  FLUMES  AND  RIVERS 

4.1  INTRODUCTION 

To  date,  most  of  the  laboratory  research  and  field 
investigations  concerned  with  sediment  transport  have  been 
confined  to  relatively  fine  materials,  namely  sand  (between 
0.062  and  2.0  mm  in  size).  Recently,  however,  the  construc¬ 
tion  of  engineering  projects  adjacent  to  coarse-bed  mountain 
rivers  has  resulted  in  a  number  of  questions  that  require 
immediate  answers .  Some  of  the  more  important  questions 
would  be : 

(1)  How  quickly  will  reservoirs  fill  up  if 
most  of  the  material  in  transport  is 
coarse? 

(2)  How  much  degradation  will  take  place  in 
the  river  channel  below  dams? 

(3)  How  much  scour  will  take  place  near 
obstacles  such  as  bridge  piers, 
abutments,  and  river  bends? 

(4)  What  is  the  best  design  approach  for  a 
stable  channel  that  is  required  to  pass 
large  quantities  of  coarse  material? 

The  answers  to  some  of  these  questions  may  exist 
for  sand-bed  rivers  in  the  form  of  empirical  equations,  but 
the  extrapolation  of  these  equations  to  rivers  having  bed- 
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material  as  large  as  four  to  six  inches  may  be  unwise.  The 
answers  depend  upon  an  intimate  knowledge  of  the  manner  and 
the  rate  of  bed-material  transport. 

The  manner  of  transport  of  coarse  bed-material  may 
differ  somewhat  from  sand  as  the  fall  velocity  law  is 
different  for  the  two  materials.  Another  important  factor 
may  be  the  effect  of  Reynold's  number  (temperature)  which 
for  sand  channels  greatly  influences  the  amount  of  material 
in  transport.  This  may  not  be  the  case  for  coarse-bed 
channels . 

Laboratory  experiments  dealing  with  coarse  bed- 
material  transport  have  been  carried  out  by  Gilbert  (1914)  , 
Meyer-Peter  and  Muller  (1948),  U.S.  Waterways  Experiment 
Station  (1935)  ,  Liu  and  Carter  (1935)  and  Bogardi  and  Yen 
(1938).  These  data  will  be  analyzed  after  examining  the 
transport  phenomena  from  the  dimensional  analysis  approach. 

A  summary  of  the  sediment  properties  is  shown  in  TABLE  6. 

The  application  of  the  presently  available  data  to 
the  transport  of  coarse  material  in  rivers  will  also  be 
discussed. 

4.2  DIMENSIONAL  ANALYSIS 

The  movement  of  coarse  bed-material  along  the  bed 
of  a  river  or  channel  is  a  two-phase  phenomenon  with  the  bed¬ 
load  transport  dependent  upon  the  flow  conditions,  the 
channel  geometry  and  the  properties  of  the  fluid  and  bed- 
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TABLE  6 

SUMMARY  OF  SEDIMENT  PROPERTIES 
COARSE  BED-LOAD  EXPERIMENTS 
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material  (see  FIGURE  23) . 

The  variables  can  be  assembled  in  the  following 

manner : 

Flow  conditions: 

qs  =  weight  of  bed-load  transported  per  unit  time  per 
unit  width 

d*  =  mean  flow  depth  equal  to  flow  area  divided  by 
water  surface  width 

S  =  slope  of  energy  gradient 

Cw  =  concentration  of  suspended  particles 

Channel  geometry: 

b  =  width  of  channel 

fc  =  factor  defining  cross-sectional  shape  of  channel 
fg  =  factor  defining  plan  geometry  of  channel 
Fluid  properties : 

p  =  density  of  fluid 
]i  =  dynamic  viscosity  of  fluid 
Bed-material  properties : 

D  =  represnetative  size  of  bed-material 
p  =  density  of  bed-material 

b 

0^  =  shape  of  sediment  distribution  curve  or  gradation 
of  bed-material 

=  shape  factor  of  bed-material 
and  g  =  acceleration  of  gravity. 
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The  variables  can  be  arranged: 

qs  =  f1(d*/s/Cw/b,fc/fg/p,y,D,ps/ab,a:b,g)  (4.1) 

It  is  possible  to  express  certain  variables  in  terms  of 
others;  g  and  S  can  be  replaced  by  y1 ,  the  submerged  specific 
weight  of  the  bed-material,  and  v* ,  the  shear  velocity,  using 
the  relationships: 

y's  =  g  (ps  '  p) 

'  and 

v*  = 

The  bed-load  transport  is  therefore  a  function  of: 
qs  =  f-L  (d*,v^,Cw,b,fc,fg,p,y,D,ps,ab,ab,Y^)  (4.2) 

The  bed-load  transport  is  dependent  upon  thirteen 
variables,  some  of  which  will  have  to  be  discarded  if  this 

complex  phenomenon  is  to  be  dealt  with.  The  concentration  of 
fine  suspended  material  has  a  significant  influence  on 
the  transport  of  sand  (Simons,  Richardson  and  Hauschild  1963) 
but  little  information  exists  on  its  influence  on  the  trans¬ 
port  of  coarse  bed-material.  For  this  reason  this  variable 
will  be  neglected. 

The  channel  width  b  has  minimal  influence  on  flow 
conditions  if  the  width  is  more  than  five  times  the  flow 
depth.  If  the  width  is  narrower  a  sidewall  correction  can 
be  applied;  this  variable  can  therefore  be  deleted. 

The  existing  data  for  coarse  material  transport  is 
confined  to  flumes  or  rivers  having  straight  alignments  with 
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reasonably  regular  cross-section  shapes.  The  variables  fc 
and  fg  are  therefore  removed. 

The  influence  of  the  gradation  and  shape  of  the  bed- 
material  have  not  been  studied  in  detail.  It  is  necessary  to 
assume  that  the  various  materials  used  in  the  flume  tests 
were  similar  in  their  gradation  and  shape. 

The  removal  of  these  variables  results  in: 
q  =  f  (d* ,v* , p ,y  ,D  ,  p  ,y')  (4.3) 

Using  the  Buckingham  ir-theorem  the  following  dimen¬ 
sionless  parameters  can  be  obtained  (choosing  p,  D,  and  v* 


as  repeating  variables) : 


Q.q  =  f  ( Pv* _ 

pv*3  X3  '  Yg  p' 


pDv 


-) 


(4.4) 


y  p 

The  relative  density  parameter  Ps/p  has  been  shown 
to  be  insignificant  for  threshold  of  motion  conditions 
(Gessler  1965)  and  the  same  conclusion  is  assumed  to  hold 
for  the  transport  of  material.  Equation  4.4  then  reduces  to 


qg  _  f  /PDv*.  p.v-iJ:  sL 

pv.  3  ”  r4  1  y  '  Y:  D'  D 


•) 


*  “  -  r-  .  g 

For  convenience  the  above  dimensionless 
will  be  defined  as  follows: 


P  = 

X  = 

Y  = 


4s/pv* 

Ppv*_ 

y 

Pv*2- 


y'  d 

'  s 


Z  =  d. 


/D 


(4.5) 

variables 
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The  most  compact,  although  not  the  most  convenient, 

form  of  presentation  would  be  a  three-dimensional  plot  as 

d . 

shown  in  FIGURE  24.  If  the  dimensionless  variable  ^  has 

a  significant  influence  on  the  transport  rate  several 

"surfaces"  will  be  evident.  It  would  also  be  more  convenient 

2 

to  change  the  variable  X  to  X  /Y  which  yields: 

X2  _  ys  D3 
Y  pv2 

In  flume  experiments,  with  water  temperatures  kept  constant, 

2 

the  above  X  /Y  variable  will  depend  only  on  the  type  and 

size  of  sediment.  The  sediment  fed  into  the  flume  can  be 

2 

controlled  at  will,  resulting  in  a  constant  value  of  X  /Y . 


.  The  final  relationship  is  therefore: 


=  fc  ( 


ysD 


Pv, 


7  Y.I  D  '  D 


) 


(4.6) 


5  pv2 

This  relationship  is  identical  to  that  of  Yalin 
(1965)  and  similar  to  that  of  Cooper  and  Peterson  (1968)  who 


arrived  at: 

=  F  ,VvgD  ,  C,  (4.7) 

gd*  1  V  D 

where:  C  =  bed- load  charge  in  parts  per  hundred 
thousand  by  weight.  Cooper  and  Peterson  (1968)  ,  however, 
considered  C  as  independent  and  imposed  on  the  system  and 


obtained  a  Froude  number  to  relate  the  intensity  of  flow 
instead  of  a  mobility  number  Y. 
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4.3  ANALYSIS  OF  FLUME  DATA 


The  data  from  the  flume  experiments  for  particles 
over  2  mm  in  size  are  tabulated  in  APPENDIX  4  and  will  serve 
as  a  basis  for  checking  the  dimensional  analysis. 

However,  slight  changes  in  the  variables  were 
necessary  for  several  reasons.  Side  wall  corrections  were 
necessary  due  to  the  difference  in  surface  roughness  between 
glass  or  steel  flume  walls  and  a  gravel  bed.  The  correction 
techniques  initially  developed  by  Einstein  (1942)  and  Johnson 
(1942)  and  modified  by  Vanoni  and  Brooks  (1957)  were  used  in 
analyzing  the  available  data.  This  sidewall  correction  re¬ 
sulted  in  a  change  in  the  hydraulic  radius  R  which  is  now 
designated  as  the  hydraulic  radius  of  the  bed  R^;  the 

friction  factor  f  becomes  the  friction  factor  of  the  bed  f, 

b 

and  the  shear  velocity  v*  now  pertains  to  the  bed  and  is 

noted  as  v*  .  The  various  dimensionless  variables  now  be- 
b 

come : 


Y 


pv*b2 

y'D 

'  S 


and 


Z  =  R^/D  (Rb  replaced  d*) 

X!  _  YS  D3 
Y  pv2 


The  first  group  of  plots  of  P  versus  Y  for  various 
values  of  X  /Y,  or  in  effect  for  different  gravel  sizes,  are 
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shown  in  FIGURES  25  to  34  and  represent  flume  tests  with  bed- 
material  of  uniform  size.  These  plots  indicate  that  the 
relative  depth  Rb/D  has  a  distinct  influence  on  the  transport 
of  the  coarse-bed  material.  Lines  having  equal  values  of  Rb/D 
ranging  from  10  to  30  were  established  on  these  figures. 

The  effect  of  increasing  the  depth  and  consequently  the 
value  of  Rb  is  to  decrease  the  rate  of  material  in  transport. 
It  is  apparent  that  doubling  the  flow  depth,  for  one  sediment 
size  and  the  same  shear  velocity,  can  result  in  a  reduction 
of  the  bed-load  transport  by  a  factor  of  four  or  five. 

These  ten  plots,  from  FIGURES  25  to  34,  all  per¬ 
tain  to  uniform  material  and  a  subsequent  plot  of  P  versus 

Y  for  a  constant  value  of  Rb/D  =  15  (FIGURE  35)  indicates 

2 

that  the  parameter  X  /Y  or  the  particle  size  is  also 
significant.  This  plot  shows  that  the  relation  between  the 
mobility  number  Y  (or  the  flow  intensity)  and  the  sediment 
transport  parameter  P  is  rather  complex;  it  would  be  diffi¬ 
cult  to  develop  equations  relating  these  parameters  since 
the  resultant  plot  shows  a  three  dimensional  curved  surface. 

A  similar  plot  (FIGURE  36)  for  a  different  relative  depth 
value  will  yield  a  surface  in  a  slightly  different  position 
from  the  first  plot  and  if  more  data  existed  a  family  of 
surfaces  could  be  obtained.  This  family  of  surfaces  should 
describe  the  phenomena  completely  for  the  sizes  of  material 
used  in  the  flume  tests,  however,  these  surfaces  may  not  be 
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applicable  to  coarse-bed  river'  since  the  range  of  particle 
sizes  found  on  the  beds  of  these  rivers  is  very  wide.  Sur¬ 
faces  developed  from  flume  tests  using  bed— material  mixtures 
would  be  required  in  this  case. 

In  this  regard ,  data  from  flume  tests  by  Meyer-Peter 
and  Muller  and  other  investigators  are  plotted  in  FIGURES  37 
to  43.  The  same  type  of  plots  with  lines  of  constant  re¬ 
lative  depth  are  obtained,  only  in  this  case  the  surfaces 
are  higher  than  in  the  previous  plots  for  uniform  materials. 
The  data  by  Liu  and  Carter,  however,  plot  rather  inconsistent¬ 
ly  -  this  is  not  surprising  as  some  of  the  tests  from  this 
series  show  the  sidewalls  of  the  flume  to  be  rougher  than 
the  bed.  These  tests  exhibit  low  shear  velocity  values 
which  would  be  the  result  of  errors  in  either  slope  or  water 
depth . 

FIGURE  44  shows  a  portion  of  a  surface  for  a  relative 
depth  equal  to  50,  but  the  data  are  sparse  yielding  a  very 
small  portion  of  the  desired  surface.  This  plot  illustrates 
conclusively  the  fact  that  many  more  flume  experiments  are 
necessary  and  shows  where  the  data  are  lacking. 

A  direct  comparison  of  the  derived  surfaces  for 
both  uniform  and  mixed  bed-material,  for  relative  depth 
values  of  20,  is  shown  in  FIGURE  45.  This  plot  indicates 
that  there  is  some  distinction  between  transport  of  bed- 
material  having  a  mixture  of  sizes  as  compared  to  uniform 
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material  when  the  flow  intensity  Y  is  identical.  The  data 
are,  however,  too  sparse  to  make  any  definite  conclusions. 

It  may  be  in  order,  at  this  stage,  to  compare  the 
derived  dimensionless  plots  to  commonly  used  bed-load  trans¬ 
port  formulas.  The  Einstein  bed-load  function  (Einstein 

1950)  is  the  most  popular  in  North  America,  and  relates  the 
bed-load  function  $  to  the  intensity  of  shear  iJj  : 


$  =  F  (40 


which  is : 


4  (4 


•)1/2  <-± 


D 


35 


n1/2  1  =  is _ 

3  *  F  Y  R/  S 


(4.8) 


ys  "s  gDJ  *  '  b 

where:  =  hydraulic  radius  with  respect  to  particle 

F  =  J2./ 


36v 


2/3  +  gD* (S-l) 


/  36vz 

gDJ (s-i) 


S  =  specific  gravity  of  particles . 
s 

Chien  (1954)  has  also  shown  that  the  Meyer-Peter  and  Muller 
bed-load  transport  equation  is  almost  identical  to  that  of 


Einstein  which  should  not  be  surprising  since  both  are  de¬ 
rived  from  essentially  the  same  flume  data. 

The  4>  function  of  Einstein  (1950)  is  the  reciprocal 
of  Y  (ip  =  i)  whereas  the  bed-load  function  <J)  =  P^,  where  w 

Tj-k  /  m  , 

is  fall  velocity  (Yalin  1966) .  The  ratio  w  is  a  function 
of  X  and  Y  which  means  that  Einstein's  relationship  took 
into  account  three  of  the  four  dimensionless  parameters,  the 


0ffgct  of  relative  depth  was  not  considered.  The  curves 
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shown  in  FIGURES  35  and  36  indicate  that  the  effect  of 
relative  depth  Z  is  significant;  the  curves  in  FIGURE  35 
define  a  surface  for  Z  =  Rb/D  =  15  which  is  somewhat  higher 
than  the  surface  for  Z  =  20  in  FIGURE  36.  Ignoring  the  Z 
parameter  and  plotting  P  versus  Y,  for  the  uniform  material, 
yields  a  graph  with  a  fair  degree  of  scatter  (FIGURE  46) . 

No  attempt  was  made  to  obtain  the  best  fit  line  through  the 
plotted  points  as  the  derived  equation  would  not  include  all  the 
significant  variables.  A  comparable  degree  of  scatter  is 
to  be  found  on  a  plot  of  Einstein's  bed-load  function 
(FIGURE  47)  as  produced  by  Laursen  and  Rouse  (Brown  1949). 
Instead  of  attempting  to  arrive  at  one  functional  relation¬ 
ship  relating  bed-load  transport  to  flow  intensity  it  would 

be  more  realistic  to  develop  design  charts  for  P  versus  Y 

X2 

taking  into  account  Z  and  the  effect  of  particle  size  (  /  )  . 

A  program  to  investigate  similar  parameters  is  underway  at 
the  University  of  Alberta  (Cooper  and  Peterson  1968,  Cooper 
1970)  using  all  available  sediment  transport  data  (sand  and 
gravel)  and  has  shown  reasonable  success  to  date. 

4 . 4  DESIGN  CURVES  FOR  SEDIMENT  TRANSPORT  OF  COARSE  MATERIAL 
Design  engineers,  in  general,  prefer  to  use  design 
curves  that  are  not  dimensionless.  The  dimensionless  plots 
discussed  previously  are  assumed  to  apply  to  all  types  of 
materials  as  well  as  fluids,  however,  if  only  coarse  mat¬ 
erial  under  the  action  of  water  is  considered  the  plots  ean 


, 


' 


4 


62 


be  simplified. 

Colby  (1964  a,b)  briefly  points  out  that  the  bed- 
material  discharge  in  natural  rivers  can  be  related  to  four 
common  measures : 

1.  Total  shear  on  the  stream  bed,  to  =  yRS 

2.  Mean  velocity,  Vm 

3.  Shear  velocity,  v*  =  /gRS ,  and 

4.  Stream  power,  tq  Vm. 

Colby  compares  these  various  measures  and  concludes  that  the 
relationship  of  bed-material  discharge  to  mean  velocity  is 
the  most  convenient  to  apply.  Using  the  mean  velocity  does 
not  require  a  knowledge  of  the  energy  slope  which  is  gen¬ 
erally  difficult  to  obtain  with  high  accuracy. 

With  this  in  mind  the  flume  data  from  the  previously 
quoted  investigators  was  used  to  obtain  FIGURES  48  to  55  for 
uniform  materials  and  FIGURES  57  to  60  for  bed-material  mix¬ 
tures.  Each  figure  shows  the  unit  bed-load  discharge  qs 
varying  with  the  mean  velocity  Vm  for  one  representative 
particle  size  and  a  range  of  relative  depths.  The  plots 
indicate  that  for  a  constant  velocity  the  bed-load  transport 
will  decrease  as  the  water  depth  (hydraulic  radius)  increases. 

FIGURE  56  shows  a  partial  three-dimensional  sur¬ 
face  with  the  relative  depth  R^/D  equal  to  15  and  the  uni¬ 
form  particle  size  being  the  third  dimension.  More  data 
are  required  to  give  a  complete  surface.  This  plot,  however. 
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indicates  that  for  a  constant  velocity  and  relative 
depth  the  amount  of  uniform  material  transported  is  depend¬ 
ent  on  the  size  of  the  material;  the  smaller  the  particles 
the  larger  the  volume  transported  in  a  definite  time  period. 
This  would  mean  that  the  transport  capacity  of  a  stream  is 
reduced  if  coarse  material  is  dumped  into  it,  even  if  all 
the  dumped  material  is  above  the  threshold  state. 

In  order  for  the  flume  data  to  be  useful  in 
estimating  bed-load  discharge  in  rivers  similar  plots  would 
have  to  be  derived  for  bed-material  mixtures.  However,  as 
shown  in  FIGURES  57  'to  60  the  range  of  data  is  insufficient; 
the  range  of  particle  sizes  used  in  the  experiments  varied 
from  2.0  to  4.4  mm. 

At  the  moment  design  curves  applicable  to  natural 
rivers  cannot  be  derived  due  to  a  lack  of  adequate  data. 
Laboratory  investigations  using  bed-material  mixtures  as 
well  as  field  investigations  covering  a  wide  range  of 
representative  particle  sizes  and  relative  depths  are  a 
necessity. 

4.5  COMPARISON  OF  EXISTING  BED-LOAD  FORMULAS  TO 
FIELD  MEASUREMENTS 

To  date  there  are  very  few  data  relating  the  bed¬ 
load  transport  of  a  coarse-bed  river  to  its  hydraulic 
characteristics.  The  procurement  of  this  data  is  extremely 
difficult  as  large  and  cumbersome  bed-load  traps  are  required 
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to  catch  moving  stones.  Recently,  however,  bar  movements 
on  the  North  Saskatchewan  River  at  Drayton  Valley  were  re¬ 
corded  (Galay  1967a)  from  which  estimates  of  bed-load  trans¬ 
port  can  be  made.  Also,  the  bed-load  transport  in  the  Elbow 
in  Alberta  was  studied  in  some  detail  (Hollingshead 
1968  a,b) .  These  cases  will  now  be  discussed  in  more  detail. 
4.5(a)  North  Saskatchewan  River  at  Drayton  Valley 

During  the  summer  of  1965  successive  soundings  of 
the  river-bed  yielded  the  bed  forms  shown  in  FIGURE  61  and 
PHOTOGRAPH  10.  Soundings  along  the  west  line  indicate  that 
several  large  bars  shifted  downstream  by  several  hundred 
feet  during  a  four-day  sounding  time  interval  (June  15  to 
June  19).  The  flow  conditions  during  these  soundings  were: 

Discharge-25,000  to  43,000  cfs,  (average  35,000  cfs) 

Mean  depth  -  10  feet 

Water  surface  width  -  775  feet 

Slope  -  0.0015 

Representative  bed-material  size  -  1.1  inches 
Considering  a  one-foot  strip  the  following  volumes 
of  coarse-bed  material  were  moved  during  the  four  day  period: 
Upstream  bar  from  x-sec.  89  -  690  cu.  ft. 

Downstream  bar  from  x-sec.  89  -  1330  cu.  ft. 

Mean  volume  moved  -  approximately  1,000  cu.  ft. 

This  volume  corresponds  to  1000  x  105  -  105,000  lbs. 
of  dry  bed-material.  The  steady  rate  of  bed-material  trans¬ 
port  required  to  supply  this  amount  of  material  over  a  four 
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day  period  would  be: 


105,000 


=  0.31  lb. /ft.  sec. 


4x24x60x60 

A  comparison  of  this  transport  rate,  which  is  in 
itself  relatively  approximate,  to  the  rate  obtained  by 
various  bed-load  formulas  is  shown  below: 


NORTH  SASKATCHEWAN  RIVER 
BED-LOAD  TRANSPORT  RATES  (LB ./FT . SEC . ) 


(x-sec.  89,  Q  = 

35,000  cfs) 

BAR  MOVEMENT 

COOPER  & 

MEYER-PETER 

TECHNIQUE 

BLENCH 

PETERSON 

& 

MULLER 

0.31 

0.48 

0.28 

0.41 

The 

Cooper-Peterson 

(1968)  computation 

yields  a 

result  relatively  close  to  that  computed  from  the  bar  move¬ 
ment,  however,  it  is  difficult  to  know  whether  the  bar  trans¬ 
port  rate  is  representative  for  the  cross-section.  Appendix 
V  shows  the  details  of  the  bed-load  computation  using  Blench ' s 
regime  equations . 

4.5(b)  Elbow  River  at  Bragg  Creek,  Alberta 

A  number  of  bed-load  discharge  measurements  were 
taken  during  the  summer  of  1967  using  basket  and  VUV  samplers 
(Hollingshead ,  1968  a,b) .  The  results  of  the  samplings  along 
with  the  hydraulic  characteristics  of  the  river  are  shown 
in  TABLE  7. 

The  bed-load  transport  qs  is  plotted  against  the 
mean  velocity  Vm  with  the  relative  depth  as  the  third 
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variable  (FIGURE  62).  Although  the  data  are  sparse  a  line 
having  d*/D  equal  to  30  was  sketched  in.  This  line  appears 
to  have  a  similar  slope  and  curvature  to  the  curves  based  on 
coarse  bed-load  transport  in  flumes.  It  is,  however, 
impossible  to  add  more  relative  depth  lines  to  this  graph  - 
more  data  is  required. 

A  comparison  of  the  measured  transport  rate  in  the 
Elbow  River  to  the  rates  computed  from  various  bed- load 
formulas  is  shown  in  TABLE  8. 

It  is  difficult  to  compare  the  computed  transport 
rates  with  the  measured  values  by  viewing  the  table.  There 
fore,  the  computed  transport  rates  were  plotted  versus  the 
mean  flow  velocity  and  curves  having  a  relative  depth  d*/D‘ 
equal  to  30  were  determined  for  each  bed- load  transport 
formulas.  These  curves  were  replotted  in  FIGURE  63  along 
with  the  measured  transport  rate.  It  is  apparent  that  there 
is  no  single  transport  formulas  that  can  be  recommended  over 
all  the  rest;  the  curves  are  based  on  flume  experiements 
using  uniform  sand  as  the  bed  material.  The  charts  deve¬ 
loped  by  Cooper  and  Peterson  (1968)  are,  however,  the 
easiest  to  use;  one  of  the  charts  is  reproduced  in 
FIGURE  64.  The  extrapolation  of  the  "sand-bed  formulas" 
to  coarse-bed  rivers  is,  however,  not  the  sole  reason  for  wide 
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TABLE  7 

BED-LOAD  TRANSPORT  IN  THE  ELBOW 
RIVER  AT  BRAGG  CREEK/  ALBERTA 


(at 

Cableway , 

Sta.  12+18 

D/S) 

°50  =  °* 

083  feet 

DATE 

Q 

cf  s 

Dw 

ft. 

d*=A/b 
ft.  w 

S 

Vm 

ft ./sec . 

qs 

lb. /ft. 
sec . 

a*/D 

1967 
June  1 

3850 

160 

2.94 

0.00745 

8.40 

0.855 

35 

June  2 

2900 

150 

2.68 

0.00745 

7.40 

0.540 

32 

June  17 

1900 

136 

2.30 

0.00745 

6.20 

0.304 

28 

June  19 

1900 

136 

2.30 

0.00745 

6.20 

0.310 

28 

June  20 

1630 

131 

2.20 

0.00745 

5.80 

0.092 

26 

June  21 

1560 

130 

2.17 

0.00745 

5.66 

0.197 

26 

June  22 

1550 

130 

2.16 

0.00745 

5.65 

0.086 

26 

June  22 

1470 

128 

2.12 

0.00745 

5.50 

0.175 

26 

June  23 

1300 

124 

2.03 

0.00745 

5.20 

- 

24 

1968 
June  8 

1490 

129 

2.13 

0.00745 

5.45 

0.196 

26 

June  9 

1400 

127 

2.09 

0.00745 

5.30 

0.014 

25 

June  10 

1370 

126 

2.05 

0.00745 

5.26 

0.023 

25 

. 
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TABLE  3 

BED-LOAD  TRANSPORT  COMPARISONS 
ELBOW  RIVER 

(Sta.  12+18) 

(lb. /ft.  sec.) 


Q 

Measured 

Meyer-Peter 

Blench 

Cooper 

cf  s 

3S 

&  Peterson 

^s 

1,000 

- 

- 

0.080 

- 

1,300 

- 

0.177 

0.104 

- 

1,470 

0.125 

0.330 

0.122 

0.011 

1,550 

0.061 

0.424 

0.126 

0.015 

1,630 

0.066 

0.540 

0.132 

0.027 

1,900 

0.217 

0.875 

0.157 

0.105 

2,900 

0.387 

2.270 

0.240 

1.090 

3,850 

0.610 

3-800 

0.314 

2.250 

5,000 

— 

5,760 

0.403 

3,570 
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divergence  of  the  curves  in  FIGURE  63.  Using  these  formulas 
on  sand-bed  rivers  results  in  curves  that  also  diverge 
markedly  from  each  other  (FIGURE  65).  Only  Blench' s  formula 
makes  an  allowance  for  the  geometry  of  the  channel  in  plan 
(fg). 

4.6  CONCLUSIONS 

The  solutions  to  important  questions  dealing  with 
. reservoir  sedimentation  and  channel  degradation  on  coarse- 
bed  rivers  depend  upon  an  accurate  assessment  of  the  bed¬ 
load  transport.  However,  dimensionless  plots  based  on 
coarse  bed-load  transport  in  flumes  indicate  that  there  is 
insufficient  data  to  develop  relationships  that  would  be 
applicable  to  natural  rivers.  The  wisdom  in  attempting  to 
arrive  at  one  relationship  that  could  be  used  in  computing 
bed-load  transport  is  questioned;  dimensional  analysis  has 
shown  that  the  transport  is  related  to  at  least  three 
dimensionless  parameters. 

A  comparison  of  computed  bed-load  transport  using 
several  common  formulas  and  charts  (Meyer-Peter  and  Muller, 
Blench,  Cooper-Peterson)  to  the  actual  measured  transport 
in  the  Elbow  River  indicates  that  no  one  formula  can  be  re¬ 
commended  over  all  others.  The  Meyer-Peter  and  Muller  formula 
is  generally  recommended  for  coarse-bed  rivers;  however, 
applying  it  to  the  Elbow  River  resulted  in  transport  values 


6  or  7  times  the  actual  measured  values.  Of  the  formulas 
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and  relationships  examined  it  was  found  to  be  the  least 

reliable.  The  relationships  developed  by  Cooper  and  Peterson 

(1968)  were  the  easiest  to  use.  Most  of  the  formulas  were 

based  on  data  from  sand-bed  channels  and  flumes  and  yield 
♦ 

widely  varying  answers  even  for  channels  in  sand.  Computing 
the  bed-load  transport,  in  the  North  Saskatchewan  River, 
from  shifting  gravel  bars  yielded  transport  rates  that 
compared  favourably  with  computations  based  on  various 
formulas . 
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4.7  RECOMMENDATIONS 

Flume  study  data  for  coarse  bed-material  mixtures 
are  virtually  non-existent.  This  data  would  be  useful  in 
understanding  the  relationship  between  transport  rate  and 
significant  variables.  The  tests  should  utilize  a  well 
graded  mixture  having  a  representative  diameter  of  at  least 
20  mm.  and  should  be  conducted  over  a  wide  range  of  depths. 
With  conventional  flumes  a  flow  depth  of  three  feet  would  be 
the  maximum  possible,  yielding  a  relative  depth  of  40. 

Larger  ratios  of  relative  depth  would  be  desirable.  Mix¬ 
tures  of  light  weight  materials  would  also  yield  information 
on  the  importance  of  the  density  of  the  transported  material. 

The  collection  of  data  from  mobile  bed  rivers  should 
coincide  with  and  complement  the  flume  studies  recommended 
above.  More  transport  measurements  should  be  conducted  on 
rivers  similar  in  size  to  the  Elbow  River  (relative  depths 
varying  from  20  to  50).  Investigations  on  larger  rivers 
such  as  the  North  Saskatchewan  and  the  Athabasca  would 
yield  information  on  rivers  having  relative  depths  up  to 
150.  The  use  of. tracer  techniques  may  be  necessary  in  these 
larger  rivers  as  the  trapping  of  coarse  material  would  be 
difficult.  Continuous  sounding  of  moving  dunes  and  major 
bars  would  yield  valuable  information  on  the  rate  of  trans¬ 
port  as  well  as  the  form  of  bed  roughness. 
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The  plotting  of  data  in  a  direct  manner  such  as  trans¬ 
port  rate  versus  mean  velocity  with  relative  depth  as  the  third 
variable  appears  promising.  Plots  of  this  type  would  be  readily 
used  by  practicing  engineers  as  they  are  direct  and  uncomplicat¬ 
ed.  Effort  should  be  directed  towards  establishing  the  most 
practical  plots  since  there  is  much  "discrepancy  among  existing 
complicated  formulas. 
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CHAPTER  V 


RESISTANCE  TO  FLOW  IN 
COARSE-BED  RIVERS  AND  FLUMES 

5 . 1  INTRODUCTION 

In  alluvial  channels  the  bed  forms  generated  during 
the  passage  of  high  flows  are  complex  and  are  continuously 
undergoing  change.  Laboratory  experiments  (Simons  and 
Richardson,  1961;  Guy,  Simons  and  Richardson,  1966)  have 
established  that  the  bed  forms  in  sand-bed  rivers  pass 
through  the  following  regimes  (see  FIGURE  66) : 

(a)  Plane  bed  with  no  sediment  motion 

(b)  Ripples 

(c)  Ripples  on  dunes 

(d)  Dunes 

(e)  Washed  out  dunes 

(f)  Plane  bed  with  motion 

(g)  Standing  waves 

(h)  Anti-dunes 

(i)  Chutes  and  pools. 

Needless  to  say,  these  changes  in  the  bed  forms 
would  produce  appreciable  changes  in  flow  resistance. 

All  the  experiments  and  observations  upon  which 
this  bed-form  classification  is  based  involve  sand  as  the 
alluvial  material.  The  extrapolation  of  this  classification 
to  coarse-bed  alluvial  channels  may,  at  the  moment,  be 
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unwarranted.  Echo-soundings  on  the  North  Saskatchewan  River 
(Galay,  1967b)  have  shown  that  dunes  in  coarse-bed  rivers 
behave  differently  from  dunes  in  sand-bed  rivers .  The  main 
differences  noted  were: 

(1)  The  wave  length  decreased  with  an  increase 
in  flow  depth  and  velocity  which  is  the  opposite 
for  dunes  in  sand-bed  channels  (see  PHOTOGRAPH  10) ; 

(2)  The  dune  amplitude  did  not  appear  to  be 
controlled  by  the  depth  of  flow.  The  amplitude 
of  dunes  in  sand-bed  channels  appears  to  increase 
directly  with  depth  (Nordin  and  Algert,  1965) . 

A  similar  phenomenon  was  noted  by  Theil  (1932)  who 
described  bed  forms  exposed  on  the  bed  of  a  stream  that 
partially  drained  a  lake  in  Minnesota.  The  wave  lengths 
varied  from  25  to  60  feet  and  were  found  to  decrease  in  length 
as  the  velocity  increased.  Theil,  however,  erroneously 
assumed  that  the  velocity  decreased  in  a  downstream  direction. 

Flume  experiments  by  many  researchers  have  also 
verified  that  coarse  bed-material  will  not  form  into  ripples. 
Simons  and  Richardson  (1966)  have  produced  a  graph,  shown  in 
FIGURE  67,  which  shows  that  ripples  do  not  develop  if  the 
bed-material  is  larger  than  0.65  mm.  Photographs  of  water 
surface  features  suggesting  the  presence  of  antidunes  in  a 
coarse-bed  river  in  New  Zealand  have  also  been  presented  by 
Thompson  (1963) .  The  range  of  bed  forms  that  may  be  generated 
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are  open  to  investigation. 

The  fact  remains  that  the  hydraulic  engineer  is 
called  upon  to  evaluate  the  roughness  or  the  resistance  to 
flow  of  a  coarse  river-bed  in  either  an  immobile  or  mobile 
state.  This  chapter  will  assess  the  resistance  of  coarse 
river  beds  in  two  distinct  states  -  when  the  river-bed  is 
immobile  and  can  be  considered  as  a  rigid  bed  and  when  the 
bed  is  fully  mobile. 

5 . 2  IMMOBILE  RIVER-BED 
5.2(a)  Dimensional  Analysis 

The  resistahce  to  flow  over  a  rigid  coarse  bed 
depends  upon  the  flow  conditions,  the  channel  geometry,  the 

i 

properties  of  the  fluid  and  the  characteristics  of  the  pro¬ 
trusions  on  the  boundary.  The  variables  involved  can  be 
assembled  in  the  following  manner  (see  FIGURE  68) : 

Flow  conditions 

Vm  =  mean  flow  velocity 
d*  =  mean  flow  depth 
S  =  slope  of  energy  gradient 
Cw  =  concentration  of  suspended  particles 
Channel  geometry 

b  =  width  of  channel 

fc  =  factor  defining  cross-sectional  shape  of  channel 
fg  =  factor  defining  plan  geometry  of  channel 


' 


76 


Fluid  properties 

p  =  density  of  fluid 
y  =  dynamic  viscosity  of  fluid 
Protrusion  characteristics 

k  =  height  of  protrusion  (usually  k5Q) 

=  shape  factor  of  protrusions 
X  =  concentration  of  protrusions 
X  =  areal  pattern  of  protrusions 
and  g  =  acceleration  of  gravity 

Setting  this  out  in  equation  form: 

fi  <va*,s,cw,b,fc,fg,p,p,k,Vx'X'S'>  =  o  (5-D 

Several  of  the  above  variables  can  be  omitted  in 
order  to  simplify  the  equation.  The  data  to  be  analyzed  have 
been  derived  from  straight  flumes  or  rivers  having  a  uniform 
wide  cross-section  (width  more  than  5d*) ;  the  variables  b, 
f  ,  and  f  are  therefore  eliminated.  The  effect  of  flow  con- 

p  g 

centration  on  resistance  involving  large  protrusions  on 
the  bed  has  not  been  investigated  and  will  be  ignored  in  this 
treatment.  The  protrusion  or  particle  shape  factor  does 
not  vary  from  river  to  river  as  shown  in  FIGURE  69  and  is 
tentatively  eliminated.  The  effect  of  the  concentration  and 
distribution  of  bed  particles  on  resistance  is  difficult  to 
assess  directly  in  natural  rivers  and  will  be  tentatively 
ignored. 
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Equation  5.1  therefore  reduces  to: 

f 2  (Vm'  d*'  S'  p'  k'  g)  =  0 


(5.2) 


According  to  the  Buchingham  tt  theorem  the  above 
variables  can  be  reduced  to  four  dimensionless  parameters. 
Choosing  p,  V  and  d.  as  repeating  variables  results  in: 


m 


f 2  •  S/  /  SLl.]  =  o 


(5.3) 


'gci*  y  k  ) 

For  any  significant  flow  in  coarse-bed  rivers  the 

flow  can  generally  be  described  as  fully  rough  turbulent  flow 

o  ^7  d 

which  makes  the  Reynolds  number  - — *-  of  secondary  import¬ 
ance.  Therefore: 


V. 


(73Sr  '  s- 


■)  =  0 


(5.4) 


This  equation  can  now  be  investigated  in  the  above 

form,  however,  replacing  S  by  tq  through  the  relationship 

V  C 

t  =  Yd.S  leads  to  a  familiar  resistance  term  — 1 11  or  -7— 
o'*  v  /g 


where  C  =  Chezy’s  C: 


f4  (Vm 


V 


m 


7gd7  /to/ 


,  d*) 


(5.5) 


or  re-arrangmg: 


V 


f4  %  '  r> 


The  Froude  number  assumes  importance  only 

when  appreciable  surface  waves  develop  (Rouse,  Koloseus  and 


Davidian  1963)  and  is  eliminated  since  most  of  the  data  to 
be  analyzed  were  for  Froude  numbers  less  than  one.  This 
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results  in: 

Vm  =  fc  (dj  (5.6) 

n  5  k 

5.2(b)  Emperical  Flow  Formulas 

Instead  of  using  the  concentration  and  pattern  para¬ 
meters  explicitly , investigators  have  taken  the  protrusion 
height  k  to  be  larger  than  the  median  size  of  the  bed- 
material  and  assumed  that  this  will  account  for  these  unknown 
parameters.  Kellerhals  (1963)  arrived  at  the  following 
formula : 


6.5 


(**) 


1/4 


(5.7) 


where:  k  =  from  a  grid  sample  analyzed  by  number. 

Kellerhals'  formula  is  almost  identical  to  the  logarithmic 
flow  formula  formulated  by  Keulegan  (1938)  for  wide  channels: 
V 

—  =  5.75  log (  £*)  +  6.00  (5.8) 

v*  k 

In  this  formula  the  k  parameter  refers  to  a  uniform  size  of 
bed-material  and  is  based  on  experiments  carried  out  by 
Nikuradse  (1932)  who  coated  the  inside  of  pipes  with  uniform 
sand  and  conducted  tests  to  assess  the  effect  of  this  sand 
roughness  on  flow.  However,  it  is  apparent  that  the  pro¬ 
trusion  height  of  the  glued  sand  particles  would  be  some 
fraction  of  the  particle  diameter,  possibly  equal  to  one-half 
the  diameter  as  the  grain  contact  points  would  be  at  D/2  from 
the  boundary . 


■ 


' 


' 

■ 


■ 


79 


It  has  also  been  noted  by  several  investigators 
(Lane  and  Carlson  1954)  that  coarse  bed-material  is  generally 
oriented  into  a  "shingle  pattern"  on  a  river-bed.  PHOTOGRAPH 
11  illustrates  this  shingle  pattern  and  FIGURE  70  shows  that 
the  particles  will  orient  themselves  so  that  the  long  "a" 
axis  is  perpendicular  to  the  flow.  It  is  apparent  that  the 
protrusion  height  of  the  bed  particles  into  the  flow  may  not 
be  the  intermediate  b  axis  as  used  by  some  investigators. 

The  question  remains,  however,  what  measure  should  be  used 
as  the  protrusion  height? 

A  "bed  roughness  meter"  was  therefore  constructed 
to  directly  measure  the  protrusion  height  of  gravel  and 
cobbles  on  a  shingled  stream  bed.  PHOTOGRAPH  12  shows  the 
"roughness  meter"  in  operation  on  the  North  Saskatchewan 
River.  A  movable  arm  is  drawn  over  a  cobble  bed  with  its 
vertical  motion  recorded  on  a  moving  chart.  The  chart  and 
arm  movements  were  geared  so  that  the  resulting  trace 
dimensions  were  exactly  one-half  the  actual  bed  pattern 
dimensions.  A  sample  of  the  resulting  trace  is  shown  in 
FIGURE  71.  Measurements  of  the  three  axes  of  cobbles  were 
also  taken  at  the  same  location  as  the  protrusion  height 
measurements . 

From  the  relatively  few  measurements  of  protrusion 
height  taken  in  1966  and  1967  a  plot  of  median  protrusion 
height  versus  median  intermediate  (b)  and  median  minor  (c) 


, 
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axis  was  obtained  (FIGURE  72).  A  tentative  conclusion  from 
this  plot  would  be: 

Protrusion  height  k  =  1.0c 

or  Protrusion  height  k  =  0.67b  or  0.67D^Q  (by  number) 

This  information  is  now  used  to  obtain  flow  formulas 
for  coarse-bed  rivers  which  would  be  based  on  the  projection 
height  of  the  bed-material  into  the  flow. 

,  Data  from  coarse-bed  rivers  having  distinctly 
immobile  beds  was  gathered  and  tabulated  in  TABLE  A-7, 

T  7—. 

APPENDIX  6.  A  plot  of  —  versus  d*/k  results  in  the 


equation  (FIGURE  73): 


3.0 


(5.9) 


The  plot  exhibits  a  fair  degree  of  scatter;  the  points  that 
deviate  markedly  are  from  investigations  reported  by  Barnes 
(1967)  and  Alberta  Water  Resources.  The  reason  for  the 
scatter  may  be  the  variation  in  residual  bed  forms  or  bars 
that  exist  when  the  bed  material  stops  moving. 

The  power  and  constant  in  equation  5.9  do  not 
agree  with  the  values  in  Kellerhals'  equation  5.7,  but  the 
protrusion  height  may  not  be  directly  related  to  the  D90 
by  number.  For  this  reason  flume  tests  on  naturally 
sorted  gravel  beds  were  carried  out.  A  description  of  the 
tests  along  with  their  results  is  presented  in  APPENDIX  7. 

The  flume  tests  with  the  natural  gravel  bed  plot  very  closely 
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to  the  derived  line  (equation  5.9)  as  shown  in  FIGURE  74. 

The  data,  however,  covers  a  limited  range  of  relative  depths. 

Further  flume  tests  were  carried  out  having  an 
artificial  cemented  gravel  bed  similar  to  Nikuradse's  sand 
pipe  tests  (see  PHOTOGRAPH  13) .  It  is  interesting  to  note 
that  the  results  from  the  cemented  gravel  bed  indicate  higher 
roughness  values  than  the  natural  sorted  gravel  bed  even 
though  the  median  projection  height  was  higher  for  the  natural' 
bed  (FIGURE  74) .  The  orientation  of  the  gravel  by  flowing 
water,  therefore,  has  a  significant  effect  on  the  resistance 
to  flow.  Therefore, ,the  extension  of  Nikuradse's  pipe  tests 
with  cemented  sand  grains  to  natural  rivers  would  be  highly 
questionable . 

The  power  formula  (equation  5.9)  can  be  replaced  by 
a  logarithmic  formula: 

—  =  8.0  log(^)  +  1.0  (5.10) 

v+  K 

* 

This  equation  is  of  the  form  proposed  by  Keulegan  (1938)  and 

Robinson  and  Albertson  (1952) : 

Vm  2.30  -|  ,d.> 

_  =  log  (j.*)  +  Cl 

where:  K  =  vonKarman  turbulence  coefficient 

constant 

The  logarithmic  formula,  5.10,  would  have  a  turbulence  co¬ 
efficient  K  equal  to  0.29  instead  of  the  commonly  accepted 
0.40.  The  turbulence  coefficient  does  range  from  0.20  to 
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to  0.40  for  sediment  laden  flows  (Chien  1956)  and  may  vary 
in  rigid  boundary  flow  due  to  such  -  factors  as  secondary 
circulation  and  boundary  roughness.  Equation  5.10  along  with 
Keulegan's  equation  5.8  are  shown  in  FIGURE  75.  The  use  of 
Keulegan's  equation  implies  the  acceptance  of  a  coefficient 
and  constant  which  are  based  on  Nikuradse's  pipe  tests  and 
results  in  a  computation  of  an  "equivalent  sand  grain  size". 

A  more  realistic  approach  is,  however,  suggested  here  -  the 
protrusion  height  that  is  actually  measured  should  be  used 
in  a  flow  formula  with  the  turbulence  coefficient  and  con¬ 
stant  adjusted  to  fit  the  observed  flow  conditions.  This 
approach,  results  in  the  use  of  actual  protrusion  heights; 
flume  experiments  utilizing  artificial  roughness  elements 
such  as  bars  (Powell  1946) ,  baffles  (Sayre  and  Albertson 
1963),  and  cubes  (O'Loughlin  and  MacDonald  1964)  all  used  the 
protrusion  heights  in  subsequent  analysis. 

A  number  of  studies  have  recently  been  conducted  on 
the  resistance  to  flow  of  streams  having  "large  bed  elements", 
(Judd  and  Peterson  1969,  Herbich  and  Shulits  1964).  Uniform 
two-dimensional  flow  was  not  always  present  for  the  tabulated 
flows  as  the  bed  elements  were  so  large  that  they  represented 
changes  in  cross-section  rather  than  surface  texture.  For 
this  reason  the  data  were  not  included  in  the  analysis  dealing 
with  immobile  channels. 
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5.2(c)  Manning's  Roughness  Coefficient 

In  North  America  the  use  of  the  Manning  equation  to 
compute  flow  in  natural  channels  is  widespread.  Many 
engineers  can  associate  the  physical  appearance  of  a  stream- 
bed  with  a  Manning’s  roughness  coefficient  n.  The  roughness 
coefficient  could  also  be  estimated  by  using  the  Strickler 
formula : 

,1/6 


n  =  0.032  D' 


with  D  being  D^q  and  measured  in  feet. 


(5.11) 


A  better  relationship  for  estimating  Manning's  n 
which  would  include  the  effects  of  flow  depth  and  make  use 
of  the  protrusion  height  of  the  bed-material  can  be  obtained 
by  using  the  Manning  equation  with  the  previously  derived 


power  flow  equation  (Equation  5.9) 


vm  =  d*2/3  s1/2 


n 


(5.12) 


and  using: 


=  3.0  (^) 


0.45 


yields : 


n 


1/6 


0.088 

(d  )0*45 
a*/  > 


(5.13) 


where  k  =  0.67  Dr_  in  feet. 

oU 


The  data  from  APPENDIX  6  have  been  plotted  with 


n/ 


1/6 


,  versus  d*  ,  in  FIGURE  76.  There  is  again  a  fair 

a*  /k 


degree  of  scatter,  but  the  above  derived  equation  does  appear 
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to  pass  through  the  points  reasonably  well;  this  relation¬ 
ship  would,  therefore,  yield  a  reasonable  estimate  of  the 
Manning , coefficient  for  a  reasonably  straight  immobile  coarse 
river-bed.  Its  use,  however,  should  be  limited  to  flows  be¬ 
low  those  that  would  cause  most  of  a  rivers 1  bed  to  be  in 
motion;  this  aspect  of  resistance  to  flow  with  a  mobile 
stream  bed  becomes  fairly  complex. 

5.3  MOBILE  RIVER-BEDS 

5.3  ( a )  Dimensional  Analysis 

As  in  the  chapter  dealing  with  the  transport  of 
coarse  bed-material  ('CHAPTER  IV)  the  resistance  to  flow  is 
related  to  the  flow  conditions,  the  channel  geometry  and  the 
properties  of  the  fluid  and  bed-material.  The  mean  flow 
velocity  is  therefore  a  function  of: 

Vm  =  f  ^  (d*  ,  S  ,C^  ,b  ,  f  c  ,  f  ,  p  ,  y  ,D ,  pg  ,  ,g)  (5.14) 

The  following  variables  C  ,  b,  f  ,  f  ,  a,  ,  and  « 

r  w  c  g  b  b 

are  now  deleted  for  the  reasons  stated  in  CHAPTER  IV,  sect¬ 
ion  4.2. 

This  results  in: 

Vm  =  f2  (d*,  S,  p,  y ,  D,  ps,  g)  (5.15) 

It  is  possible  to  express  certain  variables  in  terms 

of  others;  g  and  S  can  be  replace  by  y’  and  v*  using  the 

s 

relationships : 

y's  =  9  <ps-p) 

and 

v*  =  v/gd*S 


' 


85 


This  changes  equation  5.15  to: 

Vm  =  f2  (d*,  v*,  p,  ]i,  D ,  ps,  y')  (5.16) 

Using  the  Buckingham  tt  theorem  with  p,  D,  and  v*  as 


repeating  variables  results  in: 


Vm  _  _ 

n  =  -f3 


( 


PDV, 

y 


pv, 


ps 


7  Yg  D/  D  '  P 


) 


(5.17) 


Assuming  again  that  the  relative  density  parameter 
is  of  minor  importance,  that  d*  and  v*  are  replace  by  and 
v^  respectively  and  that  the  particle  Reynolds  number  can  be 
manipulated  leads  to: 

(5.18) 


Vm  v'  7  2 

^-  =  f4  (  sD  Pv*b  Rh  > 
v*b  H  o'  -  1  '  n 

pv  ys  D 


5.3(b)  Analysis  of  Flume  Data 

2 

Plots  of  Vm/  versus  the  mobility  number  for 

v*b  *  Yg  D 

various  values  of  Rb/p  snd  for  several  sizes  of  material  are 

shown  in  FIGURES  77  to  81.  These  plots  all  pertain  to  uniform 

material.  A  subsequent  plot  of  Vm/  versus  the  mobility  num- 

v* 

ber  for  a  constant  value  of  Rb/D  =  15  (FIGURE  82)  indicates 
that  the  size  of  material  has  little  influence  on  resistance. 
This  was,  however,  not  the  case  for  sediment  transport.  The 
amount  of  uniform  material  in  transport,  with  Rb/p  an<^  t he 
mobility  number  constant,  varies  with  the  size  of  material 
(see  FIGURE  35) ,  but  the  resistance  to  flow  is  not  greatly 
affected.  The  resistance  function  Vm y  is  reasonably  con- 

stant  at  a  value  of  10.5  for  mobility  numbers  ranging  from 


0.10  to  0.20. 
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,  The  data  for  uniform  materials  as  obtained  by 
Bogardi  and  Yen  are  shown  in  FIGURES  83  to  85,  but  the  range 
of  flow  depths  used  was  very  limited.  The  same  can  be  said 
of  the  data  from  flume  tests  using  bed-material  mixtures 
(FIGURES  86  to  89) ;  the  data  are  very  sparse  and  relative 
depth  curves  are  difficult  to  determine.  However,  the  re¬ 
sistance  for  bed-material  mixtures  is  higher  than  for  uniform 
materials,  for  the  same  mobility  number  (FIGURE  90).  These 
findings  are  contrary  to  the  conclusions  obtained  by 
Daranandana  (1962)  who  studied  the  effect  of  the  gradation  of 
sands  on  flow  phenomena.  He  found  that  the  resistance  to  flow 
was  considerably  greater  for  a  uniform  sand  as  compared  to  a 
poorly  sorted  sand  when  the  bed  forms  were  in  the  dune  phase. 
This  would  imply  that  conclusions  derived  from  flume  experi¬ 
ments  using  sands  cannot  be  extrapolated  to  coarse-bed  mat¬ 
erials  . 

5.3  ( c )  Analysis  of  Field  Data 

To  date  there  exists  little  information  on  the  re¬ 
sistance  to  flow  for  coarse-bed  rivers  under  relatively  large 
flows.  Data  from  the  Elbow  River  (Hollingshead  1968)  and 
the  North  Saskatchewan  River  (Galay  1967a)  are,  however, 
available  for  analysis. 

A  plot  of  Vm/  versus  the  mobility  number  for  the 

v*b 

two  rivers  is  shown  in  FIGURE  91,  with  the  data  presented  in 
TABLE  9.  A  curve  for  a  rigid  coarse-bed  channel,  based  on 
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equation  5.9,  is  also  shown  on  the  plot.  The  plotted  points 
indicate  that  the  resistance  increased  sharply  when  the  mat¬ 
erial  started  to  move  and  then  decreased  gradually  as  the 
mobility  number  increased.  This  would  suggest  that  the  bed- 
material  initially  formed  into  dunes  which  were  subsequently 
washed  out  as  the  flow  intensity  increased.  The  North  Sask¬ 
atchewan  River  data  falls  into  the  same  region  as  the  Elbow 
River  data  but  covers  a  smaller  range  of  mobility  numbers. 

Due  to  the  sparseness  of  the  data  it  is  difficult 
to  recommend  a  formula  or  a  procedure  to  obtain  the  resist¬ 
ance  to  flow  in  coarse  mobile  bed  rivers. 


5.4  CONCLUSIONS 


The  generation  of  complex  bed  forms  in  coarse-bed 
rivers  makes  the  prediction  of  the  resistance  to  flow  exceed¬ 
ingly  difficult.  The  problem  was  dealt  with  by  separating 
the  flow  condition  into  two  distinct  categories,  the  first 
being  flow  over  a  rigid  immobile  bed  while  the  second  con¬ 
sisting  of  flow  over  a  highly  mobile  bed. 

Dimensional  analysis  was  used  to  reduce  the  number 
of  variables  for  the  immobile  bed  condition  and  subsequent 

plotting  resulted  in  the  following  equation: 

Vm 


=  3.0  (f*)°‘45 

v*  k 


(5.9) 

The  protrusion  height  k  in  the  above  equation  is  the  actual 


protrusion  of  stones  into  the  flow,  not  the  median  of  the 
intermediate  axis.  This  value  was  determined  with  the  use 


, 
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of  a  "roughness  meter"  which  provided  a  graphical  trace  of 
the  river-bed.  The  projection  height  was  correlated  to  the 
intermediate  diameter  of  the  surface  stones  yielding: 

Projection  height,  k  =  0.67D5Q  (by  number) 


Flume  tests  were  also  conducted  to  test  equation  5.9  and  the 


concept  of  protrusion  height.  The  tests  consisted  of  water 
flowing  over  a  naturally  sorted  gravel  bed  and  the  results, 
covering  a  small  range  of  relative  depths,  agreed  well  with 
the  derived  equation. 

Equation  5.9  was  also  combined  with  the  Manning 
equation  to  yield: 


h 

<VA 


0.088 


0.45 


which  is  useful  for  estimation  of  Manning's  roughness  co¬ 
efficient. 

Definite  conclusions  regarding  resistance  to  flow 
with  highly  mobile  beds  cannot  be  made  at  this  date.  The 
data  are  too  sparse.  Flume  data,  for  uniform  materials, 
tentatively  indicate  that  the  size  of  material  in  motion  has 
little  influence  on  resistance.  The  amount  of  material  in 
transport,  with  relative  depth  and  mobility  number  constant, 
varies  with  the  size  of  material  but  the  resistance  to  flow 
is  not  greatly  affected.  A  comparison  of  uniform  materials 


and  mixtures  also  indicates  that  the  resistance  to  flow  is 


greater  for  mixtures  at  comparative  flow  conditions,  more 


. 
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flume  data  are  required  to  verify  these  findings. 

Limited  field  data  suggest  that  the  bed-material, 

once  in  motion,  forms  into  dunes  and  then  into  a  plane  bed 

as  the  flow  intensity  increases.  The  resistance  function 

Vm/  decreases  by  a  factor  of  two  (indicating  added  resist- 
v* 

ance)  as  soon  as  the  material  is  in  motion.  The  use  of  a 
rigid  bed  flow  formula  in  highly  active  coarse-bed  rivers  can 
lead  to  large  errors  in  flow  estimation.  However,  no  formula 
or  procedure  was  recommended  for  use  under  these  types  of 
flow  conditions. 

5.5  RECOMMENDATIONS 

The  recommendations  presented  in  the  chapter  dealing 
with  sediment  transport  (CHAPTER  IV)  would  also  be  applicable 
to  this  chapter.  The  material  on  the  bed  of  a  river  is  placed 
into  motion  by  the  adjacent  flow  which  then  has  its  turbulence 
characteristics  modified  by  the  moving  bed-material;  the 
movement  of  the  bed-material,  the  resistance  to  flow  and  the 
motion  of  the  fluid  are  inter-dependent. 

The  experiments  conducted  with  flows  over  a  naturally 
sorted  rigid  gravel  bed  should  be  extended  -  data  having 
relative  depths  up  to  100  would  be  valuable  in  establishing 
a  rigid  bed  flow  formula.  The  variation  of  protrusion  height 
with  the  shape  of  stones  on  the  bed  of  an  immobile  channel 
should  also  be  investigated. 


' 
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The  range  of  bed  forms  developed  in  sand-bed  channels 
undergoing  large  changes  in  magnitude  of  flow  have  been  in¬ 
vestigated  extensively.  However,  the  type  and  sequence  of 
bed  forms  generated  in  coarse-bed  channels  has  not  received 
the  same  attention;  flume  studies  using  bed-material  mixtures 
under  a  wide  range  of  relative  depths  (up  to  50)  and  mobility 
numbers  would  be  necessary.  The  variation  in  the  resistance 
function  (  /  )  with  the  types  of  bed  forms  being  generated 

^  -k 

would  be  clarified  and  relationships  for  the  prediction  of 
flow  resistance  parameters  would  be  forthcoming. 

The  collection  of  data  from  mobile  bed  rivers,  such 
as  the  Elbow  River,  should  also  continue.  Echo  soundings 
along  with  appropriate  flow  measurement  would  yield  informat¬ 
ion  on  the  development  of  bed  forms  and  their  effects  on  flow 


resistance . 


‘ 
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CHAPTER  VI 


THE  DESIGN  OF  STABLE  CHANNELS 
IN  COARSE  MATERIALS 

6.1  INTRODUCTION 

The  design  of  river  diversions  or  irrigation  canals 
in  coarse  bed-material  has  not  received  the  attention  that 
sand-bed  channels  have  to  date.  The  design  of  a  channel 
that  will  be  comparatively  stable  is  of  major  importance  in 
the  proper  layout  of  vast  river  diversions  and  irrigation 
schemes  such  as  those  contemplated  by  the  Saskatchewan- 
Nelson  Basin  Board  (T969) . 

This  chapter  will  briefly  present  aspects  of  domi¬ 
nant  and  bankfull  discharge  and  relate  bankfull  discharge 
to  stable  widths  and  depths.  Design  equations  are  presented 
for  the  design  of  stable  widths,  depths  and  slopes.  A 
tentative  chart  is  presented  from  which  an  assessment  of 
bankfull  discharge  can  be  made  knowing  flow  conditions  at 
some  low  stage  of  flow. 

6.2  DOMINANT  AND  BANKFULL  DISCHARGE 

The  hydraulic  geometry  of  the  river  channel,  namely 
the  width,  depth  and  slope,  depend  primarily  upon  the  flow 
in  the  channel  (Lacey  1929,  Blench  1969).  In  order  to 
design  channels  it  is  necessary  to  obtain  relationships 
between  the  hydraulic  geometry  and  the  flow  characteristics. 
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The  main  question  at  this  point  is  -  what  discharge  should 
be  used  in  attempting  to  arrive  at  meaningful  relationships? 
Leopold  and  Maddock  (1953)  have  used  average  annual  discharge 
while  other  investigators  have  used  dominant  (Inglis,  1940) 
and  bankfull  discharges  (Nixon,  1948) .  Inglis  (1940) 
originally  defined  dominant  discharge  as: 

the  discharge  which  controls  the  meander 
length  and  breadth.  It  appears  to  be 
slightly  in  excess  of  bankfull  stage. 

However,  it  is  somewhat  uncertain  whether  the  discharge  that 

forms  the  channel  is  similar  to  the  discharge  that  controls 

the  meander  wave-length.  Recently,  Benson  and  Thomas  (1966) 

have  defined  the  term  as : 

the  discharge  that  over  a  long  period 
transports  the  most  sediment. 

They,  however,  defined  sediment  as  suspended  load  and  assumed 
that  the  total  load  is  proportional  to  the  suspended  load. 

The  importance  of  the  suspended  load  as  compared  to  the  bed¬ 
load  in  the  formation  of  a  channel  has  been  discussed  by 
several  investigators  (Blench  1969,  Chien  1955). 

At  the  moment  there  is  no  widely  accepted  discharge 
that  is  used  by  all  investigators  in  analyzing  natural  rivers 
or  channels,  however,  in  designing  channels  the  design  flow 
can  usually  be  taken  as  the  bankfull  flow.  An  allowance 
for  freeboard  is  then  added  to  the  corresponding  bankfull 
stage  to  ensure  against  overtopping.  The  hydraulic  geometry 
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of  coarse-bed  rivers  will  therefore  be  correlated  to  bank- 
full  discharge  in  the  following  section. 

6.3  WIDTH  AND  DEPTH  RELATIONSHIPS 

Engineers  in  India  have  long  been  aware  that  the 
width,  depth  and  slope  of  stable  irrigation  canals  are 
strongly  correlated  with  the  bankfull  discharge  and  the  type 
of  sediment  through  which  the  canal  flows.  Although  many 
equations  have  been  developed  by  many  individuals  the  most 
commonly  used  are  those  of  Lacey  (1929)  and  Blench  (1969) . 
These  equations,  pertaining  to  sand-bed  canals,  are  present¬ 
ed  in  TABLE  10. 

Regime  equations  for  coarse-bed  channels  have  been 
developed  by  Kellerhals  (1967)  and  Maddock  (1969)  and  are 
presented  in  TABLE  11.  Kellerhals  (1967)  equations  were 
developed  from  his  own  field  data  as  well  as  data  collected 
by  Lane  and  Carlson  (1953) .  Maddock  (1969)  has  developed 
two  sets  of  equations  based  on  the  San  Luis  Canal  data  of 
Lane  and  Carlson  (1953);  the  first  set  of  equations,  group 
(a)  in  TABLE  11,  represent  steep  canals  having  discharges 
larger  than  500  cfs  while  the  second  set  are  for  discharges 
less  than  500  cfs. 

The  previously  quoted  data  from  Kellerhals  (1967) 
and  Lane  and  Carlson  (1953)  have  been  combined  with  data 
from  the  North  Saskatchewan  River  as  well  as  other  Alberta 
rivers  to  obtain  the  last  set  of  equations  in  TABLE  11. 
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The  width  and  depth  equations  were  obtained  from  FIGURE  92 
(data  tabulated  in  TABLE  A-ll,  APPENDIX  8).  The  data, 
although  sparse,  do  not  scatter  excessively.  Kellerhals ' 
data  have  been  modified  slightly,  the  discharges  for  the 
Cariboo  River  at  Cariboo  Lake  and  the  Chilko  River  at  Chilko 
Lake  were  reduced  from  dominant  to  bankfull. 

An  attempt  to  relate  the  slope  to  bankfull  dis¬ 
charge  met  with  little  success  (FIGURE  93)  which  is  not 
surprising  since  many  coarse-bed  rivers  have  their  slopes 
controlled  by  bed-rock  outcrops.  The  slope  equation  pre¬ 
sented  in  TABLE  11: 


S  = 


0.051  D50  °‘9° 


0.25 


(6.3) 


was  arrived  at  by  combining  the  flow  resistance  formula  de¬ 
rived  for  immobile  channels  (equation  5.9)  with  the  width 
and  depth  equations  from  FIGURE  93. 

6.4  COMPARISON  OF  REGIME  EQUATIONS 
FOR  COARSE-BED  CHANNELS 

It  may  now  be  in  order  to  compare  the  several  sets 

of  regime  equations  for  coarse-bed  channels.  Consider  a 

hypothetical  river  having  a  relatively  straight  alignment 

with  the  following  characteristics: 

Qb  =  Qd  =  20,000  cfs 

D^  ~  =  3.0  inches 

50 


5 . 5  inches 


* 
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The  widths,  depths,  and  slopes  for  this  river  as 
computed  from  the  regime  relations  are  as  follows: 


Author 

Width (bw) 
ft. 

Depth (d* *) 
ft. 

Slope 

Kellerhals 

255 

9.5 

0.00120 

Maddock 

168 

16.8 

0.00046 

This  thesis 

282 

7.6 

0.00120 

The  numerical  values  for  widths,  depths,  and  slopes 
are  almost  identical  for  Kellerhals '  equations  and  those 
developed  in  this  thesis.  However,  the  computed  cross- 
section  as  based  on  Maddock ' s  equations  is  relatively  narrow 
and  deep  while  its  slope  is  about  half  that  computed  from 
the  other  regime  equations.  The  data  used  by  Maddock  to 
develop  his  equations  are,  however,  not  as  extensive  as 
those  used  in  this  thesis. 

6.5  NON-DIMENSIONAL  RATING  CURVE 

Relatively  few  investigators  record  bankfull  flow 

* 

as  pertinent  hydraulic  data.  However,  it  appears  that  this 
information  is  most  useful  in  comparing  river  behaviour. 

Leopold,  Wolman  and  Miller  (1964)  noted  that  rating 
curves  of  various  stations  are  very  similar  in  form  and 
could  be  plotted  in  non-dimensional  form  by  plotting  the 
ratios : 

mean  depth,  d*  versus  discharge,  Q _ 

bankfull  depth, d^  bankfull  discharge,  Qb 


. 


\  ,  0  ; 
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These  non-dimensional  relationships  were  plotted 
for  various  coarse-bed  rivers  in  Alberta  (FIGURE  94)  with 
reasonable  sucess.  The  values  used  in  the  graph  are  tabulat¬ 
ed  in  TABLE  12. 

This  graph  will  be  valuable  in  the  estimation  of 
bankfull  discharge  for  coarse-bed  rivers  in  Alberta,  and 
in  assessing  stable  cross-section  dimensions.  In  order  to 
assess  the  bankfull  discharge  it  is  necessary  to  know  the 
mean  depth  at  a  definite  discharge  as  well  as  the  mean  depth 
at  bankfull  flow.  This  would  require  stream  gauging  as  well 
as  the  survey  of  several  cross-sections. 

With  the  bankfull  discharge  estimated  the  stable 
width  of  the  river  can  be  assessed  from  FIGURE  92.  This 
value  could  be  compared  to  the  measured  width  from  field 
surveys  to  establish  whether  the  river  is  in  a  stable  state. 
This  plot  can  also  be  used  to  develop  a  rating  curve  for  a 
river  provided  that  the  bankfull  discharge  and  depth  are 
known . 

It  is  interesting  to  note  that  the  recurrence 
interval  for  bankfull  flows  in  Alberta  Rivers  varies  from 
5  years  to  15  years  which  does  not  correspond  to  the  re¬ 
currence  interval  of  1.5  years  quoted  by  Leopold,  Wolman 
and  Miller  (1964)  for  rivers  in  Eastern  United  States.  The 
higher  recurrence  interval  is  also  characteristic  of  rivers 
in  British  Columbia  as  reported  by  Kellerhals  (1963) . 


■ 
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TABLE  12 


NON-DIMENSIONAL  RATING  CURVE  DATA 
COARSE-BED  CHANNELS 


River 

b 

cf  s 

Q 

cf  s 

d* 

ft. 

Q/Qb 

d*/d 

North  Sask. 

80,000 

15.0 

__ 

at  Drayton 

( 15yrs . 

7,000 

4.3 

0.09 

0.29 

Valley 

return 

10,000 

5.2 

0.13 

0.35 

(Galay,  1967a) 

period) 

20,000 

7.4 

0.25 

0.49 

n 

30,000 

9.1 

0.38 

0.61 

40,000 

10.4 

0.50 

0.70 

60,000 

13.2 

0.75 

0 .88 

* 

80,000 

15.0 

1.00 

1.00 

Red  Deer  R. 
at  Red  Deer 

30,000 

11.7 

(Qureshi , 

(10  yr. 

1,200 

2.1 

0.04 

0.18 

1962  and 

return 

3,000 

3.5 

0.10 

0.30 

Neill,  et.al 

period) 

5,000 

4.5 

0.17 

0.39 

1964) 

8,000 

5.7 

0.27 

0.49 

A 

10,000 

6.4 

0.33 

0.55 

20,000 

8.9 

0.67 

0.76 

30,000 

11.7 

1.00 

1.00 

Highwood  R. 
(Alta . W. Res . ) 

6,000 

5.3 

0 

(4  yr. 

1,000 

2.0 

0.17 

0.38 

return 

2,000 

3 . 0 

0.33 

0.57 

period) 

3,000 

3.7 

0.50 

0.70 

4,000 

4.3 

0.67 

0.81 

5,000 

4.8 

0.83 

0.91 

6,000 

5.3 

1.00 

1.00 

100 


- 
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6.6  CONCLUSIONS 


Regime  equations  developed  by  Kellerhals ,  Maddock 
and  the  writer  have  been  compared  and  those  developed  by 
Kellerhals  were  found  to  compare  closely  with  those  deve¬ 
loped  in  this  thesis.  The  equations  developed  in  this 


thesis,  for  coarse-bed  rivers  are: 

1/2 


b  =  2.0  Q 

w  b 


d*  0.32 


0.32 


(6.1) 

(6.2) 


and 


S  = 


0.051  °50 

Qb  0.25 


0.90 


(6.3) 


These  equations  are  based  on  bankfull  discharge  and  the  50% 
bed-material  size  as  determined  from  a  grid  sample  analyzed 
by  number.  They  are  tentative  as  they  are  based  on  rather 
sparse  data  and  on  rivers  that  may  not  have  been  completely 
stable . 


A  chart  having  non-dimensional  parameters  d^/d^ 

versus  Q/ q  has  been  developed  from  which  a  rating  curve 
can  be  deduced..  it  was  noted  that  the  recurrence  interval 
for  bankfull  flows  for  Canadian  rivers  flowing  from  the 
Rocky  Mountains  ranged  from  3  to  15  years  which  is  somewhat 
higher  than  the  recurrence  interval  for  rivers  in  Eastern 


United  States. 


■ 
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6.7  RECOMMENDATIONS 

The  collection  and  analysis  of  more  data  from  stable 
coarse-bed  irrigation  canals  would  be  desirable  in  order  that 
regime-type  equations  can  be  derived  with  more  confidence. 

The  hydraulic  geometry  of  coarse-bed  channels  at 
one  location,  as  well  as  along  a  channel,  should  be  studied 
in  more  detail.  It  is  possible  to  develop  a  rating  curve 
for  a  channel,  provided  that  the  bankfull  discharge  and  depth 
are  known,  from  the  non-dimensional  rating  curve.  However, 
if  the  cross-section  characteristics  of  coarse-bed  channels 
were  to  be  known  in  more  detail  dimensionless  relationships 
between  the  flow  characteristics  and  channel  characteristics 
may  be  developed. 


. 
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CHAPTER  VII 


SCOUR  AT  BENDS  OF 
COARSE-BED  RIVERS 


7.1  INTRODUCTION 

In  locating  highway  and  pipeline  crossings  it  is 
important  to  have  a  reasonable  estimate  of  the  maximum  depth 
of  scour  which  usually  occurs  at  bends  and  constrict¬ 
ions  . 

This  chapter  will  present  depths  of  scour  at  bends 
of  coarse-bed  rivers  in  Alberta  and  relate  them  to  the 
hydraulic  characteristics  of  the  river.  A  discussion  of 
scour  and  fill  at  river  bends,  which  appears  to  be  absent 
in  coarse-bed  rivers,  is  also  presented. 

7.2  SCOUR  AT  BENDS  OF  ALLUVIAL  CHANNELS 

Alluvial  channels  have  been  defined  as: 

open  channels  formed  in  granular,  noncohesive 
material  that  has  been  transported  to  its 
present  site  by  flow  in  the  channel 

(Task  Force  on  Bed  Forms  in  Alluvial  Channels,  1966).  With 

this  definition  it  would  be  in  order  to  describe  a  coarse- 

bed  river  having  sand,  gravel,  cobbles,  and  boulders  on  the 

bottom  as  an  alluvial  river.  The  phenomenon  of  scour  and 

fill  at  bends  of  sand  rivers  is  now  generally  familiar  to 

river  engineers  and  was  well  documented  by  Lane  and  Borland 

(1954)  who  stated  that: 
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•at  high-water  stages,  therefore,  the 
bends  usually  scour  out  and  deposition 
occurs  at  the  crossings;  whereas  during 
low-water  stages,  the  crossings  scour 
out  and  the  bends  fill  . 

This  fact  has  been  supported  by  echo  soundings  during  the 
passage  of  a  flood  on  the  Beaver  River  in  Alberta  (Neill, 
1965)  (see  FIGURE  95)  and  by  soil  borings  after  a  flood  on 
the  Umfolozi  River  in  South  Africa  (Zwamborn,  1966)  (see 
FIGURE  96) .  The  scour  and  fill  process  on  the  Umfolozi 
River  was  also  accurately  reproduced  in  a  hydraulic  model 
with  coal  as  the  bed-material.  It  should  be  noted  that  in 
both  of  these  reported  cases  the  bed-material  was  fine  or 
medium  sand. 

During  the  summer  months  of  1965  and  1966  the 
writer  was  engaged  on  river  regime  surveys  of  coarse-bed 
rivers  in  Alberta.  One  of  the  rivers  investigated  in  some 
detail  was  the  North  Saskatchewan  at  Drayton  Valley  (Galay, 
1967a)  (see  FIGURE  1  for  location)  which  experienced  fairly 
high  flows  in  1965  and  1966  (see  FIGURE  97) .  Successive 
cross-section  soundings  at  various  bends  before,  during  and 
after  the  peak  flow  revealed  several  interesting  facts  (see 
FIGURES  98  and  99) .  Most  of  the  bends  did  not  experience 
any  lowering  of  their  bed  levels  during  peak  flow.  There 
was  some  scour  at  cross-sections  100  and  102  but  this  appears 
to  be  due  to  a  shift  of  the  channel  in  conjunction  with  a 


\ 
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reduction  of  the  bend  radius.  FIGURE  100  illustrates  the 

( 

change  in  the  channel  pattern  which  took  place  during  the 
passage  of  the  1965  peak  flow.  Soundings  at  the  bridge 
(FIGURE  101)  during  low  flow  in  1967  indicate  that  the  scour 
hole  has  not  filled  in. 

Some  of  the  other  cross-sections  also  show  that  the 
scour  holes  did  not  fill  in  as  the  flow  receded. 

Subsequent  to  this  river  survey  at  Drayton  Valley 
a  continuous  longitudinal  sounding  along  the  rivers '  thalweg 
from  Drayton  Valley  to  Edmonton  was  carried  out.  This  sound¬ 
ing  was  conducted  when  the  river  was  relatively  low  and 
indicated  deep  scour  holes  at  river  bends.  Again  it 
appeared  that  there  was  little  filling  of  the  scour  holes 
during  recession  of  peak  flows  (FIGURE  102) .. 

Longitudinal  soundings  along  the  midstream  of  a 
sand-bed  river  (Neill,  1965)  indicate  that  the  depth  at 
bends  was  approximately  equal  to  twice  the  mean  depth  of 
flow,  even  for  low  flow  conditions  (FIGURE  103).  However, 
for  the  North  Saskatchewan  River  the  depth  at  bends  was  some 
three  to  five  times  the  mean  flow  depth  as  shown  by  FIGURE 
102. 

In  the  summers  of  1967  and  1968  a  river  depth  pro¬ 
filing  program  was  carried  out  on  several  Alberta  rivers. 

The  purpose  of  this  program  was  (Hollingshead  and  Schultz, 
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1968) : 

to  obtain  continuous  longitudinal  depth 
profiles  along  major  rivers,  and  conduct 
detailed  depth  surveys  at  bridges,  sharp 
bends  and  other  locations  where  scour 
holes  occur  . 

The  rivers  surveyed  in  this  program  were  the  North  Saskatch¬ 
ewan  downstream  of  Edmonton,  the  Oldman  and  the  Athabasca 
(see  FIGURE  1  for  location) . 

The  survey  methods  used  were: 

1 .  Longitudinal  Soundings 

Depth  profiles  were  obtained  by  sounding 
along  the  river's  thalweg  with  fix  points  being 
located  from  topographic  maps  or  airphoto  mosaics 
and  the  transverse  position  by  a  range  finder. 

2 .  Cross-Sections 

Soundings  of  the  river-bed  and  ground  surveys 
of  the  banks  were  carried  out  at  right  angles  to 
the  thalweg.  At  each  cross-section  information 
as  to  the  bank  and  bed  material,  vegetation, 
erosion,  etc.  was  recorded. 

3 .  Scour  Holes 

At  bends  showing  deep  scour  holes  additional 
cross-sections  and  longitudinal  soundings  were 
taken  in  order  to  determine  the  extent  of  the 


scour . 
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All  these  coarse-bed  rivers  exhibited  relatively 
deep  scour  holes  at  their  bends  as  shown  by  the  soundings 
taken  on  the  Oldman  River  (FIGURE  104) .  An  important 
question  now  arises  -  why  don't  the  scour  holes  at  bends  of 
coarse-bed  rivers  fill  in  during  the  recession  of  peak  flows? 

The  answer  to  this  may  lie  in  the  phenomenon  known 
as  "bed  paving  or  armouring".  This  phenomenon  has  received 
attention  in  rivers  directly  below  dams.  The  bed-load  in 
this  case  is  trapped  in  the  reservoir  resulting  in  a 
deficiency  in  bed-load  downstream  of  the  dam.  The  river 
then  erodes  its  bed  and/or  banks  in  order  to  pick  up  more 
sediment.  This  erosion  of  the  bed  continues  until  the  bed 
is  covered  with  stones  that  are  too  large  to  move.  This 
phenomenon  will  occur  in  rivers  having  sand  beds  with  small 
quantities  of  gravel  as  noted  by  Livesey  (1963)  reporting 
on  the  Missouri  River.  The  rivers  investigated  in  Alberta, 
however,  have  bed-material  which  is  predominantly  larger 
than  1/2  inch  in  size  resulting  in  an  armoured  bed  occurring 
more  frequently.  All  coarse-bed  rivers  display  a  definite 
shingled  or  armoured  effect  during  low-water  stages  (PHOTO¬ 
GRAPH  11) .  This  armouring  must  have  taken  place  at  some 
stage  of  flow  higher  than  the  reference  low  stage  and  may 
account  for  the  fact  that  the  scour  holes  do  not  fill  in 
with  receding  flows.  Let  us  now  examine  what  may  happen  to 
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the  bed  of  a  river  during  the  passage  of  a  peak  flow  (see 
FIGURE  105) . 

As  a  flood  wave  passes  a  definite  river  cross- 
section  the  bed  will  remain  immobile  until  the  mean  velocity 
is  large  enough  to  cause  a  reasonable  number  of  the  smaller 
particles  to  start  moving  (case  b,  FIGURE  105).  As  the  dis¬ 
charge  and  velocity  increase  the  larger  particles  will  be 
moved  in  a  rolling  or  saltating  manner  (Durand,  1951  and 
Thompson,  1965)  with  the  finer  particles  being  carried  off 
in  a  rush  (case  c,  FIGURE  105) .  As  the  flow  recedes  the 
larger  particles  will  segregate  from  the  finer  ones  and  drop 
to  the  bed  to  form  a  paved  or  semi-paved  bed.  Harrison, 

(1950)  reported  that: 

a  complete  layer  of  non-moving  particles 
is  not  necessary  in  preventing  scour  . . . 
only  50  percent  completeness  is  probably 
required  . 

The  finer  bed-material  in  motion  will  continue  to  move  with 
the  flow,  but  there  will  not  be  any  further  supply  of  fine 
particles  from  the  river-bed  due  to  the  paving;  the  flow 
would  be  capable  of  carrying  more  material  along  the  bed 
than  would  be  available  from  the  bed.  Therefore,  the  part¬ 
icles  still  in  motion  would  be  comparatively  small  and  easy 
to  transport  resulting  in  little  or  no  deposition  along  the 
river's  thalweg.  Strong  helicoidal  currents  would  be  present 
at  the  bends  of  rivers  which  would  sweep  out  any  fine  gravel 
or  sand  that  may  be  transported  into  the  scour  holes.  The 
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strength  of  these  helicoidal  currents  can  be  attested  to  by 
the  fact  that  several  120  lb.  concrete  "blocks''  were  deposited 
on  the  top  of  a  high  point  bar  at  the  inside  of  a  bend  on  the 
North  Saskatchewan  River  (PHOTOGRAPH  9) . 

The  hypothesis  that  scour  holes  at  bends  do  not  fill 
in  during  low  stages  makes  it  possible  for  river-bend  surveys 
to  be  carried  out  at  low  stages  on  coarse-bed  rivers.  It 
.should  be  possible  to  correlate  the  scour  data  to  pertinent 
hydraulic  variables. 

7 . 3  TYPES  OF  RIVER  BENDS 

Prior  to  setting  out  the  variables  that  would  in¬ 
fluence  the  depth  of  scour  it  may  be  in  order  to  consider 
the  types  of  river  bends  that  exist  in  nature. 

From  field  and  airphoto  investigations  it  is  apparent 
that  there  are  at  least  two  significantly  different  types  of 
bends : 

1.  Free  bends 

2 .  Forced  bends . 

A  free  bend  is  formed  by  a  river  flowing  on  a  flood  plain 
which  is  free  to  erode  its  banks  and  to  migrate  laterally  or 
longitudinally  (FIGURE  106) .  A  forced  bend,  however,  usually 
impinges  on  a  valley  wall  resulting  in  an  abrupt  change  in 
flow  direction.  The  radius  of  curvature  is  usually  larger 
for  free  bends,  but  the  depth  of  scour  is  usually  smaller 
because  the  scouring  action  is  directed  at  the  bank  as  well 
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as  the  bed  of  the  stream. 

Rzhanitsyn,  (1960)  has  also  classified  river  bends 
into  several  distinct  categories  namely  free  bends,  limited 
bends  and  forced  bends;  this  was',  however,  for  sand-bed 
channels.  This  discussion  will  consider  only  the  two  categor¬ 
ies  originally  mentioned. 

7.4  DIMENSIONAL  ANALYSIS 

Dimensional  analysis  of  the  depth  of  scour  in  re¬ 
lation  to  the  channel  properties  and  the  pertinent  hydraulic 
variables  will  now  be  carried  out.  This  analysis  will  per¬ 
tain  to  a  forced  bend' flowing  at  bankfull  stage  (see  FIGURE 
107)  . 

It  would  appear  that  the  depth  of  scour  d  ,  which 

s 

is  measured  from  the  water  surface,  is  dependent  upon  the 
following  variables : 

Fluid  properties 

p  =  density 
y  =  dynamic  viscosity 
Flow  properties 

Qb  =  bankfull  discharge 
S  =  water  surface  slope 

d*  =  mean  flow  depth  upstream  of  bend  at 
bankfull  stage 


. 
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Channel  properties 

b  =  water-surface  width  at  bend 
w 

r  =  radius  of  curvature 

6  =  internal  angle (in  radians) 

Bed-material  properties 

D  =  typical  particle  size 

p  =  density 
s 

=  gradation  of  bed-material 
Bank  material  properties 

M  =  type  of  bank  material 


and 

g  =  acceleration  of  gravity 
In  equation  form: 

d^  —  f^  ( P  i  y  /  /  S  ,  d^, r  b^ ,  r ,  0  ,  D ,  p  3  ,  CT]-, ,  M ,  g )  (7.1) 

The  bankfull  discharge  can  be  replaced  by  Vm 
through  the  continuity  equation.  Using  the  Buckingham  tt 
theorem  with  p,  Vm  and  b^  as  repeating  variables  results  in 


_s  _  ,  pVmbw  Vm 

,  1  y  7  /gbw' 


S,  d*  r  Q  D  Pc  (7.2) 


bw'  bw'  6/  bT  ' 


•/  a,  M) 


w 


W  P 

Some  of  the  terms  can  be  modified  by  combination 


with  the  d./,  term: 

Dw 


s  _  ^  ,pVmd*  Vm  ^  d*  r  Q  D  ps  _ 

K - f9  (-^-T : /  ~rT->  s'  '  yT“/  0z  t  /  0,  M)  (7.3) 

t>w  ^  P  /gd*  b  d*  p 


The  terms  in  this  equation  can  be  reduced  by  con¬ 


sidering  the  following  assumptions: 


' 


* 
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(1)  Fully  rough  turbulent  flow  conditions  will 
prevail  making  the  Reynolds  number  of  secondary 
importance . 

(2)  The  bed-material  is  similar  in  size,  gradat¬ 
ion  and  density  from  one  river  to  the  next  making 
it  possible  to  discard  f3s/p  and  a. 

(3)  The  bank  material  for  fixed  bends  is  assumed 
to  be  of  resistant  bedrock  varying  little  from 
one  bend  to  another;  the  term  M  is  deleted. 


These  deletions  result  in: 


ds  =  f  /  Vm  q  r_  D_. 

bw  *3  ygd '  '  b  /  bT,'  6'  d.  } 


(7.4) 


'w  ~w  * 

•  At  this  stage  it  is  difficult  to  reduce  the  equat¬ 
ion  to  fewer  dimensionless  terms. 

The  significant  variables  affecting  scour  would  be 
the  radius  of  curvature  r  and  the  internal  angle  6;  the 
depth  of  scour  would  definitely  be  larger  for  a  river  with 
a  severe  bend  as  compared  to  a  gradually  winding  river  (see 


FIGURE  108).  The  importance  of  the  Froude  number 


Vm 

7gd~ 


the 


slope  S,  and  the  ratios  */bw  and  D/^  are  difficult  to 
assess;  insufficient  data  regarding  bankfull  flow  in  coarse- 
bed  rivers  makes  it  difficult  to  graphically  or  statistically 
check  their  influence.  In  the  following  analysis  it  will  be 
assumed  that  the  effect  of  radius  of  curvature  r  and  the 
internal  angle  0  and  the  width  bw  are  dominant  resulting  in: 


. 
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w 


(7.5) 


7.5  ANALYSIS  OF  SCOUR  DATA  FOR 
COARSE-BED  RIVERS  IN  ALBERTA 

« 

The  previously  mentioned  dimensionless  variables 
are  tabulated  in  TABLE  13  and  have  been  obtained  from  field 
surveys  of  Alberta  rivers. 

A  log-log  plot  of  _s  vs  r_  with  0  as  the  third 

.  bw ’  bw 

variable  was  found  to  have  a  large  degree  of  scatter  as 

shown  in  FIGURE  109.  It  is  somewhat  difficult  to  establish 

lines  having  distinct  0  values,  however,  lines  having  0 

equal  to  1.5  and  2.0  are  tentatively  sketched  on. 

d  r 

A  subsequent  plot  of  ^  vs  0  for  in  the  range 

w  d  w 

2.5  to  3.5  (FIGURE  110)  indicates  that  — S-  is  related  linearly 

bw 

to  the  internal  angle  of  the  river  bend.  Since  0  effects 

the  depth  of  scour  in  a  linear  manner  it  was  arbitrarily 

2T  d  2T 

combined  with  the  term  /b  and  a  plot  of  s/b  versus  , — - 

w  w  bw6 

for  fixed  bends  is  shown  in  FIGURE  111.  A  "design  line"  was 

fitted  through  the  top  of  the  plotted  points  as  engineers 

would  be  primarily  interested  in  the  maximum  possible  scour 

at  a  river  bend.  The  relationship  derived  from  this  plot  is: 

ds/b  =  0.14  (7.6) 

w  bw0 

The  data  for  free  bends  was  plotted  in  the  same 
manner  (FIGURE  112)  yielding  another  design  line.  A  compari¬ 
son  of  the  two  "design  lines"  shows  the  depth  of  scour  for 
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TABLE  13 
SCOUR  AT  RIVER  BENDS 

HYDRAULIC  PROPERTIES  OF  BENDS  AT  BANKFULL  STAGE 
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the  same  value  of  g-  will  be  larger  for  a  forced  bend  (see 

w 

FIGURE  113) .  This  confirms  the  previous  statement  that  a 
free  bend  will  have  caving  banks  resulting  in  shallower 
scour  holes. 

The  fact  that  the  "design  plots"  (FIGURES  111  and 
112  possess  a  fair  degree  of  scatter  may  be  attributed  to 
the  neglect  of  several  variables  as  discussed  previously  as 
well  as  the  following  factors: 

(1)  The  assessment  of  the  bankfull  stage  for 
the  Oldman  and  the  Athabasca  Rivers  was  ten¬ 
tative.  Plotting  of  all  the  available  cross- 
sections  on  the  slope  profiles  would  result 
in  more  accurately  established  stages. 

(2)  The  radius  of  curvature  r  and  the  internal 
angle  0  were  obtained  from  aerial  photos  for 
some  rivers  and  from  topographic  maps  for  parts 
of  the  Oldman  River. 

(3)  The  survey  dates  of  the  rivers  investigated 
did  not  coincide  with  the  dates  on  the  air- 
photos.  In  some  instances  it  was  evident  that 
the  surveyed  channel  had  shifted  somewhat  since 
the  photographs  were  taken. 

(4)  The  arc  of  curvature  was  taken  as  being 
parallel  and  midway  between  the  river  banks 
which  in  many  instances  would  not  coincide 
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with  the  thalweg  of  the  river.  This  is 
clearly  indicated  in  FIGURE  114  which  shows 
the  thalweg  as  obtained  from  detailed  bed- 
contours  as  well  the  centre-line  of  the 
channel . 

(5)  Some  bend  cross-sections  were  of  rather 
complex  shapes  requiring  a  rather  arbitrary 
division  of  the  section  into  an  "active" 

and  a  "passive"  portion  (See  FIGURE  98) .  The 
flow  over  the  shallow  portion  at  the  inside 
of  the  bend  would  be  retarded  much  more  than 
the  flow  in  the  "active"  flow  portion  and 
would  probably  have  little  influence  on  the 
scouring  process  at  the  bend. 

(6)  The  scour  depth  was  measured  below  a  bank- 
full  stage,  however,  some  rivers  may  experience 
bankfull  flow  conditions  so  infrequently  that  the 
scour  holes  at  bends  have  not  been  fully  developed. 

(7)  The  geology  from  river  to  river,  as  well  as 
along  the  course  of  an  individual  river,  may  vary^ 
some  bends  may  be  located  at  or  near  subsurface 
rock  outcrops  which  would  limit  the  full  develop¬ 
ment  of  the  scour  hole. 

(8)  The  rivers  may  have  distinctly  different 
channel  patterns,  such  as  braided  and  meandering, 
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which  would  influence'  the  development  of 

secondary  helicoidal  currents. 

(9)  The  bank-material  may  vary  widely; 

some  bends  that  have  been  classified  as  free 

may  have  fairly  resistant  banks. 

7.6  PROCEDURE  FOR  ASSESSING  SCOUR  AT  RIVER  BENDS 

d  r 

The  plots  of  s/u  versus  - — -  for  forced  and  free 

w  b  0 

w 

•river  bends  can  now  be  used  to  assess  the  maximum  scour  depth. 
The  use  of  these  plots  requires  the  following  data: 

(1)  Topographic  maps  and/or  airphotos  in 
order  to  obtain  the  radius  of  curvature 
r  and  the  bend  internal  angle  0 . 

(2)  Cross-section  surveys  in  order  to  assess 
the  bankfull  width  and  to  obtain  a  bankfull 
reference  stage  below  which  the  depth  of 
scour  ds  is  to  be  measured. 

The  procedure  for  establishing  the  bankfull  stage 
has  previously  been  outlined  by  Neill  and  Galay,  (1967)  and 
proceeds  as  follows : 

(a)  Locate  and  survey  a  cross-section  at  the 
bend  in  question  as  well  as  two  cross-sections 
upstream  and  two  downstream  of  this  bend; 

(b)  Plot  the  cross-sections  and  draw  a  tentative 
bankfull  stage  at  each  cross-section,  without 
reference  to  the  other  sections; 
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(c)  Plot  a  tentative  bankfull  stage  at  each 
cross-section  on  the  slope  profile; 

(d)  Draw  an  average  bankfull  stage  line  on 

the  slope  profile  and  use  this  as  the  reference 

* 

plane  below  which  the  depth  of  scour  is  to 
be  measured. 

Therefore,  with  - — ^  computed  one  obtains  a  value 
d  bwD 

for  — S-  and  subsequently  ds .  To  this  depth  of  scour  it  may 
dw 

be  advisable  to  add  several  feet  to  account  for  the  develop¬ 
ment  of  dunes  and  the  fact  that  the  bed  of  a  "gravel"  river 
becomes  very  active  at  high  flows  (Henderson,  1966) . 

In  the  training  of  rivers  many  free  bends  will  be- 

m 

come  forced  bends  with  the  construction  of  revetments  -  the 

relationship  for  forced  bends  should  therefore  be  used  in 

these  instances  to  estimate  the  depth  of  scour  that  one 

would  expect  in  the  future . 

7.7  COMPARISON  OF  DERIVED  SCOUR  DEPTH 
RELATIONSHIPS  TO  REGIME  EQUATIONS 

The  depth  of  scour  for  sand-bed  rivers  can  be 

estimated  by  using  the  regime  relationships  of  Lacey  (1929) 

and  Blench  (1969)  and  these  same  relationships  will  now  be 

assumed  applicable  to  coarse-bed  rivers.  In  order  to  compare 

relationships,  computations  of  scour  depth  were  carried  out 

for  Prairie  Creek  near  Rocky  Mountain  House  (Alberta  Water 

Resources,  1967) .  The  computed  depths  of  scour,  for  a  free 
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bend,  are  as  follows: 

This  Thesis  Actual 

LACEY  BLENCH  (Figure  112)  Scour 

Depth 

of 

Scour 

(ft.)  6  10  11  9 

The  detailed  computation  steps  are  presented  in  APPENDIX  9. 

The  depth  of  scour  as  obtained  from  Blench' s  relationships 
and  that  derived  from  FIGURE  112  are  very  close  to  the  actual 
scour  while  the  results  from  Lacey's  relationships  under¬ 
estimate  the  scour  depth.  However,  this  comparison  needs  to 
be  carried  out  for  a  large  number  of  bends  before  more  definite 
conclusions  can  be  made. 


7.8  CONCLUSIONS 

Longitudinal  thalweg  soundings  on  a  limited  number 
of  coarse-bed  rivers  in  Alberta  have  revealed  that  scour 
holes  at  bends  do  not  scour  and  fill  during  the  passing  of 
high  flows.  This  phenomenon  is  unique  to  coarse-bed  rivers 
as  sand-bed  rivers  have  their  scour  holes  filled  in  during 
the  recession  of  high  flows.  However,  more  detailed  field 
investigations  are  necessary  in  order  to  conclusively  state 
that  scour  holes  do  not  scour  and  fill. 

Dimensional  analysis  was  used  in  an  attempt  to  re¬ 
late  the  depth  of  scour  to  pertinent  channel  and  flow  variables, 
however,  a  relationship  between  the  depth  of  scour  and  sever¬ 
al  dimensionless  variables  was  not  possible  as  too  many 


variables  were  involved.  Treating  only  the  variables  of 
major  significance,  namely  width  of  river  bw,  radius  of 
curvature  r,  and  internal  angle  0  resulted  in: 
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f  -  f 4  e) 

^w  Dw 


(7.5) 


Using  data  from  Alberta  Rivers,  and  classifying  bends  into 
forced  and  free  bends  resulted  in  the  following  tentative 
design  equations: 

forced  bends : 


'w 

and  free  bends : 
ds 


f  -  °-14 

b™  bw 


-0.48 


(7.6) 


-0.38 


=  0.11  ( 


'w 


b  e 

w 


) 


(7.7) 


The  derived  design  formula  for  forced  bends  was 
applied  to  a  typical  coarse-bed  river.  Prairie  Creek  near 
Rocky  Mountain  House,  and  the  derived  scour  depth  compared 
closely  to  the  actual  depth  as  well  as  to  the  depth  obtained 
by  using  Blench 's  equation  for  zero  flood  depth. 


7.9  RECOMMENDATIONS 

The  basic  phenomenon  of  scour  at  bends  of  coarse- 
bed  rivers  needs  far  more  study.  Field  surveys  using  echo 
sounders,  undertaken  during  the  passage  of  flood  flows,  are 
necessary  to  conclusively  state  that  scour  holes  do  not 
experience  scour  and  fill.  The  recording  of  the  modification 
of  bed  forms  as  they  pass  through  scour  holes  would  be  of 
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interest.  Surveys  should  also  be  conducted  at  constrictions; 
the  secondary  currents  at  constrictions  may  differ  distinctly 
from  those  at  bends  resulting  in  a  different  scour  and  fill 
process . 

Laboratory  research  on  the  transport  and  segregation 
of  bed-material  mixtures  under  unsteady  flow  conditions  could 
clarify  the  bed-paving  or  sorting  process  that  is  so  distinct 
in  natural  rivers.  Also,  in  hydraulic  models  of  coarse-bed 
rivers  it  has  been  common  to  use  a  coarse  light-weight  mat¬ 
erial  having  only  one  size  as  the  bed  material  in  the  model. 
This  uniform  size  of  material  may  severely  restrict  the  model 
from  reproducing  the  scour  process  as  it  occurs  in  nature  if 
this  process  is  highly  dependent  on  the  sorting  of  the  various 
sizes.  This  sorting  phenomenon  will  not  occur  in  the  model  - 
one  may  find  that  the  scour  holes  will  fill  in  and  that  the 
losses  would  not  be  reproduced  correctly.  Model  studies  of 
scour  at  river  bends,  using  various  types  of  materials,  would 
be  most  informative. 
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CHAPTER  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS 

8.1  introduction 

.  At  the  outset  of  this  study  the  original  purpose  was 
to  extend  the  knowledge  of  coarse  bed-material  transport  by 
investigating  the  movement  of  bed-material  in  the  North 
Saskatchewan  River.  However,  it  soon  became  apparent  that 
the  behaviour  of  rivers  having  coarse  beds  was  extremely 
complex.  The  study  then  expanded  into  an  evaluation  of  the 
more  critical  problems  confronting  engineers  dealing  with 
these  rivers.  The  problems  considered  dealt  with  bed-material 
sampling  techniques,  threshold  of  motion,  bed-material  trans¬ 
port,  resistance  to  flow,  stable  channel  design  and  scour  at 
river  bends.  After  analyzing  the  available  information  in 
each  of  these  topics,  and  they  are  by  no  means  completely 
distinct  from  each  other,  conclusions, and  recommendations 
regarding  further  research  were  formulated. 

8.2  CONCLUSIONS 

1.  The  techniques  used  in  sampling  the  river-bed  are 
directly  related  to  the  purpose  of  sampling  which  ususally 
involves  the  threshold  of  motion,  the  volume  of  material  in 
transport  or  the  resistance  that  a  coarse  bed  offers  to  the 
flow.  Sampling  the  bed  in  the  following  two  ways  would 
generally  serve  the  mentioned  purposes: 
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(a)  Scoop  sample  at  several  depths  below 
the  surface,  analyzed  by  weight. 

(b)  Grid  sample  by  a  paced  or  taped  grid, 
analyzed  by  number. 

2.  In  order  to  assess  the  threshold  of  motion  for 
coarse-bed  rivers  having  a  range  of  stone  sizes  field  invest¬ 
igations  were  conducted  on  the  North  Saskatchewan  River  and 
Wilson  Creek.  The  data  obtained  were  combined  with  data  from 
other  investigators  and  analyzed  in  several  different  ways , 
namely : 

(a)  Mean  critical  velocity  related  to  stone  size. 

(b)  Critical  shear  stress  related  to  stone  size. 

(c)  Shields  mobility  number  related  to  particle 
Reynolds  number. 

(d)  Modified  mobility  number  related  to  relative  depth 
and 

(e)  Blench' s  zero  bed  factor  related  to  relative  depth. 
After  examining  the  various  relationships  the  relationship 
between  mean  critical  velocity  and  stone  size: 

Vmc  =  8.0  D1/3 

was  recommended  for  assessing  when  large  stones  in  a  mixture 
will  start  to  move.  The  velocity  in  this  relationship  was  the 
mean  velocity  in  a  cross-section,  not  the  mean  in  a  vertical 
distribution.  For  this  reason  it  was  difficult  to  directly 
compare  this  relationship  to  similar  relationships  developed 
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by  other  investigators. 

Of  the  various  dimensionless  relationships  invest¬ 
igated  that  of  Blench: 

F, 


bo 


=  29  (D/d*) 


1/2 


appeared  to  be  reasonable. 

3.  Existing  sediment  transport  data  utilizing  coarse 
materials  were  analyzed  through  the  use  of  dimensional 
analysis.  This  analysis  indicated  that  there  are  insufficient 
data  to  develop  relationships  that  would  be  applicable  to 
natural  rivers.  A  subsequent  replotting  of  the  data  in  a 
simplified  form  having  the  unit  bed-load  discharge  versus 
mean  velocity  with  relative  depth  as  the  third  variable 
proved  promising.  The  plotting  of  data  in  this  form  would 

be  more  acceptable  to  practicing  engineers . 

A  comparison  of  the  computed  bed-load  transport, 
using  several  common  formulas,  to  the  actual  measured  trans¬ 
port  in  the  Elbow  River  indicated  that  no  one  formula  can 
be  recommended  over  all  the  others.  Computing  the  bed-load 
transport,  in  the  North  Saskatchewan  River,  from  shifting 
gravel  bars  yielded  transport  rates  that  compared  favourably 
with  computations  based  on  various  formulas. 

4.  The  measurement  of  the  direct  protrusion  height 
of  stones  on  a  river  bed,  using  a  roughness  meter,  resulted 
in  the  following  resistance  formula  being  developed  for 
rigid  coarse-bed  channels: 
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=  0.67  D^q  (by  number).  The 
instead  of  a  computed 
formula  would  make  it 
from  flume  studies  with 

convert  the  preceding 
roughness  coefficient  n  and 


5 .  For  mobile  coarse-bed  channels  no  equations  or 
charts  have  been  developed  to  assess  flow  resistance;  as  in 
sediment  transport,  the  data  are  insufficient.  Several 
interesting  observations,  from  flume  experiments  and  field 
studies,  are,  however,  worth  noting.  Flume  data,  for  uniform 
materials,  tentatively  indicate  that  the  size  of  material  in 
motion  has  little  influence  on  resistance.  Limited  field 
data  suggests  that  the  bed-material,  once  in  motion,  forms 
into  dunes  and  then  into  a  plane  bed  as  the  flow  intensity 
increases. 

6.  For  the  design  of  stable  channels  in  coarse  bed- 
material  several  equations  have  been  compared  and  those 


^7=  3-°  <ih> 


0.45 


In  this  formula  the  protrusion  k 
use  of  direct  protrusion  height, 
equivalent  grain  size,  in  a  flow 
comparable  to  formulas  developed 
artificial  roughness  elements. 

Engineers  may  prefer  to 
equation  into  terms  of  Manning's 
this  would  result  in: 


n 
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0.088 

"(d*/k) 


0.45 


where  k  =  0.67  D^q  (by  number) 
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developed  in  this  thesis  are  recommended 

0.50 


bw  =  2‘° 


d.  =  0.32  Q. 


0.32 


and 


S  = 


0 


nc i  D  °-90 

.05!  u 5 ^ 


Q, 


0.25 


These  equations  are  based  on  bankfull  discharge  and 
the  D5q  bed-material  size  as  determined  from  a  grid  sample 
analyzed  by  number.  The  data  used  in  developing  these 
equations  are  from  rivers  in  North  America  and  the  derived 
equations  may  not  apply  to  coarse-bed  rivers  elsewhere. 


7.  A  chart  having  non-dimensional  parameters 

versus  has  been  developed  from  which  a  rating  curve  for 

coarse-bed  rivers  can  be  deduced.  Bankfull  discharges  can 
also  be  estimated  if  cross-section  and  flow  data  at  some  low 
flow  are  available. 

8.  During  river  regime  surveys,  echo  soundings  at 
river  bends  having  coarse  beds  revealed  that  scour  holes, 
once  established,  are  relatively  stable.  They  do  not  under¬ 
go  a  scour  and  fill  process  as  is  the  case  for  sand-bed  rivers 
during  the  passage  of  a  flood. 

From  an  analysis  of  the  dominant  variables  influenc¬ 
ing  scour  at  river  bends  the  following  relationships  for  the 
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maximum  scour  depth  have  been  developed 
Fixed  bends : 

-0.48 

^  =  0.14  (=r_) 

bw  Dw  e 


Free  bends : 


ds 

b 


-0.38 


=  0.11  <g^) 

W  W 


These  equations  are  tentative  as  the  available  data  is  re¬ 
latively  inconsistent  and  sparse. 


8.3  RECOMMENDATIONS 

After  treating  a  number  of  problems  associated  with 
coarse-bed  rivers  it  becomes  apparent  that  adequate  field 
as  well  as  experimental  data  are  lacking.  The  following  re¬ 
commendations  aim  at  overcoming  this  deficiency. 

8.3(a)  Field  Investigations 

Field  investigations  should  be  organized  to  obtain 
a  wide  range  of  useful  data  from  one  river  or  canal. 

Several  of  the  problems  dealt  with,  namely  the  threshold  of 
motion,  the  transport  of  sediment,  the  resistance  to  flow 
and  the  design  of  stable  channels  can  be  investigated 
simultaneously  on  one  or  more  rivers.  These  investigations 
should  be  carried  out  on  the  following  types  of  channels : 

1.  Straight  irrigation  canal  or  river  having 
a  relatively  low . flow  (bankfull  discharge  up 
to  5000  cfs) . 
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2.  Straight  stretch  of  river  having  a  re¬ 
latively  high  bankfull  flow  (up  to  10,000  cfs) . 

The  field  program  should  be  organized  so  that  the 
following  information  is  obtained: 

1.  Size  distribution  of  surface  and  subsurface 
bed-material,  as  well  as  bank  material. 

2.  Variation  of  bed-material  sizes  with  their 
location  in  the  channel. 

3.  The  hydraulic  conditions  under  which  a 
reasonable  portion  of  the  bed  is  in  motion. 

4.  The  amount  of  bed-material  moving  along 
the  channel  and  the  portion  of  the  channel 
bottom  that  is  mobile. 

5.  The  resistance  to  flow  that  a  rigid  bed 
offers  to  the  flow  and  the  change  in  this 
resistance  as  the  bed  becomes  highly  mobile. 

A  continuous  record  of  the  bed  forms  gen¬ 
erated  would  be  useful  in  correlating  the 
bed  form  dimensions  to  the  corresponding 
resistance . 

6.  The  values  of  width,  depth,  and  slope  at 
bankfull  flow  conditions. 

A  detailed  field  investigation  on  the  behaviour  of 
scour  holes  at  bends  should  also  be  undertaken.  A  sounding 
program  to  obtain  detailed  bed  contours  before,  during  and 
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after  the  passage  of  flood  flows  would  be  necessary. 

8.3(b)  Laboratory  Investigations 

At  the  outset  it  is  stressed  that  laboratory  invest¬ 
igations  should  deal  primarily  with  coarse  bed-material  mix¬ 
tures  instead  of  uniform  material. 

In  studying  the  threshold  of  motion  the  bed-material 
should  initially  be  sorted  by  comparatively  high  flows.  The 
resulting  data  would  therefore  be  applicable  to  natural  rivers 
the  placing  of  material  in  the  bottom  of  a  flume  with  the  use 
of  screeds  may  result  in  bed  particle  orientations  that  are 
not  representative  of  field  conditions.  The  threshold  of 
motion  state  may  have  to  be  assessed  by  hydrophone  or  by 
radioactive  labelling  of  particles. 

As  the  intensity  of  flow  increases  the  bed-material 
becomes  mobile;  the  continuous  measurement  of  transport  rates 
as  well  as  development  and  movement  of  bed  forms  would  yield 
valuable  data  on  sediment  transport  and  resistance  to  flow. 

The  investigations  should  be  conducted  with  mixtures  having 
a  median  size  of  at  least  20  mm  and  a  relative  depth  up  to 
40.  These  investigations  should  also  be  carried  out  with 
various  materials  having  a  range  of  densities. 

Laboratory  investigations  dealing  with  the  develop¬ 
ment  of  scour  holes  at  river  bends  having  coarse  materials 
may  clarify  field  observations  and  indicate  the  more  signifi¬ 
cant  variables  involved  in  the  scour  process.  These  studies 
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would  require  a  fairly  large  model  and  various  types  of  light¬ 
weight  materials  to  simulate  bed-material  mixtures. 
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APPENDICES 


APPENDIX  1 


LIST  OF  SYMBOLS 


SYMBOLS 


A 

a 

a 

b 

b 

b 

w 

C 

Cl 

cw 

c 

D 

Di 


Dm 


Fb0 

Fs 


Flow  area  (sq.  ft.) 

Major  axis  of  bed-material  particle 
Constant 

Channel  width  (ft) 

Intermediate  axis  of  bed-material  particle 
Water  surface  width  (ft) 

Bed  load  charge  (pt.  per  1000,000  by  weight) 
Constant 

Concentration  of  suspended  particles 

Minor  axis  of  bed-material  particle 

Particle  size  (ft)  of  bed-material 

The  particle  size  corresponding  to  the  ith 
percentile  on  the  particle  size  distribution 
curve  for  bed-material  (ft) 

Geometric  mean  particle  size 

Representative  particle  size  (ft) 

Mean  flow  depth  =  bw  (ft) 

Bankfull  depth  (ft) 

Depth  of  scour  (ft) 

Bed  factor  (Blench) 

Zero  bed  factor  (Blench) 

Side  Factor  (Blench 


\ 


' 


.  ■ 


f  Function 

f  Lacey's  silt  factor 

fc  Factor  defining  cross-sectional  shape  of  channel 

fg  Factor  defining  plan  geometry  of  channel 

2 

g  Acceleration  of  gravity  (ft/sec  ) 

k  von  Karman  turbulence  coefficient 

ks  Equivalent  grain  size  (often  D^) 

k  Protrusion  height  of  bed-material 

M  Type  of  bank  material  parameter 

n  Manning's  roughness  coefficient 

P  Dimensionless  bed-load  transport  parameter  (-^s-  -) 

pv*3 

P  Wetted  Perimeter  (ft) 

3 

Q  Discharge  (ft  /sec) 

3 

Qb  Bankfull  discharge  (ft  /sec) 

3 

Qd  Dominant  discharge  (ft  /sec) 

qs  Weight  of  bed-load  transported  per  unit 

time  per  unit  width  (lb/ft. sec) 

R  Hydraulic  radius  (ft) 

R^  Hydraulic  radius  of  the  bed  (ft) 

R^  Hydraulic  radius  with  respect  to  particle  (ft.) 

Rn  Reynolds  number 

r  Radius  of  curvature  of  river  bend  (ft) 

S  Slope  of  energy  grade  line 

Vm  Mean  flow  velocity  based  (ft/sec) 

Vmc  •  Critical  mean  flow  velocity  for  the  threshold 
of  motion  of  bed-material  (ft/sec) 


i  • 


. 


2 . 


. 


V 

* 

Wp 

w 

X 

Y 
Z 


Shear  velocity  (ft/sec) 

Weight  of  particle  (lb) 

Fall  velocity  of  particle 

Particle  Reynolds  number 

2 

Mobility  number 

Relative  depth  (d^ ) 

D 


(ft/sec) 

(pDv*) 

y 


ab 

Y 

y' 

s 

e 

x 

y 

V 

p 

ps 

tb 

ag 

TO 

TC 

4> 

X 


Shape  factor  of  bed-material 
Specific  weight  of  water  (lb/^^3) 

Bouyant  weight  of  bed-material  (lb/^  3) 

•Internal  angle  of  curvature 

Concentration  of  bed-material  or  protrusions 

2 

Dynamic  viscosity  of  water  sediment  mixture  (lb  sec/ft  ) 

2 

Kinematic  viscosity  of  water  sediment  mixture  (ft  /sec) 

2  4 

Mass  density  of  fluid  (lb  sec  /ft  ) 

2  4 

Mass  density  of  bed-material  (lb  sec  /ft  ) 

Gradation  of  bed-material 

Geometric  standard  deviation 

Shear  stress  acting  on  the  bed  (lb/^2) 

Critical  shear  stress  acting  on  the  bed  for 
the  threshold  of  motion  of  bed-material  (lb/^^2) 

Einstein ’ s , bed-load  function 

Areal  pattern  of  bed-material  or  protrusions 
Einstein's  intensity  of  shear  parameter 
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TABLE  A-l 

NORTH  SASKATCHEWAN  RIVER 
NEAR  DRAYTON  VAUEY 
CHANNEL  CROSS  SECTION  DIMENSIONS 


River 

Stretch 

Cross 

Section 

No. 

Height 

Pef . 
to 

Peak-flow 

Stage 

(ft.) 

RE 

Cross 

Section 

Area 

(ft.*) 

FERENCE  STAG 

Water- 

Surface 

Width 

(ft.) 

£ 

• 

Kean 

Depth 

(ft.) 

A/Ws 

P 

Cross 

Section 

Area 

(ft.*) 

EAK-EIOW  STA 

Water- 

Surface 

Width 

(ft.) 

GE  (1565-66) 

Kean 

Depth 

(ft.) 

A/Ws  • 

Maximum 

Depth 

(ft.) 

S2 

9.5 

960 

570 

1.5 

6800 

638 

10.7 

14.0 

S3 

9.5 

2160 

513 

4.2 

7900 

624 

12.7 

15.0 

54 

9.5 

1700 

515 

3.3 

8020 

775 

10.4 

15.7 

55 

9.0 

2200 

500 

4.4 

7240 

590 

12.3 

15.0 

A 

56 

9.0 

2120 

407 

5.2 

7100 

710 

10.0 

16.5 

57 

9.0 

2030 

520 

4.0 

7180 

600 

12.0 

16.2 

58 

8.5 

1200 

325 

3.7 

5840 

820 

7.1 

13.3 

59 

e.5 

ie40 

377 

4.9 

6030 

670 

9.1 

15.0 

60 

8.0 

1500 

363 

3.9 

6760  ' 

845 

8.0 

13.5 

61 

7.5 

920 

280 

3.3 

8120 

1020 

8.0 

13.2 

62 

8.5 

1320 

640 

2.1 

7100 

830 

8.6 

12.5 

63 

8.5 

1340 

450 

3.0  * 

7120 

710 

10.0 

12.7 

64 

9.0 

560 

190 

3.1 

7860 

840 

9.4 

13.6 

65 

9.3 

1260 

440 

2.9 

7420 

740 

10.0 

15.6 

B 

66 

10.0 

1000 

590 

1.7 

8320 

855 

9.7 

13.0 

67 

10.0 

1400 

325 

4.3 

5720 

520 

11.0 

18.0 

68 

9.5 

1100 

5C0 

2.0 

8600 

8S0 

9.8 

14.0 

69 

9.3 

1440 

3S0 

3.8 

8440 

920 

9.2 

14.8 

. 

70 

9.0 

7C0 

370 

1.9 

7520 

930 

9.1 

12.2 

71 

9.0 

1520 

250 

6.1 

4460 

450 

9.3 

19.3 

72 

9.0 

2100 

450 

4.7 

12280 

1390 

8.8 

15.8 

73 

9.0 

2140 

530 

4.0 

12140 

1214 

10.0 

13.9 

74 

8.5 

2420 

620 

3.9 

135C0 

1315 

10.3 

16.2 

75 

8.5 

20C0 

440 

4.6 

6660 

640 

10.4* 

14.8 

76 

8.3 

2760 

450 

6.2 

6600 

510 

12.9 

16.5 

77 

8.0 

2140 

375 

5.7 

5900 

560 

10.5 

16.3 

78 

8.0 

2200 

420 

5.3 

6200 

553 

11.2 

19.1 

79 

8.0 

1840 

500 

3.7 

7660 

725 

11  .6 

14.2 

80 

7.5 

1500 

490 

3.9 

6620 

685 

9.7 

15.0 

81 

7.5 

1340 

430 

3.1 

6760 

720 

9.4 

14.1 

82 

8.0 

440 

225 

1.5 

6100 

1170 

5.2 

10.1 

83 

8.0 

1722 

290 

5.9 

5142 

650 

7.9 

17.0 

84 

8.0 

1480 

590 

2.5 

9100 

845 

10.8 

13.4 

es 

8.0 

1950 

400 

4.9 

7620 

850 

9.0 

14.0 

86 

8.0 

1840 

360 

5.1 

6340 

720 

8.9 

17.9 

87 

6.0 

1440 

503 

2.9 

5700 

545 

10.5 

16.0 

88 

8.0 

2050 

330 

6.2 

6800 

E80 

10.0 

19.0 

89 

8.3 

1200 

365 

3.3 

7460 

810 

9.2 

13.0 

90 

8.5 

2460 

650 

3.8 

8100 

745 

10.9 

15.8 

c 

91 

8.5 

1940 

475 

4.1 

5160 

725 

7.1 

16.2 

92 

8.5 

2200 

685 

3.2 

8520 

760 

11.2 

14.8 

93 

8.5 

1000 

460 

2.2 

114SO 

1440 

8.0 

13.3 

94 

8.5 

I860 

430 

4.3 

54C0 

1310 

7.2 

15.9 

95 

8.0 

1600 

460 

3.5 

7300 

740 

9.9 

13.0 

96 

8.0 

1940 

490 

4.0 

8920 

890 

10.0 

17.3 

97 

8.0 

1220 

360 

3.4 

9380 

1040 

9.0 

12.5 

98 

8.3 

1360 

229 

6.0 

6220 

780 

8.0 

20.3 

99 

8.3 

2650 

415 

6.4 

9140 

960 

9.5 

17.5 

100 

8.0 

2540 

290  ' 

8.8 

7280 

716 

10.2 

30.5 

101 

8.0 

ECO 

470 

1.6 

9580 

1013 

9.5 

21  .5 

102 

8.0 

3700 

510 

7.2 

7720 

570 

13.5 

24.8 

103 

8.0 

2240 

420 

5.3 

10420 

1055 

9.9 

18. 0 

104 

8.0 

2120 

460 

4.6 

12660 

1840 

7.0 

15.3 

105  " 

6.0 

960 

416 

2.3 

esso 

1182 

7.6 

12.0 

105 

7.S 

1580 

540 

2.9 

7220 

770 

9.4 

12.2 

107 

7.6 

760 

240 

3.2 

5200 

732 

6.7 

12.5 

108 

7.6 

14  =  0 

600 

2.5 

6960 

740 

9.4 

12.5 

1C9 

7.5 

14C0 

530 

2.6 

7040 

810 

8.7 

13.5 

110 

7.5 

1940 

400 

4.9 

8140 

1125 

7.2 

18.0 

111 

7.5 

1680 

350 

4.8 

104S0 

1333 

7.8 

14.5 

112 

- 

- 

- 

- 

- 

- 

* 

- 

113 

8.0 

1120 

330 

3.4 

10029 

1003 

10.1 

16.5 

114 

7.6 

2050 

515 

4.0 

8060 

eso 

9.2 

14.0 

115 

7.8 

3060 

410 

7.5 

10060 

1030 

9.8 

19.0 

116 

7.8 

2000 

370 

5.4 

7780 

645 

9.2 

16.2 

117 

7.5 

2050 

600 

3.4 

7020 

730 

9.6 

14.0 

o 

116 

6.0 

1650 

450 

3.8 

6820 

550 

10.6 

17.5 

119 

7.8 

14S0 

550 

2.7 

9430 

1130 

8.4 

11.8 

120 

7.8 

2640 

505 

5.2 
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TABLE  A- 3 


NORTH  SASKATCHEWAN  RIVER 
NEAR  DRAYTON  VALLEY 
ANALYSES  OF  COARSE  BED-MATERIAL 


Cross 

Section 

LOCATION' 

Part  of 
Channel 

Median  Size 
(in) 

ANAL 

Type  of 
Analysis 

-YSES 

Mean  Size 
(in) 

Sorting  (Geon. 
Standard  dev.) 

SAMPLING 

TECHNIQUE 

es  d/s 

Right  point  bar 

2.2 

\  by  nur.be  r 

2.3 

1.50 

Line  (tape)  -  2  ft. 

87  D/S 

Right  point  bar 

2.3 

\  by  number 

2.5 

1.78 

Line  (tape)  -  2  ft. 

If 

•  1 

2.4 

V  by  number 

2.6 

1.67 

Line  (pace)  -  3  ft. 

87  D/S 

In  channel 

2.4 

%  by  number 

2.6 

1.71 

Line  (pace)  -  3  ft. 

90  U/S 

Right  point  b3r 

1.4 

%  by  number 

1.5 

1.71 

Line  (tape)  -  1  ft. 

If 

II 

1.1 

\  by  weight 

- 

- 

Sub-surface  scoop 

90 

Right  point  bar 

1.1 

\  by  nurber 

1.2 

1.81 

Line  (tape)  -  1  ft. 

«! 

It 

0.7 

%  by  weight 

- 

- 

Sub-surface  scoop 

94  U/S 

Right  point  bar 

'  1.2 

\  by  number 

1.3 

1.66 

Grid  (lines) 

If 

II 

0.8 

%  by  weight 

- 

- 

Sub-surface  scoop 

94 

It 

1.6 

\  by  number 

1.8 

1  .87 

Line  (tape)  -  1  ft. 

II 

•  1 

0.9 

\  by  weight 

- 

- 

Sub-surface  scoop 

100 

Left  point  bar 

2.0 

%  by  number 

2.1 

1.50 

Line  (pace)  -  3  ft. 

ft 

II 

1.7 

\  by  number 

1.9 

1.60 

Line  (tape)  -  2  ft. 

100  U/S 

Left  Point  bar 

1.1 

V  by  lumber 

1.2 

1.55 

Line  (tape)  -  1  ft. 

If 

If 

1.4 

%  by  weight 

- 

- 

Sub-surface  scoop 

100 

Left  point  bar 

1.9 

%  by  ntnber 

2.2 

1.8S 

Line  (tape)  -  1  ft. 

II 

•  I 

1.1 

4  by  weight 

- 

- 

Sub-surface  scoop 

100  D/S 

Left  point  bar 

1.4 

%  by  nur.be  r 

1.5 

1 .70 

Line  (tape)  -  1  ft. 

If 

II 

1.0 

\  by  weight 

- 

- 

Sub-surface  scoop 

101 

Right  point  bar 

2.1 

V  by  number 

2.2 

1.60 

Line  (tape)  -  1  ft . 

II 

•  I 

1.3 

%  by  weight 

- 

- 

Sub-surface  scoop 

101  D/S 

Right  point  bar 

1 .9 

\  by  nur.be  r 

- 

- 

Line  (tape)  -  1  ft. 

It 

It 

l.S 

%  by  weight 

- 

- 

Sub-surface  scoop 

104 

Right  bar 

1.5 

V  by  number 

1 .6 

1.85 

« 

Line  (tape)  -  1  ft. 

II 

If 

1.1 

%  by  weight 

- 

- 

Sub-surface  scoop 

104  P/S 

Right  bar 

1.1 

%  by  number 

1.2 

1.55 

Line  (tape)  -  1  ft. 

*• 

II 

0.7 

%  by  weight 

* 

Sub-surface  scoop 

. 

’ 


• 

TABLE  A- 4 


ANALYSES  OF  BANK-MATERIAL  SAMPLES 


LOCATION 

ANALYSES 

REMARKS 

Cross 

Section 

Part 
of  Bank 

Percentages  of 
Gravel  Sand  Silt 

Clay 

59 

Upper  part  of 
right  point  bar 

- 

80 

15 

5 

Recent  deposition 

92 

Right  eroding 
bank 

1 

44 

- 

55  - 

Flood  plain  deposit, 
one  foot  below  ground 

92 

Right  flood  plain 

- 

80 

14 

6 

Four  feet  below  ground 

55 

Left  bank 

- 

77 

— 

23  - 

Four  feet  above  water, 
recent  deposit 

APPENDIX  3 


RIVER  DATA 

THRESHOLD  OF  MOTION  OF  BED-MATERIAL 
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APPENDIX  A 
FLU HE  DATA 


COARSE  BED  MATERIAL 


START 


V 


7 


READ 
INTERPOLATION 
DATA  SETS 
_ / 


WRITE  / 
GENERAL  / 

^heading/ 


SKIP  TO 
TOP  OF  NEW 
PAGE 


-  -  SP¬ 


READ 

DATA 


\  WRITE  / 
CALCULATED 
\  DATA  j 


CALCULATIONS 

V  =  discharge/ (breadth*depth) 
m 

P  =  breadth  +  2*depth 
R  =  breadth*depth/P 
SLOPE  =  SLOPE/lOO . 


V*  =  32 . 172*R*SLOPE 

F  =  8*(V*/Vm)2 

R  =  4*V  *R/Kin  Vis 
n  m 

RNF  =  Rn//F 


T 

ISN 

RNF  WITHIN  N 
INTERPOLATION 
DATA  SET^ 

V 

yes 

J 

CALL  SUBROUTINE 
TO  INTERPOLATE 
FOR  FW 


1/2 


no 


-7 

RNF  =  RNFX 10 
R  =  R  *10"6' 

n  n  _  o 

X  =  X*10 

RBD50  =  R.  /D50 

13  3  2 

X  =  y's  *D50  /(p*Kin  Vis  ) 

PARTI  =  PART/(p*V*3) 

PART  =  charge/1600* (discharge 
/breadth) 

Y  =  p*V*2/(y  '  s*D50 ) 


D50  =  D50/304 . 8 

V*b  =  (32 .172*Rb*SLOPE) 

R  =  R*F  /F 

b  ®  1/2 

C  =  (8/Fb)±/Z 

CALCULATIONS 


1/2 


F,  =  F+2*depth* (F-FW) 
/breadth 


NOTE 


T  t 

* 


means  multiply  except  where  it 
is  used  as  a  subscript 


COMPUTER  PROGRAM-  BED-LOAD  TRANSPORT 

FIGURE  A- IQ 
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TA8LF  A  -  6 

ANALYSIS  OP  CDAPSF  BED- MATER  I AL  Ft  IJMF  DATA 
GILBERT  -  G p  A D E  E ,  1.7.3  MR,  UNIFORM 


T  p  ST 

V  M  c  ^  N| 

F 

C 

RB 

V*P 

Y 

PARTI 

P  P  /  D  6  0 

3-  1 

2o  2  3 

0©  05  8 

10.66 1 

0.13  1 

0.  2C7 

0o  141 

4.412 

1  9.  8 

1-  2 

2.^6 

0,051 

11,486 

0o  3.03 

0o  205 

0o  1  70 

40  5?4 

1  P,  4 

1-  3 

2.  4;*) 

Oo  0 55 

1 1. 032 

0o  102 

0.217 

0.3  55 

4.958 

18,7 

1  _  u 

* 

7.  43 

0,  056 

1  0o  84  7 

Oo  103 

0e  721 

0,  1  61 

4,  9?q 

1  P.  3 

1-  5 

2o  4,2 

0  ©  0  5  8 

10, 6« 3 

Go  102 

0.226 

0.  3  69 

4.515 

ip.  2 

3-  6 

.  ?o  “'5 

Oo  044 

11 O  701 

0o  165 

0o  2  35 

0o  3.  8? 

4,709 

2  c,4 

3  -  7 

2o  70 

Oo  0  52 

10o 63? 

0.167 

0o  263 

0.225 

4.  36  7 

2  6.  8 

3-  8 

2  «  7  5 

Go  05  3 

1 0o 470 

0o  3.70 

0.263 

0.227 

4.499 

70.  3 

1  -  Q 

2e  97 

Oo  045 

11.  7  3  5 

0  ©  1  5  5 

0.2  53 

0o  2  3  1 

5.122 

2  7  0  6 

3  -3  0 

20oo 

OoOP? 

515 

0o  ?P8 

0.334 

0.3  25 

2o  64  5 

5  9,  1 

3-13 

2.  77 

Oo  063 

8.,  612 

0.318 

0.3  34 

0.367 

2,275 

5  6.  6 

1-1  ? 

7.  07 

Oo  03  3 

11. 627 

0.374 

0o  2  59 

0.  2  ?  2 

4,  A  59 

6  6,6 

1-13 

1 , 94 

0.060 

10.939 

0o  088 

Oo  177 

0.103 

3.  912 

1  c.  6 

3-3  4 

2o  Ca 

Oo  040 

1.  2,  7  76 

0o  0P1 

0.  3.  70 

Oo  C  96 

4  0  3  P  4 

14,4 

1-15 

?  o  2  5 

0.064 

10, 677 

0o  076 

0.  210 

0,  3  4  6 

5o  5  36 

13,6 

3-36 

7  o  3  4 

Oo  0 ^  5 

3.1  ©  5  45 

0,07? 

0.205 

Go  3  3P 

6.217 

3  2,  9 

1-17 

lo  9  6 

Oo  3  07 

So  306 

0o  089 

0.  236 

0o  1  83 

4.348 

1  5.9 

3  -18 

1,0? 

0.076 

3.7,  0  67 

0,  2^9 

0  0  0  7  8 

0.020 

0o  0?p 

4  7,  Q 

1-20 

lo  3  3 

Oo  06  0 

9o  966 

0.239 

Go  131 

0o057 

1.234 

4  2,6 

1-21 

lo  50 

Oo  047 

3  3.686 

0,  205 

0o  1  23 

0.054 

1  0  1 1 1 

36.5 

3  -22 

1  c  6  9 

0  o  0  6 1 

10.163 

0  0  3. 9 1 

0.  3  66 

0.091 

1.786 

3  4,  0 

1  _  7  ** 

l  Po  3  5 

Oo  003 

Oo  ooo 

0a  000 

0.000 

0,000 

0o  000 

Co0 

3-24 

2o  i  7 

Oo  0q7 

10. S  09 

0o  151 

C.  201 

0.  133 

7.  207 

2  6.8 

1-25 

23 

0o04^ 

1  2  ©  8  9  7 

Oo  103 

Go  ?49 

0o  204 

7,  3  **  3 

3  8,3 

1-26 

3.  16 

Go  055 

11. 295 

Go  1  06 

0.279 

0.257 

8.  037 

19.0 

1-27 

lo1G 

0  o  0  ?  0 

13.  730 

0.  33  2 

0.  063 

0.013 

0.  0Q-) 

55,7 

3-28 

1.31 

0o  0  28 

3 4. 468 

0o  362 

0.090 

0.027 

0.032 

6  40  5 

3-20 

1  o  6  4 

Oo  043 

1 1 o  ?4? 

Go  34  9 

0„  3  46 

0o  070 

0.691 

6  2  n  2 

1  -30 

3  0^4 

Oo  073 

8,  076 

0.457 

Oo  167 

0.092 

0.460 

83.5 

1  -71 

lo  P-5 

Oo  0  5  2 

1 0o 168 

0.332 

Oo  179 

Go  3  06 

1,127 

5  9.  3 

1-32 

lo  30 

0.051 

10.  3  10 

0o  334 

0o  175 

0o?  01 

I.  81 1 

5C,5 

i  -37 

2.0? 

Oo  059 

°o  606 

0.  309 

0.  210 

0o  3  46 

2.  485 

r  5,  7 

1-34 

2.4  2 

Oo  040 

12. 373 

0o  ?4» 

0o  3  96 

0.  3  26 

3. OPS 

44,2 

1-35 

2  ©  44 

Oo  044 

1 1. 6  75 

0o  251 

0.  2C9 

0o  144 

3.  3. 3  4 

44,3 

1-36 

2©  27 

Oo  0 7 2 

8.735 

0.786 

0o  259 

0.223 

2.793 

5  lo  0 

3-37 

2 o  25 

Oo  074 

8o  5  36 

0,  289 

0o  26? 

0.226 

2.767 

51.5 

3-38 

2  ©69 

Go  052 

10.  703 

0.  276 

0.251 

0o  2  OP 

3.  297 

4  2,1 

1-39 

3  o  1  1 

0.-0  51 

11.053 

0.  2  07 

0o  281 

0.261 

4.457 

3  6.9 

1-40 

30  29 

Oo  044 

1 2o 051 

0o  197 

0.273 

0o  2  46 

5.454 

?4 ,4 

1  -41 

a  Q3 

-  o  -  — ■ 

0o  0  74 

1 4, 279 

0.160 

0.  ?76 

0o  2  50 

7,  451 

2  8.  4 

1-42 

*3  p  3 

-  © 

0.0  35 

1 3. 995 

0,  3.62 

0 . 2 7T 

0,2  54 

7.  746 

78,3 

1-4-3 

2  o  94 

0.033 

1 3  a  2  3  3 

0,294 

0.22? 

0o  1  67 

p.  06? 

5  2.4 

1-44 

2©  7  6 

0.04  5 

11.000 

0  0  9  7 1 

0,251 

0.  207 

7.  5C6 

59.0 

1-45 

2.  3  0 

0.043 

13.  057 

0o  3  20 

0o  3  61 

0.085 

4,3  57 

2  3,4 

1-46 

2  a  1  7 

0.045 

12.774 

0o  13  9 

Oo  ]  67 

0.092 

3.  71  7 

?l.i 

1-4*7 

2,  34 

0.048 

12.403 

0,  091 

0.229 

0.177 

6.433 

16,3 

3  -48 

2o  3 3 

0o  0  50 

1 1,466 

0.  23  3 

0o  2C3 

0.  1  36 

3.409 

38.0 

}  —  4° 

2©  33 

0.0  55 

Uo  097 

0o  1  PI 

0o  2  53 

0.211 

4.  576 

3  2 . 0 

ro  iv>  S\  \.k>  \J\  ivj  >.n  JO  O  O  Jl  O'  DO  -f)  O'  ^  O  O'  h  O  v>J 
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ANALYSIS  GT  COARSE  BED-MATERIAL  FLUME  DATA 
GILBERT  -  GRADE  E,  1.7].  MV,  UNIFORM 


T  F  S  T 

V  V  E  A  *vj 

F 

C 

PE 

V*B 

Y 

PARTI 

P  8/n50 

1-EO 

2©  64 

Go  066 

°,Q4T 

Go  lc5 

Oo  266 

Oo  231 

3®  991 

34,  7 

1  _5i 

2  o  6  ^ 

0,047 

1  1  O 41 7 

0,281 

0®  2 2 T 

0®  177 

3.  426 

5  Co  2 

1-82 

2c  67 

Go  04  8 

I  i.o?°! 

0o  276 

0®  22 5 

Oo  1  8  2 

3,314 

4  9  o  7 

AM4LYS I S 

OP  C0ARS 

E  REP- 

MATERIAL 

FLUME 

DATA 

CILBFRT  -  G p  A D E  F, 

3  0  1 7  M  8  t 

UNIFORM 

T  R  ST 

V  M  T  A  \j 

c 

c 

R  3 

Y*B 

Y 

PARTI 

R 8/ 050 

3-  1 

2  o  3  4 

0.053 

11,230 

0,  1  04 

0.  203 

0,  077 

2,  100 

1  0,  0 

1  -  ? 

2  o  3  4 

0o  063 

3,  673 

0,13  6 

0,221 

0  0  0  8  7 

1,836 

11.1 

1-  3 

2©  ^0 

Do  0  69 

9.  842 

Oo  093 

0.  275 

0,134 

4,587 

9,0 

1-  4 

?o  60 

0,  076 

°o  288 

0,093 

0.  2  30 

0,  1  39 

4.  ??6 

7,4 

3.  -  5 

?0  66 

0  0  0  5  2 

10,531 

Oo  3.75 

0,251 

0,3  12 

2,837 

1 6 , 3 

1-  6 

2  o  70 

0,  0  50 

1  0,  79 1 

0,  3*7? 

0,  250 

Oo  3  3  3 

2,  874 

16,5 

1  -  7 

3o  3.3 

0o057 

10,235 

Oo  153 

0,  3C5 

0,  3.66 

40  5  45 

14.  7 

>— • 
l 

JO 

3o  24 

0o  054 

10,  6  03 

0.  3.48 

0,305 

Oo  1  6  c 

4 q  3  04 

1 4  „  ? 

1-  7 

3o  21 

0o  041 

3.1©  3.72 

0,  265 

0,233 

0  0  147 

2,  603 

2  5.  5 

1-30 

3 o  27 

0o  040 

11, 2° 3 

0,261 

0,  290 

Co  .1  47 

2,  5n4 

2  5o  0 

1-11 

?o  02 

0©  065 

1.0,  53  2 

0,  C85 

0,  1  7? 

0,066 

1,087 

80  1 

1  - 1 2 

2n7? 

0o  070 

9,  554 

0,074 

0,  244 

Oo  106 

3.  243. 

7©  2 

1-33 

2  o  27 

0o0S7 

9,116 

0,077 

0,2  50 

0,  1 11 

3,332 

7,4  . 

3  -14 

2o  ?7 

0o0^2 

3  1 0  A5 1 

0,  144 

0,  198 

0,0  70 

1  0  4  2  9 

1  3  0  8 

1-15 

2c  2  5 

0  ©  0  5  6 

3.0,  Q?3 

Go  145 

0, 2C6 

Oo  076 

1  0  2  7 1 

1  40  0 

3  -16 

2o  34 

Co  0  54 

3  3  0  ?  5  0 

0,1  3  7 

0,  250 

0,1 1.1 

3,577 

3  1,3 

1  —  17 

20  79 

0o  058 

1  0©  0 1  3 

Oo  120 

0,256 

0,1 16 

3o  46 3 

1 1 , 5 

1-3  8 

2 ©  5° 

Co  0  70 

9,  784 

0,  130 

0,  26  5 

0,3  25 

3o  267 

12,  5 

1-19 

2,  74 

0o  046 

3  3  ©  507 

0,229 

0,233 

0  0  1  0 1 

2,  157 

2  2.0 

1-20 

2  ->  77 

0©  04  5 

?.l,77Q 

Oo  225 

0S  236 

0,  070 

2,  204 

21.  7 

1-21 

3©  3  8 

0,057 

11 ©034 

Oo  1°2 

0,307 

Oo  167 

3,981 

18.5 

1-22 

3©  4  .9 

0c  0^0 

33,215 

0,137 

0©  310 

0.  3  7i 

4,  3  60 

3  8,0 

3-23 

3 0  3° 

0,036 

1 2, 975 

0,26  5 

0,261 

0,  121 

3,  35? 

2  5,  5 

3  -24 

3 ©  26 

0o0i  2 

3.1  ©73  7 

0. 2«3 

0 

0 

‘Y> 

-J 

CO 

0,i  37 

2,87]. 

27,  ? 

1-25 

2.  37 

0,052 

3  lo  783 

Oo  108 

0,2  01 

Oo  072 

2 .219 

10.4 

1-26 

2  0  *  ? 

Go  057 

11,222 

0  0  1 1 3 

0.203 

0,077 

2,  013  • 

1  0,  7 

1-27 

2o  41 

0,  Ohr 

3,  937 

0,  3  00 

0,263  . 

0,1.28 

7,07). 

10,  5 

3-23 

?o  55 

0.  076 

9,  716 

Oo  103 

Oo  267 

Go  1  22 

3 0  396 

7 

1  -29 

2  0  6  2 

0o050 

13 ,609 

0  ©  1  9 1 

0,2  26 

Oo  0 9 1 

U  95  7 

1  8,4 

3-30 

?o  66 

0,  049 

l.io720 

0,188 

0,  227 

0,0  92 

1,960 

18.1 

1-31 

2  ©5° 

0,0  54 

11,136 

0,  195 

0©23? 

0,  096 

1.  831. 

3.  8.  7 

3-3? 

?o  16 

0,053 

1  1,330 

0,  3.62 

0,273 

0,1  37 

4,185 

15.6 

1-33 

3©  0^ 

0.  060 

10,  609 

Oo  167 

Oo  291 

0,1  51 

3.664 

16.0 

3  -34 

3,03 

0,  0?,Q 

1 2, °74 

Oo  2?  7 

Oo  23  8 

Oo  100 

3,  46? 

2  2,8 

1-35 

2o  94 

0,050 

U,  216 

0,255 

0,262 

0,3  22 

2,569 

2  4  0  5 

,  •' 

' 
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ANALYSIS  nc  GOARSF  BFn-MATF  R I AL  FLUMF  DATA 
0IL8FPT  -  GD  AOF  Gt  4 0  9 A  M M  ,  UNIFORM 


T  F  ST 

V^EAN 

r 

C 

PD 

V*8 

Y 

PAPT1 

R8/D50 

1-  1 

2«  7 3 

0©  046 

11.452. 

Oo  1  65 

0.243 

0.067 

1.26  6 

10.2 

1-  2 

2o 

0.  0C0 

10.  FI  6 

0,  1  c8 

0.  271 

0.084 

2.14  0 

0  ,  8 

}  -  ? 

2o  34 

0„oa6 

10. 2  OR 

0 .16  5 

0.  231 

Go  0  FQ 

2.086 

10,  2 

1  -  4 

3  4 

0„  054 

10o  5^0 

0.157 

0.29  3 

0.3  01 

3.  184 

9,4 

i  -  q 

3o  1 4 

Go  0C7 

10.  ?  9  9 

0.152 

0.  3  05 

0.  1  06 

3.  030 

9,4 

3-  6  ♦ 

3.44 

0o  060 

10.069 

0.  14? 

0.  341 

0.  133 

4.  334 

8,  7 

1-  7 

3,  48 

Oo  061 

9.  °  76 

0.  140 

0.  34  9 

0o  1  39 

4.1  7  0 

8,7 

1-  8 

?  o  9  3 

Oo  034 

12.311 

0.  268 

0.2  47 

0o  067 

1.122 

1  6.  5 

1-  9 

3.06 

Oo  038 

12.617 

0.260 

0.2  42 

0.067 

1.  274 

16.  1 

3-IO 

3  o  2  5 

Oo  039 

1 1, 55  3 

0©  750 

0o  2  31 

0.090 

1.976 

1  6.0 

1-11 

7  o  43 

0.04? 

1  1  .  T  ?  3 

C.  251 

0.  31  o 

0.110 

2.  820 

1  5.  5 

1-1? 

3  O  4  ? 

Oo  04? 

11.015 

0.254 

0.31  2 

0J  1  1 

2.783 

1  5,7 

1-13 

3.  3  3 

Oo  047 

10,  5  42 

0.  237 

0.  3  63 

0.149 

?o  836 

1  4,  6 

1-3  4 

4  o  ?  A 

Oo  0  50 

10.326 

0.218 

0.413 

0.  1  94 

5.  043 

13.4 

1-^5 

4o  21 

Oo  0  5  0 

in.  33  3 

0.  220 

0.403 

0.  1  90 

5.446 

13.6 

3-16 

3,04 

Oo  0^6 

1 0o  °7  8 

0.384 

0.  277 

0.  0  87 

0,  QQ4 

23,  7 

3.-1  7 

3.6  9 

0.076 

11.434 

0,  330 

0.  322 

0.1  18 

2.  5 77 

2  0,  3 

l-iq 

%  74 

0.  0^5 

1 1  ,  5  3  6 

0,  7?4 

0o  3  26 

0.1  21 

2.580 

7  0.0 

1  -3,9 

4.10 

Go  038 

11.135 

0.  309 

0.  366 

0.153 

3o  405 

19.  1 

3.-20 

4  e  1  ? 

Go  0  37 

i 1 ,  493 

0.303 

0.359 

0.147 

3.930 

l  8,  7 

1-21 

4.5  2 

0.043 

10,  520 

0.  293 

0.431 
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T  A  B  L  *  A 
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ANALYSTS  Oc  CHARSF  BFP-MA  TER  I  Al.  FLUME  DATA 
MFYFR— PETPRRMULLCR-  4.40MM»  MIXTURE 


T  F  ST 

V  M  F  A  N 

F 

c 

•  RB 

Y*R 

Y 

PART? 

RR/050 

3-  3 

3o  03 

0.060 

11.171 

Oo  3  77 

Oo  271 

Oo  09  3 

1.958 

1  3.9 

3-  2 

?o  3C 

0o  0 F9 

n,  ?  4  ? 

Oo  344. 

0.293 

Oo  1  13 

2  o  7  3  3 

23.  6 

3-  3 

3o  40 

0C  OF 7 

1 1. ^65 

Oo  3  3  8 

Oo  297 

0.112 

2.679 

2\2 

3-  4 

3o  «7 

0o  044 

1 3o 0P1 

Oo  3/*^ 

Oo  29ft 

Oo  1  11 

3.  34  2 

23.0 

■*-  5 

3  o  3  ^ 

OoOFQ 

11  7  *}  o 

0  o  3 1  6 

0  o  ?  P  3 

0.  1  06 

2.  36  5 

2  ^0° 

3-  6 

2o  6  3 

0.076 

9C  QftO 

Go  275 

Oo  ?ft  9 

0.09], 

lo  53  0 

3  8.6 

3-  7 

2  o  41 

Oo  075 

lOo  09.1 

Oo  219 

0.239 

Oo  070 

Oo  943 

14.  5 

* 

t  \m  i=  a 


-  5 

ANALYSIS  OF  COARSE  BFD- MATERIAL  FLUME  DATA 
MFYER-'°FTCR& MULLER-  3.30MM,  MIXTURE 


t  r  st 

V  M  E  A  N 

F 

C 

PR 

Y*P 

Y 

PARTI 

R  8/6* 

A—  ] 

2,  22 

Oo  05], 

1 Oo 151 

Oo  557 

Oo  2  1  3 

0,081 

Oo  1  94 

5  1 ,0 

4-  7 

2„  1  5 

0e  052 

10, 146 

0,  51  7 

0,  ?  1  2 

0,  077 

0,  1  6  8 

47,  8 

4-  7 

2  o  07 

0  o  0  5  5 

10,020 

0,484 

0.207 

0,073 

0.  U  9 

44,  7 

4-  4 

0,  63 

0,  004 

1,776 

i  o  446 

0.  ?  57 

Oo  2  1  8 

0  o  006 

13?  ,9 

4-  5 

1  e  30 

0o  0  60 

Q0 

0,  371 

0,  181 

0.  0  56 

Oo  01  2 

3  4.  2 

4-  6  . 

Oo  0  7? 

6,095 

0  ,  ?  ?  6 

Go  1  79 

0,056 

0,003 

31.0 

4-  7 

1,^5 

Oo  075 

Oo  ?41 

0,  278 

Oo  157 

0,  042 

Oo  0 

2  6  0  7 

cc 

1 

>r 

lo  2  3 

0  o  0  3  3 

8,  86  5 

0  O  2?  7 

Oo  145 

0,036 

Oo  0 

21,6 

4-  o 

2.  14 

Oo  0  5  3 

10, 076 

0,516 

0,21.4 

0,078 

0.145 

47.6 

4- IQ 

?.o  4  3 

Oo  047 

1  0,  2  0° 

Oo  654 

0„  240 

0,098 

0,885 

60,  4 

4-n 

?o  16 

Oo  05? 

TO, 003 

Oo  51  6 

0,214 

0,078 

Oo  31  1 

4  70  7 

4-1? 

2.46 

Oo  048 

10,207 

0,662 

0.  241 

0,069 

1.078 

6  0.2 

4-13 

2  o  4? 

Oo  04  3 

10, 203 

0, 624 

0.  2  3  8 

0,  0  97 

1,055 

57,6 

4-14 

2 o  32 

Oo  05 0 

1  0, 2  3  ? 

Oo  579 

0.  27  7 

Oo  0  8  3 

Oo  797 

53,5 

4-15 

?„  7 3 

0,052 

10, 078 

0,  5^2 

0.  221 

0,0  P4 

Oo  5  9? 

51.0 

4-16 

?,  1  5 

0,  053 

10, i 09 

0,517 

0,  213 

0,  077 

0.  496 

47,  7 

4-  I"7 

2, 01 

Oo  0  56 

6,607 

0,463 

0,201 

0,069 

0.234 

4?  „  3 

> 

l 

JO 

1  O  P  7 

Oo  060 

6,  881 

0,  406 

Oo  1  89 

0,061 

0,119 

37,5 

4-1  9 

2o  16 

0,062 

10,234 

C  o  5 1 3 

0,211 

0  o  0  7  6 

0.  450 

47,4 

4—70 

Ioq7 

0,  0  56 

10,  2  30 

0,402 

0,183 

0,074 

0.  382 

4  5,4 

i— £ 

(NJ 

1 

>4' 

2  ,  1  7 

Oo  05? 

10,  19  5 

Oo  522 

Go  21  3 

0  o  0  64 

0,  164 

3  9,  5 

4-2? 

.  lcq? 

0  o  0  5  6 

10, 1 13 

0,421 

0,  lc9 

0,0  50 

Oo  017 

?  1 , 6 

4-2? 

2o  4  9 

0„  047 

10,  2  3 2 

0,6^6 

0,  ? 44 

0  o  0  8  3 

0,618 

46  0  6 

4-?4 

2o  01 

0  o  0  5  5 

10,041 

Oo  464 

Oo  200 

0,056 

Oo  0  30 

3  5.i 

* 


■ 

v.n  vn 
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ANALYST  S  OF  CQAPSF  0-  ma  T  p  d  I  \  L  FLUME  OATA 

MFY  FR-pFTFP ^MULLER-  2.00MM,  M I  XTIJP  F 


TFST 

V.MF  4N 

F 

C 

•  P  8 

v*e 

Y 

PART! 

a B/050 

5-  1 

2.20 

0o  042 

13.353 

0.351 

0.  165 

0.090 

2.  119 

63.0 

5-  2 

2  o  l  0 

0o07F 

1.4-0  944 

0.  246 

0.  14\ 

0,  070 

0.  95c 

47,  2 

C_  "3 

2.09 

0  o  0  7  6 

1.4.616 

0.24P 

0.143 

0.070 

1.055 

4  6.4 

5“  4 

2.  04 

0o  0^ 

1 60 706 

0.  1  75 

0.122 

0.050 

0.  F47 

3  2.4 

5_  a 

?  .  0  8 

0  o  0  4  1 

i 3. 626 

0.  293 

0.  157 

0.  0  80 

1.437 

1  5.  4 

5-  A 

2.47 

Go  044 

12.945 

0.  466 

0.19 1 

0. 1  08 

2.459 

75,7 

5-  7 

2  a  40 

0o042 

1 30 277 

0.305 

0.181 

0.  099 

2.441 

6  3.1 

5-  8 

2.  8° 

0  o  0  47 

12 o  ?  46 

0.  531 

0.234 

0.157 

0.  249 

3  90  3 

5-  o 

2,  56 

0,  050 

3  lo 824 

0.632 

0.225 

0.143 

1.794 

96.3 

5-10 

2 o  66 

0©  04Q 

1 2o 00O 

0.  631 

0.221 

0.1  4  0 

2.1  20 

9  7.  7 

5-3  l 

2  o  6  7 

0.050 

11.310 

0.640 

0.222 

0.  3  40 

2.  0  59 

9  Pa  6 

5-1.2 

2*  7? 

0©  040 

1 2. 104 

0.41  5 

0.225 

0.3  40 

7.427 

93.5 

5-17 

2  ©  74 

0c  0  49 

1 2. 064 

0.613 

0.  22  7 

0©  3  41 

2.  1  90 

9  0.  7 

5-14 

2  o  5  8 

0.045 

51  2.724 

0.514 

0.203 

0.117 

2,484 

7  9,  7 

5-15 

2o  51 

0.044 

12.  379 

0.  527 

0.  195 

0.108 

7.060 

P0.fi 

5-1 A 

2©  14 

0.  046 

12.819 

0.  362 

0.  1.67 

0.085 

3.333 

S9«  Q 

5-  l7 

2©  71 

0,  044 

1 2.920 

0.422 

0.179 

0.03° 

1.726 

6  3.  1 

. 

■ 

’ 
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ANALYSTS  Op  CHAPSF  P^P-MATF R I AL  FLUME  DATA 
BOG  AR  DIE.  YEN-  1  0  «  0  M  m  t  UNIFHPM 


TEST 

VMR  AN 

c 

c. 

PR 

Y*B 

Y 

PARTI 

RB/D5 

1-  1. 

3o  l.a 

Go  000 

015 

Oo  166 

0.352 

Oo  0  70 

0.  034 

5.1 

1-  2 

2,7=; 

0J  4] 

7  o  126 

Oo  150 

0.  b  85 

0.  0  84 

0. 0°4 

c.  8 

1-  3 

2.3? 

Oo  3  3B 

7  o  1°  3 

0.194 

0.393 

0.087 

0.13  9 

5.9 

1-  4 

2.77 

Go  1  2  7 

7.  66  0 

0.172 

0.366 

Oo  076 

0. 0°° 

5.  3 

1-5 

2 o  73 

0o  097 

Bo  40? 

0o247 

0.322 

Go  0  58 

0.  04  5 

7.  5 

1-  6 

2  o  7  5 

Go  1  67 

6o  35  8 

0.  2  84 

0o4?3 

Ool  05 

0.077 

8,7 

1-  7 

7®  70 

Oo  1  66 

6  o  4  ?  4 

Oo  265 

0.421 

0,  1  00 

Oo  030 

Pol 

1-  P 

2  o  67 

Go  1  49 

60  746 

0  o  2  76 

Oo  3 °5 

0.0  88 

Oo  035 

P.4 

1  -  9 

2o  74 

G„  1  38 

7o  06  0 

Oo  265 

Oo  38  8 

0  o  0  P  5 

0.  030 

Pol 

1-10 

2  oft? 

0o  1  63 

6,434 

Oo  2C4 

Go  4C7 

Oo  093 

0.  027 

c,  0 

1-11 

7o  71 

Go  1  48 

6.776 

Oo  28? 

0.3  9Q 

0.0  90 

0.054 

8.6 

3-12 

2o  71. 

0o  1  64 

60  604 

Oo  295 

0.41  0 

0.095 

0.076 

CoO 

1-13 

2  o  B  1 

0  o  1  39 

60  974 

Oo  286 

Oo  402 

0.091 

0.  059 

8.  7 

1-14 

2o  75 

Go  1  28 

7,43  4 

Oo  22  5 

0.3  71, 

0  o  0  77 

0.039 

6,9 

1-3  5 

2.P7 

0o  1  28 

7 o  3«6 

Oo  24] 

Oo  3  83 

0.0  85 

0.  1  1  7 

7,4 

1-16 

3  o  01 

0.11° 

7o  74  0 

Oo  201 

0.375 

0.079 

Oo  154 

6.  1 

1-17 

30  4  4 

0o  049 

0o  36  9 

Oo  2 1  0 

0.367 

0.076 

0.5  30 

6.4 

1  - 1  P 

2o  5^ 

Oo  1  0° 

7  o  1  4  3 

0.  255 

0.  356 

0.  072 

0.  07  5 

7,  8 

1-19 

3.  25 

Oo  0  56 

10. 047 

0.274 

0.  3  24 

Oo  05° 

0.  014 

P.  4 

1-20 

3.35 

Oo  069 

Bo  43  6 

0.371 

Oo  3 97 

0.089 

0.031 

11.3 

,  ' 
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ANALYSIS  C04RSF  BED-MATFR I Al  FLUME  OATA 
POGARfVI  8YEN-  6.PVM,  UNIFORM 


TEST 

VMEAN 

F 

c 

R  3 

Y*R 

Y 

RAPT  1 

33/^50 

2-  1 

2  o  3  3 

0,  1  00 

8  ,  6  5? 

0*  162 

0,273 

0,  064- 

0,02  0 

7,3 

?-  2 

2 , 2  3 

0o0n« 

9,141 

0,  129 

0.  244 

0,04° 

0,  CC5 

5,  8 

2-  3 

2o  73 

0,  098 

55  3 

0,200 

0o  3  1 9 

0,084 

0,23? 

9.0 

?-  4 

?.o  54 

0,  09* 

p  o  3 1  2 

0,1*0 

0,  2PP 

0 o  069 

0,  048 

6.  7 

?-  5 

2o  70 

0,  0°8 

3,630 

Go  17] 

0,31? 

0,0  31 

0,  215 

7,  7 

?-  6 

?o  51 

O,  09° 

3,  71  -> 

0o  1  29 

0,2  83 

0,068 

0,045 

5,3 

2-  ~f 

c.  O  -  C. 

0,095 

3  »  34  3 

0,128 

0©  2  8  5 

0  o  0  6  7 

0,  040 

5,  7 

2-  ? 

2  o  °-5 

0o0°4 

Bo  ^54 

0,160 

0,322 

0,0  P6 

0,4^.0 

7,  2 

2-  0 

2  e  55 

0o0°l 

9.  079 

0,  124 

0,781 

0,066 

0, 04  ]. 

5.6 

2-1 0 

2 , 5  R 

0, 1  04 

Po  47Q 

0o  126 

0,3  05 

0,  077 

0  ,  2  c  3 

5.  7 

2-11 

2 c  36 

0,  1  09 

Po  2C  3 

0,113 

0,285 

0,067 

0,057 

5,1 

2-1? 

2,  51 

0o  3  01 

3  ,  6  4  5 

0,  11  5 

0,  2  90 

0,  07Q 

0,08° 

5,1 

2-13 

2  o  6  3 

0,09  3 

3,  570 

0,  1°6 

0,313 

0,0  81 

0,  064 

8,  8 

2-14 

2o  64 

0o  1  1  0 

7 „  091 

0,235 

0,331 

0,091 

0.210 

10.5 

2-15 

2*5  3 

0,  0°9 

8,  56  0 

0,  177 

0,  29  5 

0,0  72 

0.  024 

7,9 

2-14 

2,75 

0,  096 

P  ,  6  6  7 

0, 1  °6 

0o  3  1  7 

0 o  0P3 

0,  06R 

P,  « 

2-17 

3o  4  2 

0,  044 

11 , 748 

0,26  1 

0,201 

0,070 

0.  3.4Q 

10,8 

2- IP 

2, 0-7 

0  o  0  5  6 

10, 734 

0,  16  3 

0,  276 

0,  063 

0,  101 

7,  C’ 

. 

*• 

/ 
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ANALYSTS  OF  C9APSF  BFD-MATFRIAL  PLUME  DATA 
BOGARni&YEM-  1 5  o  0  M  M  t  UNIFORM 


T  PST 

VME  AN 

F 

C 

RB 

V*B 

Y 

PARTI 

R 8/950 

■*  _ 

1 

30  61 

0a  063 

8o  087 

0*  543 

0*  4 46 

0o  0  75 

0o  01  2 

lloO 

3- 

3 

°o74 

Oo  062 

8o  3  74 

Oo  3P5 

0„  421 

0o  067 

0.  0C5 

7o8 

o_ 

-j 

40  0’ 

0o  047 

10*382 

0,^74 

Oo  3  83 

0o  056 

Oo  0 3 1 

70  6 

-J— 

4 

4,  C? 

Co  040 

Uo  60S 

Oo  344 

Oo  3K? 

0,047 

0.072 

7.0 

3- 

5 

A  0  03 

0  o  0  5  3 

10o 4? 1 

Oo  2 90 

Oo  38  7 

Oo  0  56 

Oo  061 

6,9 

3  — 

6 

4o  2° 

Oo  0  56 

9,752 

0  ©  ?"?  4 

0*440 

0  o  0  73 

0,088 

6,6 

3- 

7 

3o  95 

Oo  055 

1  0o  2  3 .2 

Oo  247 

0.385 

0o  0 5 6 

Oo  065 

6.0 

3- 

p, 

4  o  4  5 

0  o  0  5  1 

10o 545 

Oo  263 

0.  4?5 

0o  068 

0.  41  8 

6,  5 

3- 

Q 

3c  94 

Oo  055 

10,  741. 

Oo  260 

0  *  3  8  5 

Oo  0  56 

0.  126 

5  o  3 

'i  — 

10 

70  73 

Oo  05  0 

3  Oo  ?  7  0 

Oo  21  2 

0o  365 

0o  050 

Oo  01  6 

4,3 

' 

T  A  ft  L  r  A 


-  6 

ANALYSIS  OF  COARSE  BED-MATERIAL  F  LUME  DATA 
USFWS  -  4, 1  mm,  MIXTURE 


TEST 

V  MR  AN 

c 

C 

Rp 

V*R 

Y 

°  AP  T 1 

PB/09Q 

9-  1 

lo  86 

0.  045 

1  ?o  65  3 

0o  22  3 

0.147 

0.030 

0.005 

16.6 

9-  2 

1  o  QQ 

Go  0  43 

l.?o  ?9  5 

Go  247 

0.154 

0. 033 

0.023 

1  °o  4 

o_  -> 

2o  04 

Go  045 

12, 45  ° 

0o  2Tfi 

0.164 

0.037 

0.022 

20,7 

9-  4 

M4 

0.  044 

l?o 2° ft 

0o  ?0Q 

0.  173 

0.041 

0.  04A 

2  3.0 

q_  c 

?o?? 

Co  045 

l 2o  2 69 

0.340 

0.  1  81 

0.045 

0.036 

?c,  3 

9-  6 

2.  15 

Go  0r  3 

1  Go  9  ft  6 

0 . 40 1 

0.197 

0.053 

0.  3  0  2 

2°.  ft 

9-  7 

1.85 

0.  0  56 

1 1  o  2  5  ft 

0.211 

0.  165 

0.037 

0.  027 

1  5.6 

9-  ft 

2.00 

Go  05  3 

1 lo  494 

0.235 

0.  174 

0.042 

0.  062 

1  7,5 

o-  9 

2.1? 

0o  091 

11.722 

0.  259 

0.  5  81 

0.045 

0.39R> 

18.9 

9-10 

?o  ?? 

0c  0  50 

1 1 o  7?  6 

0o  2  7ft 

0,  1  ft  9 

0.  0  49 

0.  510 

20.  7 

o-i  I 

2  .  2  5 

0o  052 

1 lo  460 

0  o  3  0 1 

0.  19  7 

0.05-2 

0.619 

2  2.4 

9-1?. 

1  o  03 

0  o  0  c»  6 

11.334 

0.200 

0.170 

0.040 

0.013 

14,3 

9-13 

2  o  0  7 

0  o  0  5  2. 

11.66° 

0.217 

0.177 

0.043 

0.  086 

16.1  • 

Q—  ].4 

2o  C ft 

Go  055 

1  lo  267 

0.237 

0.185 

0.047 

0,180 

1  7.6 

9-15 

?o  1  9 

0o  095 

llo 262 

0.262 

0.195 

0o  0  52 

0,4  25 

19,  5 

9-1  6 

2o  30 

0o  0  54 

11. '21  4 

0.290 

0.  2 C- 5 

0.058 

0.961 

ft1, 6 

f-H 

i 

0 

20  40 

0.  094 

11 o 0°9 

0.  323 

0.216 

0,  064 

0.741 

24,0 

9  - 1  ft 

2b  42 

Go  059 

1  Go  42 7 

0.  36ft 

0.  2  31 

0.  073 

0.  510 

2  7.  3 

■o  IP 


>  r* 
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ANALYSTS  Pc  CHAP  SC  BET-MATERIAL  FLUME  DATA 
LIU  6  CARTER-  4©?MM.  wiXTiJPF 


T  c  ST 

y  M  rr  a.  \j 

F 

c 

•  R3 

1-  1 

2©  84 

0o  0 1  6 

21  ©  005 

Oo  ?9P 

I-  2 

3o  03 

0©  019 

19,  67R 

Oo  7 2  R 

1-  3 

3©  12 

0,01  7 

20©67  7 

Oo  R2° 

1-  4 

7©  72 

0©  016 

22©  404 

0©  241 

1-  5 

7  ©  37 

0o018 

20©  82 7 

0©  24  1. 

1  -  6, 

2  0  88 

0  o  0  !  7 

21©  (V  7 

0,272 

1-  7 

2*37 

0©  016 

21,664 

0©  273 

I-  B 

3o  01 

0©  0  1 6 

22© 097 

0©  275 

1-  o 

Ml 

OoO1  3 

0©  000 

0©000 

1-10 

3  *  1  3 

0,016 

7 2©  01.3 

0,288 

? -11 

3  0  2  3 

OoO?) 

1  8  ©  3  6  7 

0,321 

1-1? 

3,05 

0©  01  4 

0©  0  00 

OoOOO 

1-13 

3  ©  4  ° 

0,  014 

0©  000 

0©  000 

1-14 

3,47 

0©  OOO 

n0ooo 

0,000 

1— 15 

3,  3  P 

0©  OOO 

0©  00  0 

Oo  000 

1-16 

3  ©  7  R 

0©  007 

0,000 

OoOOO 

1  -?  7 

\  50 

0o  007 

OoOOO 

OoOOO 

1-13 

3  o  6 1 

0©  003 

0©  000 

0,000 

1-19 

Bo  4 0 

0,  010 

0,000 

0,  000 

1-20 

3o  42 

0©  006 

Oo  000 

0,000 

1-21 

30  4R 

OoOll 

0©  000 

0©  0  00 

1  -22 

3,  30 

0o004 

0,000 

OoOOO 

1-23 

3,  54 

0©  003 

Oo  000 

OoOOO 

1-24 

70  5^ 

0©  011 

0,000 

OoOOO 

1  -25 

3*31 

0,006 

OoOOO 

0©  000 

1-26 

30  12 

0©  005 

0©  000 

0,  000 

1-27 

3,73 

0  ©003 

Oo  000 

0©  000 

1-28 

3,  TJ 

0,010 

0©  000 

0©  000 

1  -?9 

2,4  7 

0,024 

17©  P61 

0©  148 

1  -30 

2  ©  5  8 

0o  0  23 

1 p©  32  7 

0©  154 

1-31 

2  o  6  3 

0©  0  23 

18©  077 

0  ©  16  5 

1-32 

2  o  79 

0o021 

19©  2  7  8 

0©  17? 

1  -33 

2 ©  83 

0©0  19 

20©  23  5 

0©  18  0 

1-34 

3©  C 2 

0©  024 

1  7  0  446 

0©  1.  °5 

1-35 

7©  00 

0©  023 

i  7©  873 

0,205 

1-36 

3  a  0  3 

Oo  024 

17© 643 

.  0©  20  5 

1-37 

3©  03 

0©  024 

17©  6  37 

0.218 

1  -38 

3  ©  3  0 

0,022 

1 3© 475 

0©  220 

1-30 

3 o  53 

Oo  0 1 8 

20,41 2 

0©  722 

1-40 

3  ©40 

0©  07Q 

1°. 250 

0©  2?  R 

1-41 

3©  63 

Oo  02  0 

19,54? 

0,  238 

1  -42 

3,20 

OoO  24 

1  8,  1  3  3 

0©  QM 

1-43 

3  ">7 

0©018 

21©  324 

0©  097 

i  -44 

7©  44 

0,015 

OoOOO 

OoOOO 

1-45 

3©  46 

0,  O'!  7 

21© 674 

0©  106 

1-46 

3  O  40 

0,018 

70©  777 

0©  1 1  4 

1  -47 

3  ©  65 

0©01? 

0,000 

OoOOO 

1-43 

.  3.41 

0,010 

0©  000 

0©  000 

Y*B 

Y 

PART] 

R 0/050 

0,130 

OoO?? 

0,017 

21©  1 

0,154 

0©  031 

0.01  1 

23,2 

0,151 

0,030 

0.  01  3 

2  7©  3 

0,  171 

0.01° 

0 , 0 1  c 

1  7,1 

0,138 

0©  025 

0,  012 

1  8©  5 

0,  1  37 

Oo  0?  5 

0  ©  0 1  3 

1  9©  3 

0,  133 

0©  023 

0,014 

19.3 

0©  13  6 

0©  0  24 

0,014 

I  Co  R 

OoOOO 

OoOOO 

OoOOO 

0,0 

0©  144 

0©  027 

0,  06  6 

2  Oo  4 

0©  179 

0©  04? 

0©  04  6 

;  2.7 

0,000 

OoOOO 

OoOOO 

0,0 

0,  000 

0©  000 

OoOOO 

0©  0 

OoOOO 

OoOOO 

0©  000 

0©  0 

OoOOO 

OoOOO 

0.000 

0,0 

0,000 

0©  0  00 

0©  000 

0,0 

OoOOO 

OoOOO 

0,000 

0.0 

0.  000 

OoOOO 

OoOOO 

0.0 

0,000 

0,0  00 

0,  000 

Oo  0 

OoOOO 

OoOOO 

0.000 

0,0 

0,  000 

OoOOO 

0©  000 

OoO 

OoOOO 

OoOOO 

0©  000 

0,0 

0,000 

OoOOO 

0,000 

0,0 

0,  000 

0©  000 

Oo  coo 

0©  0 

0,000 

0,0  00 

Oo  000 

OoO 

0,  000 

OoOOO 

OoOOO 

OoO 

0©  000 

0©  000 

0©  000 

0©  0 

OoOOO 

0.000 

0,000 

0,0 

0©  133 

Oo  075 

Oo  02 5 

10.5 

0,  1  41 

0,0  26 

0,021 

10©  9 

0,146 

0©  0  28 

0  ©  036 

11,7 

0©  145 

0©  027 

0©  02  7 

1  ?.  2 

0.14? 

0  ©  0  27 

0  o  0  9  5 

12©  7 

0  ,  1  7  3 

0©  039 

0©  07o 

13.8 

0,  163 

0©  037 

0.181 

1  4  ©  5 

0,  174 

0,040 

0,216 

14©  5 

0,  17? 

0©  039 

0,2  38 

15.5 

0.  179 

0©  04? 

0.  62  9 

1  5.  6 

0©  173 

0,039 

0,826 

15,7 

0©  18? 

Oo  043 

0. 767 

16©  1 

0©  186 

0,045 

1  o  33  5 

16,8 

0,176 

0.041 

0,045 

6.9 

0.15  3 

0.037 

0©  182 

6©  9 

0,000 

0,0  00 

0,  000 

C.  0 

0,  1  6  0 

0 o  0  34 

0©  845 

7,5 

0,  16  4 

0©  035 

0,  800 

8©  1 

OoOOO 

0,000 

OoOOO 

0,0 

0©  000 

0.  0  00 

OoOOO 

0.0  • 
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F 

r 

RP 

Y*R 
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PARTI 

RB/050 

S— • 

I 

-IS 

O 

30  5R 

Go  01  3 

20.962 

0,  122 

Oo  171 

OoO  38 

0.9  88 

Po  7 

1  -SO 

3  o  TO 

OoOl  6 

2  2  •  ?  ?  2 

Oo  121 

Oo  1  69 

0.037 

0.  835 

8,6 

1  -51 

20  SO 

0,010 

0,000 

0.000 

Oo  000 

Oo  000 

0.  000 

0,0 

1-52 

7,  50 

0o034 

15,1^0 

Oo  094 

0  0  1 7  2 

Oo  0  39 

0.010 

(r  ,6 

1-50 

2  o  6  7 

Oo  037 

14. 350 

Oo  1.04 

Oo  1  83 

Oo  04^ 

0.044 

7.4 

1-54 

2  o  RO 

Go  0  31 

1.  5o  507 

0  0  10  5 

Oo  1  79 

0.047 

0.  C  87 

7,  5 

1-55 

2o  34 

Oo  0?? 

15.474. 

0.107 

0.1  84 

0,044 

Oo  214 

7  0  6 

1-56 

2o  5? 

0  o  0  ?.  2 

19.  01  0 

Oo  106 

Oo  149 

Oo  029 

Oo  310 

7,  5 

1-57 

2  a  Q  5 

Go  0  3  2 

15,352 

Oo  109 

0  0  1  8  5 

Co045 

0.722 

7,7 

1  -KR 

?o  90 

Go  051 

1  5o  t 0 1 

Oo  11  2 

Go  1  90 

0o047 

0.461 

PoO 

1-59 

?.  1.7 

Go  0  79 

16, 

0  0  1  2  2 

Oo  I  94 

Go  049 

Oo  773 

8.  6 

1-60 

M3 

0o030 

15,007 

Oo  138 

Oo  200 

0o057 

1.286 

1-61 

3o  53 

0o02.5 

17.  519 

Oo  148 

Oo  205 

0,055 

lo  352 

1  0.5 
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F 
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Y 
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P  8/ 05  0 

2-  1 

2o  3  2 

0o  03  2 

0,  000 

0,000 

0,000 

0,0  00 

0.000 

0,0 

2-  2 

2  o  4  2 

0,0:  4 

0,  000 

0,000 

0©  000 

0,  0  00 

0,  0  00 

0.  0 

2-  3 

2 ©  47 

0,034 

0,000 

OoOOO 

0,000 

0,000 

0.  oco 

0,0 

?-  4 

2o  5  5 

n0  014 

0,  000 

Do  000 

0,000 

OoOOO 

0,000 

0,0 

2-  5 

2«  72 

0,010 

Oo  000 

o,oco 

0,  000 

Oo  000 

Oo  coo 

0,0 

2-  6 

2 o  64 

0,01  3 

Go  090 

0,000 

0,000 

0,000 

0,  000 

OoO 

2-  7 

2  o  66 

0©017 

21,470 

0o  272 

0,  124 

Oo  027 

0,191 

2  5.5 

2-  3 

2  o  79 

0,01  5 

23,005 

0, 26.9 

0.121 

0,  0  26 

0,  109 

2  5,  3 

2-  o 

2,  34 

0,014 

23, 504 

0,269 

0,121 

0o025 

0,150 

2 5,  2 

2-10 

2c  97 

0,0  3  6 

21,793 

0,  28R 

0.  136 

Oo  0^2 

Oo  1  P  7 

27,0 

2-U 

2oQ4 

0  0  0 1  8 

19,987 

0,304 

0,  143 

0.0  38 

0,262 

2  8,  5 

2-12 

30  08 

0,017 

20o  674 

Go  306 

0©  149 

0,038 

0,  ?99 

78,7 

2-13 

3  c  14 

0,017 

2  1 0  06  Q 

0.,  30  6 

0,  149 

Oo  033 

0.  255 

2  8,  7  ■ 

2-14 

2c  32 

0,021 

18,91.7 

0, 1.87 

0,  1.2  3 

0,026 

0,026 

1  7,  6 

2-15 

2o  36 

0,  017 

21,397 

0,  188 

0,  1  10 

0,021 

0,  026 

17,7 

2-16 

2  o  5  7 

0  ©  0  1 7 

71  o  74  3 

0,  19P 

0,113 

0,  0?4 

0,  047 

1  P.  6 

2-17 

20  60 

0,017 

21,63 8 

0,204 

0.  1  20 

0,025 

0,084 

IP.  2 

2-18 

2,66 

0o  015 

0,  000 

0o  000 

0©  0  00 

0,  000 

OoOOO 
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2 -IP 

2  o  7  6 

0,  017 

2 1 0  6  2  9 

0,  220 

0,  123 

0,023 

0,292 

2  0,  7 

2-20 

2o  32 

0,01  3 
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0,  000 

0,000 

0,000 

OoOOO 

0,0 

2-21 

2  o  3  5 

Gobi  6 

22,  56  4 

0,725 

Oo  1  26 

0,02  3 

0.  624 

21.1 

2-22 

2  o  3  9 

0,015 

27,646 

0,  230 

0©  1  23 

0,02  3 

0,  67 C 

2  1.6 

2-23 

2o  36 

0,  01  8 

2.0,  5  00 

0,  242 

0,  139 

0,034 

0,628 

22,7 

2-24 

3o  00 

0o  0 16 

22,203 

0©  737 

0,  13  5 

Oo  032 

0,  775 

22.2 

2-25 

3  c  1  4 

0,015 

22,843 

0.744 

0,  137 

0o033 

1,032 

2  2,9 

2-26 

3o  1 6 

0,  01 .3 

o,  000 

0,  000 

0,  000 

0,0  00 

0,000 

0.0 

2-27 

3  o  2. 5 

0,012 

0,  000 

0,000 

0,  000 

0,  000 

0,000 

0.  0 

2-28 

3,  R0 

0.015 

9?o  506 

0,  267 

0,3  47 

0,037 

1,907 

2  5,  1 

2-29 

3o  64 

0,  01  3 

0o  000 

0,  000 

0,  000 

0,000 

Oo  000 

0,  0 

2-30 

2  o  5  3 

0,025 

17, 359 

0,  130 

0,  1  44 

Oo  0  36 

0,  047 

12,  2 

2-31 

2,  23 

0,  029 

16, 271 

0,  153 

0,140 

0,034 

0  0  0 1  3 

14,  3 

2-32 

2©  43 

0,023 

16, 430 

0, 1.47 

0o  1.46 

Oo  037 

Oo  012 

1  3,  S 

2-33 

2,6  8 

0,024 

17,934 

0,138 

0,  149 

Oo  039 

0,167 

13.0 

2-34 

7,63 

0,  026 

17*  6  3Q 

0,  145 

0,  153 

0,041 

0,27  8 

13.6 

2-35 

2.4  3 

0o  026 

1 60  94  3 

0©  1  co 

0©  146 

0o037 

0,  01  8 

14.9 

2-36 

2.72 

0,0  25 

1 7  e  6  3 1 

0,3  48 

0©  154 

0,041 

0,333 

13.9 

2-37 

2,  72 

0,02  5 

17,  36  5 

0,159 

0,  157 

0,  0  43- 

Oo  51  6 

14.9 

2-38  • 

7.80 

0,024 

1 7, 636 

0,163 

0,  1  59 

Go  044 

0.613 

1  5,  3 

2-30 

2,89 

0,024 

17, 972 

0,165 

0,161 

0,045 

0,672 

19,5 

2-40 

7,99 

0,023 

1  8,  05  0 

0,172 

0o  1  65 

0,  047 

0.944 

1  f  0 1 

2-41 

3  ©  09 

0,022 

1 3a4Q6 

0, 1.76 

0,  i  67 

Oo  049 

1.  17  4 

1  6,  5 

2-42 

3.  71 

0,  021 

1.9,  1  73 

0,179 

0o  167 

0,048 

1.41  4 

1  i  ,  3 

2-43 

3,35 

0,  0  20 

19,75  0 

0,187 

0o  3.  70 

Oo  050 

1.  48  4 

1  7,  5 

2-A.4 

3,43 

0,019 

30, I 40 

0,  3.9A 

0,  170 

0,0  50 

.1,  851 

18,  4 

2-45 

2,  n 

0,  040 

13, 847 

0,  077 

0o  I  52 

0,040 

0,025 

7  0  2 

2-46 

2,11 

0,043 

1 3, 334 

0,  0«3 

0,  153 

0,043 

0.  024 

7.  7 

2-47 

2,77 

0,03  6 

14,669 

0,082 

0,355 

0,042 

0,199 

7.7 

2-48 

O  9  0 

r.  O  /. 

0o  038 

14,161 

0,084 

Oo  1  57 

0,043 

0,232 

7,9 
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RB 
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Y 
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2-49 

2 a  70 

Oo  0  76 
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Oo  084 

Oo  15A 

0.042 

0.105 

7.9 

2-50 

2c  76 

Oo  077 

1 4„  474 

Oo  058 

0.  1  67 

0.  046 

0.  457 

P.2 

2-51 

?•  76 

0.040 

1 3o 506 

0.094 

Go  1  70 

0.0  50 

0.  37? 

8.  9 

2-52 

\46 

Ooms 

1  4.  160 

Go  C°8 

0  c  1  7  3 

0,052 

Oo  5° 2 

2-  ^ 

2o  56 

0o07? 

1  5  0  473 

0.095 

Oo  165 

Oo  047 

1.058 

Bo  9 

2-54. 

?  0  5  a 

Oo  032 

1  5  0  5  ?  7 

0.094 

0  c  1  6  3 

0.047 

0.  606 

Bo  9 

2-55 

2,  6i 

Oo  032 

1  5o  3  04 

Oo  0°9 

Oo  1  69 

0.0  50 

Oo  55° 

c  .3 

2-c6 

?n  6-? 

Oo  0  32 

1 5, 576 

0,  103 

Oo  163 

Oo  049 

1.059 

c«  6 

2-57 

2.  3  7 

Oo  0  27 

r,  07  7 

Oo  057 

Oc  1  66 

0.043 

1.461 

9.  1 

2-5R 

7.67 

0  0  0  7  3 

1  50  1  0  7 

Oc  106 

Oo  177 

Oo  0  54 

1.258 

°  0  9 

2-59 

2  0  F  7 

Oo  02° 

16. 153 

Oo  107 

Oo  173 

Oo  0  5  5 

I.  370 

10,0 

2-50 

?o  51 

0  *>  0  3 1 

15.  571 

Oo  1.1? 

Go  IPO 

0.05  6 

1.259 

10.  5 
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C 

RP 
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R R/050 

3-  1 
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0©  01  * 

.?  Oo  40  3 

0.  1.8  9 

0©  099 

0.  024 

0.  015 
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3-  2 

2.’0 

0©  021 

1  0 . 2  4  ° 

0.200 

0.109 

0.0  30 

0©  013 

2  7  0  0 

3-  3 

2o  1  9 

0.01S 

?0o  °3  7 

Oo  ?05 

0©  1  04 

0,037 

0.  046 

27.7 

3-  4 

2,  l  6 

0  o  0 1  3 

?n, 614 

0©  2.30 

0©  1  C5 

0.027 

0.  1  80 

2C,7 

3-  5 

2.34 

0©  0?0 

19. 35  0 

0.212 

0.13? 

0.037 

0.  1  28 

31.3 

3-  6  . 

2c  4  1 

0.019 

2  0.  007 

0©  22° 

0.126 

0.03° 

0.371 

^  0  -) 

/  0  £ 

3-  7 

2.50 

0©  022 

13©  514 

0.252 

0.13  5 

0.045 

0.  307 

3  4,  0 

3-  P 

3.53 

0.008 

0©  00  0 

0.000 

0.000 

O.O'OO 

0.000 

0.0 

3-  9 

1.33 

0o  090 

7  ©  7  5  0 

0.  567 

0.  1  79 

0.080 

0.  3A9 

7  6.4 

3-10 

2  o  5  5 

OoOlS 

?0, 525 

0®  273 

0.  124 

0  ©  0  3  8 

0©  7?  3 

36 . 9 

3-11 

2,77 

0©  017 

2 1  ©  43 1 

0.267 

0  ©  1 2  9 

0.042 

1.147 

35.9 

w  U 

2  o  S3 

0©  015 

22© S3  4 

0.  267 

0©  126 

0.040 

1.687 

36©  0 

3-]  3 

3.  02 

0©  017 

21  © 17? 

Oo  281 

0.143 

0.051 

1.607 

3  7,  9 

3-14 

oo 

0,  0.14 

2  3©  8 7 4 

Oo  271 

Oo  l  2  5 

0.03° 

2  ©98  7 

36.6 

3-15 

3.0? 

0o01  3 

20. 6  7 4 

0.294 

0©  146 

0.  05  3 

2.  51  2 

3  7 

3-16 

2, 01 

0.025 

17.446 

0.158 

0.115 

0.033 

0.015 

21,3 
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APPENDIX  5 


COMPUTATION  OF  BED- LOAD  TRANSPORT 
USING  BLENCH'S  REGIME  EQUATIONS 


NORTH  SASKATCHEWAN  RIVER 


AT 

DRAYTON  VALLEY 


Computatio n  of  Bed-Load  Discharge 

Using  Blench 's  Regime  Equations 

The  following  slope  equation  was  used: 

S  =  *  Fho11/12  f11X(c) 

K  bX/6  Q1/12 

where:  (for  cross-section  89) 
k  =  2.00  (transverse  bars) 

S  -  0.0015 

Q  =  35,000  -cf s 

K  =  vitI  =  1950 

b  =  775  ft. 

w 

F,  =  4.6  from  Dc A (number )  =  0.09  ft. 
bo  50 

or  Dj-q  (weight)  =  0.16  ft. 
and  FIGURE  7.3  (Blench  1969) 

f111  (c)  =  0.0015  x  1950  x  (775)  1//6  x  ( 35 , 000 )  1/12=  2.75 

2.00  x  (4. 6)  11//12 

From  FIGURE  7.2  (Blench  1969) 

C  =  17  parts  per  100,000  by  weight. 

Therefore 


CxQ 
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1600 
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=  0.48  lb/ft. sec. 
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RESISTANCE  TO  FLOW  DATA 
IMMOBILE  CHANNELS 
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APPENDIX  7 


RESISTANCE  TO  FLOW 

FLUME  TESTS  WITH  ARTIFICIAL  CEMENTED 
BED  AND  NATURAL  SORTED  BED 


DESCRIPTION  OF  FLUME 


The  one-foot  wide  wooden  tilting  flume  at  the 
University  of  Manitoba  Hydraulics  Laboratory  is  shown  in 
PHOTOGRAPH  A-l.  This  flume  is  37  feet  long  and  is  mounted 
on  a  truss,  whose  slope  can  be  adjusted  by  a  hand  cranked 
shaft.  Along  the  top  of  the  flume  are  mounted  steel  angles 
which  support  a  movable  point-gauge  carraige.  The  point- 
gauge  can  measure  vertical  elevations  relative  to  the  flume 
to  the  nearest  0.001  foot.  A  slotted  end-gate  is  used  to 
control  the  water  surface  profile  in  the  flume.  The 
water  supply  pipe  at  the  inlet  is  fitted  with  a  calibrated 
orifice  plate  for  discharge  measurement. 

PREPARATION  OF  THE  ARTIFICIAL  CEMENTED  BED 

The  material  used  in  the  flume  tests  consisted 
of  gravel  from  1/4  to  3/4  inches  in  size.  Sizes  smaller 
than  1/4  inches  were  excluded  from  the  mixture. 

The  mixture  was  cemented  to  plexiglass  sheets 
which  were  placed  in  a  wooden  form.  The  sheets  were 
coated  with  a  contact  cement  (3M  Ten  Bond)  and  a  four 
inch  layer  of  gravel  was  tamped  into  place  over  the 
sheets.  After  the  cement  hardened  (usually  24  hours)  the 
loose  stones  were  removed  and  the  plexiglass  sheets  were 
placed  into  the  flume.  In  order  to  prevent  movement  of 
the  sheets  small  holes  were  drilled  through  the  sheets 


and  these  were  nailed  to  the  wooden  bottom  of  the  flume. 


PREPARATION  OF  THE  NATURAL  GRAVEL  BED 

The  gravel  mixture  was  placed  in  the  flume  to  a 
depth  of  approximately  2  inches.  This  gravel  layer  was 
levelled  and  tamped  into  place.  Water,  having  a  fairly 
high  velocity,  was  then  introduced  into  the  flume  to 
allow  the  gravel  to  move  into  a  natural  position.  Some 
gravel  was  removed  from  the  flume  by  the  flow,  which  took 
place  for  2  hours,  but  at  no  location  was  the  thickness 
of  the  remaining  mixture  less  than  1  inch.  The  surface 
gravel  particles  would  then  be  sorted  by  the  flowing 
water,  in  fact,  many  of  the  particles  were  found  to  have 
their  long  axis  perpendicular  to  the  direction  of  flow  as 
in  actual  coarse-bed  rivers. 


' 


TABLE  A- 8 


RESISTANCE  TO  FLOW- ARTIFICAL  CEMENTED  BED 

Flume  Width  =  1.0  ft.,  'v)  =  1.41  x  10  ^ft/sec.^ 
Protrusion  Height  k  =  0.024  ft. 


Vm/ 

v 


TEST 

SLOPE 

Q 

d 

R 

Vm 

fb 

b 

Rb 

ft  3£ 

f  t . 

f  t . 

f  t  /  s 

ft. 

A- 1 

0 . 00334 

0 .290 

0.271 

0 .176 

1.07 

0 . 203 

6 . 26 

0.248 

A- 2 

0.00334 

0 . 190 

0.232 

0 . 158 

0.82 

0.279 

5 . 35 

0 . 218 

A-  3 

0 .00334 

0 .265 

0.249 

0.167 

1.06 

0 .171 

6 . 82 

0.227 

A- 4 

0 .00334 

0 . 152 

0 .203 

0.144 

0 . 75 

0 . 296 

5 .17 

0.194 

A-5 

0.00334 

0 .092 

0 . 165 

0 . 124 

0.56 

- 

- 

— 

B-l 

0.00244 

0 .260 

0 .285 

0.181 

0 .91 

0 .216 

6 .05 

0 .258 

B-2 

0.00244 

0.210 

0 . 250 

0.167 

0 . 84 

0 . 169 

6 . 84 

0.228 

B-3 

0.00244 

0 . 152 

0 .229 

0 . 148 

0.66 

0 .286 

5 .26 

0.222 

B-4 

0 . 00244 

0.171 

0 .222 

0 . 154 

0.  77 

0 .317 

5.01 

0.210 

B-5 

0 . 00244 

0.176 

0.233 

0 . 159 

0  .  76 

0 .239 

5  .  76 

0 .217 

B-6 

0 .00244 

0  .  152 

0.222 

0.154 

0 . 69 

0 . 280 

5'.  32 

0.210 

B-7 

0 .00244 

0 . 052 

0 . 148 

0.114 

0 . 35 

- 

— 

— 

C-l 

0 .00183 

0.272 

0 .286 

0 . 182 

0.95 

0 . 130 

7 . 80 

0.250 

C-2 

0 .00183 

0 .240 

0.280 

0.180 

0.  86 

0 . 160 

7 .07 

0 .250 

C-3 

0 .00183 

0  .  206 

0 .265 

0.173 

0.78 

0 .186 

6 .54 

0 . 240 

C-4 

0 . 00183 

0 .162 

0 . 241 

0.163 

0.67 

0 . 230 

5 . 87 

0 .222 

C-5 

0 .00183 

0 . 100 

0 . 199 

0 . 142 

0 .50 

- 

- 

- 

C-6 

0.00183 

0.058 

0 . 149 

0.115 

0.39 

— 

— 

— 

D-l 

0 . 00122 

0.276 

0.327 

0.198 

0 . 84 

0 . 123 

8 .02 

0.297 

D-2 

0 . 00122 

0 . 220 

0. 343 

0.203 

0.64 

0 . 236 

5 .78 

0 . 308 

D-3 

0 .00122 

0.270 

0 . 339 

0.202 

0 . 80 

0 . 145 

7 .38 

0.292 

D-4 

0 .00122 

0.225 

0. 315 

0.193 

0.71 

0 . 171 

6 . 83 

0.276 

D-5 

0 .00122 

0 . 191 

0.268 

0 .175 

0  .  71 

0 . 140 

7.51 

0.236 

D-6 

0 . 00122 

0 . 151 

0 .265 

0.173 

0.57 

0 .234 

5 .81 

0 .241 

D-7 

0 .00122 

0 . 135 

0.240 

0. 162 

0 .56 

0 .217 

6.06 

0 .218 

D-8 

0 .00122 

0 .058 

0 . 191 

0.138 

0 . 30 

— 

— 

— 

E-l 

0 . 00061 

0 .264 

0 . 350 

0 .206 

0  .  76 

0 .076 

10 . 03 

0.274 

E-2 

0 . 00061 

0.234 

0 . 356 

0 . 208 

0.66 

0 . 105 

8.62 

0.238 

E-3 

0 .00061 

0 . 192 

0. 341 

0 .202 

0 .56 

0 . 143 

7 . 40 

0.246 

E-4 

0 . 00061 

0 . 160 

0 . 306 

0.192 

0.47 

0 . 140 

7.51 

0.248 

E-3 

0 . 00061 

0 . 120 

0.271 

0.176 

0.44 

0 . 195 

6 . 39 

0.244 

E-6 

0 . 00061 

0 .062 

0 .219 

0.153 

0.28 

- 

- 

- 

E-7 

0.00061 

0.275 

0.378 

0.215 

0.73 

0 .089 

9 . 40 

0 . 300 

Rb/k 


11 . 8 
10 . 4 
10 . 8 
9 . 2 


12  .  3 
18. 8 
9  .  6 
10 . 0 

10 . 4 
10 . 0 


11.9 
11.9 
11 . 4 
10 . 6 


14 . 1 
14  .  7 
13.9 

13.2 

11 . 3 

11 . 5 

10 . 4 


13. 1 
11. 4 
11.  7 
11 . 8 

11 . 6 

14 . 3 


♦ 

• 

1 

TABLE  A-9 


RESISTANCE  TO  FLOW-NATURAL  SORTED  BED 

Flume  Width  «  1.0  ft.,  V  a  1.21  x  10  ^  ft2/sec. 
Protrusion  Height  k  =  0.030  ft. 


TEST 

SLOPE 

'  Q 

d 

R 

Vm 

fb 

v*b 

Rb 

ft3/s 

f  t . 

f  t . 

f  t  /  s 

f  t  . 

F-l 

0 .00334 

0 .260 

0.254 

0 .169 

1.02 

0 . 194 

6.40 

0 .234 

7 . 8 

F-2 

0 .00334 

0.245 

0.250 

0.167 

0 .98 

0.207 

6 . 25 

0.230 

7.  7 

F-3 

0.00334 

0 . 202 

0.240 

0 . 162 

0.84 

0 .272 

5 . 40 

0 .224 

7 . 5 

F-4 

0.00334 

0 . 158 

0.208 

0 . 147 

0 . 76 

0 . 293 

5 . 20 

0 . 196 

6 . 5 

F-5 

0 . 00334 

0 .178 

0 .219 

0 . 152 

0 . 81 

0.267 

5 .44 

0.180 

6 . 0 

F-6 

0.00334 

0 . 108 

0 . 183 

0.134 

0 .59 

- 

- 

— 

— 

G-l 

0 . 00244 

0.273 

0.273 

0 .176 

1.00 

0 . 154 

7 . 20 

0 .244 

8.  1 

G-2 

0.00244 

0 .264 

0 .269 

0 .175 

0 .95 

0 . 170 

6  .  81 

0 . 245 

8.2 

G-3 

0 . 00244 

0 .241 

0 .266 

0 .173 

0 .91 

0 .183 

6.60 

0 . 240 

8.0 

G-4 

0 . 00244 

0.282 

0 .254 

0 . 169 

0 . 80 

0 . 234 

5 . 83 

0.236 

7 . 9 

G-5 

0.00244 

0 . 181 

0 .237 

0 . 161 

0 . 76 

0 . 236 

5 . 80 

0.220 

7  .  3 

G-6 

0 . 00244 

0 . 159 

0.223 

0 . 154 

0 . 71 

0 . 255 

5.57 

0.208 

6 . 9 

G-  7 

0 .00244 

0  .  115 

0.189 

0 .137 

0.61 

0 . 299 

5 . 16 

0 .178 

6 . 0 

H-l 

0.00183 

0.272 

0.293 

0 .185 

0.93 

0 . 139 

7 .51 

0.266 

8.9 

H-2 

0.00183 

0 . 260 

0.288 

0.183 

0 .90 

0 . 148 

7 . 34 

0.258 

8 . 6 

H-3 

0 . 00183 

0.233 

0.283 

0 . 181 

0 . 82 

0 .176 

6  .  73 

0 .256 

8.6 

H-4 

0 . 00183 

0 . 192 

0 .262 

0 .172 

0.73 

0 . 209 

6 . 16 

0.239 

8.0 

H-5 

0 . 00183 

0 . 152 

0.233 

0.158 

0 .65 

0 . 236 

5 . 80 

0 .215 

7  .  2 

H-  6 

0 .00183 

0 . 122 

0.210 

0.148 

0 .58 

0 .275 

5 .38 

0.197 

6 . 6 

H-7 

0 . 00183 

0 .073 

0.193 

0 .139 

0 . 38 

— 

— 

— 

— 

1-1 

0 . 00122 

0.268 

0 . 329 

0.198 

0 . 82 

0 . 132 

7  .77 

0 .281 

9 . 4 

1-2 

0.00122 

0 .251 

0.  318 

0 . 194 

0 . 79 

0 . 138 

7 .56 

0 .276 

9  .  2 

1-3 

0 . 00122 

0 .213 

0 . 301 

0 .188 

0 . 71 

0.167 

6 .90 

0.268 

9 . 0 

1-4 

0.00122 

0.178 

0 . 260 

0.171 

0 .69 

0 . 154 

7 .17 

0.231 

7  .  7 

1-5 

0 . 00122 

0.178 

0.270 

0 .175 

0.66 

0 . 167 

6 .90 

0.233 

7  .  8 

1-6 

0.00122 

0.140 

0.253 

0.168 

0 . 55 

0 . 237 

5 . 80 

0 . 195 

6 . 5 

1-7 

0 . 00122 

0.128 

0. 238 

0 . 161 

0 .54 

0 .238 

5  .  80 

0.220 

7.  3 

J-l 

0 .00061 

0.278 

0  .  369 

0 .212 

0 . 76 

0.079 

10.00 

0.287 

9 . 6 

J-2 

0 . 00060 

0.259 

0.362 

0 .211 

0 . 72 

0 .089 

9 . 40 

0 .290 

9  .  7 

J-3 

0 . 00061 

0.214 

0 . 337 

0 . 201 

0 .64 

0 . 108 

8.57 

0.280 

9  .  3 

J-4 

0 . 00061 

0 . 199 

0. 330 

0.198 

0 .60 

0 . 118 

8 . 20 

0.272 

9 . 1 

J-5 

0.00061 

0. 178 

0. 324 

0.197 

0 .50 

0 . 145 

7 . 40 

0.272 

9 . 1 

J-6 

0 .00061 

0 . 150 

0.294 

0 .185 

0 .51 

0 . 154 

7 . 16 

0.268 

8.9 

J-7 

0 . 00061 

0.110 

0.240 

0 .162 

0 .46 

0 . 158 

7.07 

0 .213 

7 . 1 

. 

* 

. 

' 

TABLE  A-10 


RESISTANCE  TO  FLOW  NATURAL  SORTED  BED 
MANNING'S  n  COMPUTATION 

Flume  Width  =  1.0  ft.,  v  =  1.21  x  10~5  ft2/s 
Protrusion  Height  k  =  0.030  ft. 


TEST 

:  q 

Rb 

f  t 3  /s 

f  t . 

F- 1 

0 . 260 

0.234 

F-  2 

0.245 

0.230 

F-3 

0 .202 

0.224 

F-4 

0. 158 

0. 196 

F-  5 

0 . 178 

0 .180 

F-6 

0 . 108 

- 

G- 1 

0.273 

0.244 

G-2 

0 .264 

0.245 

G-3 

0.241 

0 . 240 

G-4 

0 . 202 

0.236 

G-5 

0.181 

0.220 

G-6 

0 . 159 

0.208 

G-7 

0. 115 

0.178 

H-l 

0.272 

0.266 

H-2 

0.260 

0.258 

H-3 

0.233 

0. 256 

H-4 

0. 192 

0.239 

H-5 

0 . 152 

0 .215 

H-6 

0.122 

0.197 

H-7 

0 .073 

— 

1-1 

0.268 

0.281 

1-2 

0 .251 

0.276 

1-3 

0.213 

0 . 268 

1-4 

0.178 

0.231 

1-5 

0.178 

0.233 

1-6 

0.140 

0 . 195 

1-7 

0.128 

0 .220 

J-l 

0.278 

0 .287 

J-2 

0 .259 

0 .290 

J-3 

0 .214 

0 .280 

J-4 

0.199 

0.272 

J-5 

0.178 

0.272 

J-6 

0.150 

0 .268 

J-7 

0 . 110 

0.213 

Rb/k 

n 

n/Rb1/6 

7 . 8 

0 .032 

0 .041 

7  .  7 

0 .033 

0 .043 

7 . 5 

0 .037 

0 .048 

6 . 5 

0 .037 

0 .050 

6 . 0 

0 .046 

0.062 

8.1 

0 .029 

0 . 036 

8.2 

0 .030 

0 .038 

8.0 

0 .031 

0  i  0 40 

7.9 

0 .035 

0.045 

7 . 3 

0 .035 

0.045 

6 . 9 

0 .036 

0 .047 

6 . 0 

0 .039 

0 . 051 

8.9 

0 .029 

0 . 035 

8 . 6 

0 .029 

0 .036 

8.6 

0.031 

0 .039 

8.6 

0 .033 

0 .042 

7.2 

0 . 035 

0 .045 

6 . 6 

0 .038 

0 .049 

9 . 4 

0 .028 

0 .034 

9.2 

0 .028 

0 .034 

9 . 0 

0 .031 

0 .037 

7.  7 

0 .029 

0 .036 

7 . 8 

0 .030 

0 .038 

6 . 5 

0 .032 

0.042 

7.3 

0 .035 

0.045 

9.6 

0 . 021 

0.025 

9  .  7 

0 .022 

0.027 

9.3 

0 .024 

0 .030 

9 . 1 

0.025 

0.032 

9 . 1 

0 .028 

0 .034 

8.9 

0 .030 

0.037 

7 . 1 

0.029 

0.037 

- 

* 

PHOTOGPAPH  A-l 

FLUME  USED  FOR  RESISTANCE  TO 
FLOW  STUDIES 


APPENDIX  8 


BANKFULL  STAGE  DATA 
COARSE-BED  CHANNELS 


. 


TABLE  A- 11 


BANKFULL  STAGE  DATA 
coarse-bed'  CHANNELS 

Q  Vrr,  s  d  ws  D50  DgQ 

River _ cfs  ft/sec  xlO  ft.  ft.  ft.  ft.  ft. 

Kellerhals 

(1963) 

O 

Chilko  R.  at 


Henry ' s  Crs . 

5600 

8.1 

5.03 

5.5 

150 

4.6 

.47 

.83 

Taseko  R. 
below  Taseko 
Lake 

6400 

.  8.0 

2.9 

6.2 

150 

5.3 

.49 

1.0 

Chilko  R.  at 
outlet  of 
Chilko  Lake 

5000 

3.4 

.92 

270 

5.5 

.37 

.58 

Cariboo  R. 
at  Quesnel 
Forks 

12000 

9.6 

4.19 

8.5 

190 

6.6 

.88 

1.5 

Quesnel  R. 
at  Lawless 
Creek 

14800 

11.4 

6.33 

8.0 

200 

6.5 

.71 

1.3 

Cariboo  R. 
at  outlet 
of  Cariboo 
Lake 

8000 

4.7 

1.92 

260 

6.6 

.64 

1.1 

Galay  (1967a) 

North  Sask. 
River  ^ 

80000 

7.0 

1.50 

800 

14.3 

.10 

Alta . W. Res . 
(1967) 

13 

Highwood  R. 

6000 

5 . 5 

1.55 

205 

5.3 

.18 

Prairie  R. 

2900 

6.5 

3.20 

100 

4.5 

.15 

James  R,. 

7000 

2.0 

4.00 

105 

3.3 

.12 

Sheep  R. 

6000 

6.3 

3.80 

210 

4.5 

.11 

• 

’ 

TABLE  A- 11 


BANKFULL  STAGE  DATA 
COARSE-BED  CHANNELS 


River 

Location 

Q 

cf  s 

ft/sec 

S-3 

xlO 

d 
ft . 

W 

ft? 

d 

ft* 

?§° 

San 

1 

1500 

5.88 

2.80 

4 . 87 

73 

3.50 

.  32 

.70 

Luis 

2 

668 

5.83 

3.76 

2.81 

55 

2.08 

.26 

.50 

Valley 

4 

768 

6.53 

3.59 

3.11 

48 

2.54 

.26 

.  53 

Canals 

5 

448 

5.82 

3.68 

2.50 

40 

1.92 

.18 

.28 

6 

159 

4.59 

2.95 

1.88 

21.7 

1.60 

.14 

.  20 

V  ijanc  — 

95. 

6  4.36 

2.90 

1.73 

15.9 

1.38 

.14 

.18 

&  Lai  lbUIl  q 

1953)  r 

46 

3.00 

3.16 

.96 

19.2 

.  80 

.  13 

.  26 

y 

— 

— 

— 

— 

— 

— 

.  13 

.  19 

X 

10 

16. 

6  2.90 

9 .65 

.60 

11.1 

.  515 

.21 

.40 

n 

203 

3.88 

2.35 

1.88 

32.3 

1.62 

.17 

.24 

12 

128 

4.00 

2 . 43 

1.77 

21.9 

1.46 

.12 

.17 

14 

110 

3.29 

1.36 

2.00 

21.4 

1.56 

.07 

.12 

15 

477 

4.84 

1.99 

3.05 

39.4 

2.50 

.17 

.28 

16 

- 

- 

- 

- 

- 

- 

.17 

.  32 

17 

531 

5.51 

2.74 

2.60 

41 

2.35 

.13 

.17 

18 

235 

3.80 

.  80 

2.94 

25 

2.48 

.073 

.10 

APPENDIX  9 


DEPTH  OF  SCOUR  COMPUTATIONS 
PRAIRIE  CREEK 


DEPTH  OF  SCOUR  COMPUTATIONS 


Prairie  Creek  Near  Rocky  Mountain  House 
(x-sec.  7+60  and  8+60  U/S) 


Properties  of  Flow  and  Channel 


=  2900  cfs 

bw  =  110  ft. 


r  =  240  ft. 

0  =  90° 

ds  =  9  ft. 

D  =  2.3  in.  (grid  sample  by  weight) 


Lacey  Equations  (1929) 

a*  =  o.47  (Qb/f) 1/3 


f 

a* 


=  /64D 
=  0.47 


=  12.2 


(Questionable  whether 
applied  to  coarse-bed 


(2900) !/3 


12 .2 


=  2.9  ft. 


this  can  be 
rivers) . 


ds  =  2.9  x  factor  =  2.9  x  2  =  6  ft. 


Blench  Equations  (see  p.  122,  1969) 

=  qf2/3 


fo 


f  o 


F,  1/3  '  Fbo  5,0 
bo 


=  ( 


2900/110 


) 


2/3 


=  5.1  ft. 


ds 


(5)  1/3 

5.1  x  factor  =  5.1  x  2  =  10  ft 


FIGURE  112  (free  bend) 


r/ 


bw0 


240 


=  1.4 


110  x  1.57 


From  Figure  112,  dsA  =  0.10 

Dw 

Therefore : 

ds  =  110  x  0.10  =  11  ft. 


^  _ 
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